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ABSTRACT 


Cellular spaces, regular networks of Similar 
interacting components, constitute an important class of 
models in many disciplines, from automata theory to 
biological systems. The primary goal of this research is to 
design a comprehensive interactive programming system as a 
research tool for the simulation of cellular spaces at the 
University of Alberta. An important subgoal of the research 
is an examination of the fundamental requirements of an 
interactive simulation package. 

The proposed simulation package consists of a 
Simulation language, especially designed for describing and 
programming cellular models, and an interactive run time 
System with which the user can carry out his simulation. 
Simulation programming techniques are discussed; in 
Particuiar, the implications of using an interactive 
approach to running a Simulation are examined and a set of 
guidelines for designing interactive simulation systems 
suggested. 

A general description of the nature and 
characteristics of cellular spaces is provided and some 
examples of the types of models for which the system is 


designed given. The basis and motivation for this research, 
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R. F. Brender's system for simulating cellular spaces, is 
described, and its value illustrated by examples of cellular 
research for which it has been used. A detailed description 
of the proposed cellular space simulation package is given. 
The proposed package is evaluated in terms of how well it 
meets the criteria developed for interactive simulation 


systems and how it compares with Brender's system. 
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CHAPTER I 


INTRODUCTION 


The primary goal of this research is to design a 
comprehensive interactive programming system for use in the 
Simulation of cellular spaces (variously called cellular 
structures, tesselation structures, iterative arrays, 
tesselation structures, and tesselation automata). The 
proposed system is meant as a research tool for use at the 
University of Alberta and therefore is dependent on the 
hardware facilities available to users at this installation. 
An important subgoal arising from the research is an 
examination of the fundamental requirements of an 
interactive simulation package. 

The field of research with which the proposed system 
is primarily concerned is cellular automata theory. This is 
a fast growing branch of computing science which deals with 
the study of changes in the state of systems composed of 
regular arrays of components of the same type. This area of 
research promises valuable applications to the design of 


almost every type of automaton, from parallel computer 
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systems to models for the growth of living organisms. The 
proposed system has potential applications in other fields 
as well. Cellular models can be used for solving partial 
differential equations by means of relaxation methods {4], 
and for studying the growth of patterns (44, 38]. 

The motivation and basis for this research was the 
work done by R. F. Brender on a system used for cellular 
research at the Logic of Computers Group at the University 
of Michigan (4, 15, 16]. He and his associates designed and 
implemented a simulation system with which to study cellular 
automata. Following his lead, work was undertaken to look at 
providing a programming package which could be used for the 
Same purposes at the University of Alberta. The proposed 
system has drawn on Brender‘s design for ideas regarding the 
types of facilities which should be included in the systen. 
The vast differences in the hardware available, however, 
required different approaches to some aspects of the design. 
In addition, a set of fundamental capabilities valuable for 
interactive simulation was proposed and these were 
incorporated into the design of the package. As a result, 
the design presented in this paper is a more powerful and 
generalized system for cellular space modelling than that 
proposed by Brender. For example, facilities were included 
in the package to support run time user/program interaction, 
and to allow multiple definition of the behavioral 


Characteristics of the model in order to permit run time 
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changes in model behavior. The package consists of a 
Simulation language, especially designed for describing and 
programming cellular models, and an interactive run time 
system with which the user can carry out his simulation. The 
implementation of the system has been started but, due to 
the large amount of programming effort involved, has been 
curtailed and left for furthur development by others. 

The next chapter briefly examines Simulation 
programming techniques and describes seme of the problems 
which must be handled when representing real world phenomena 
on a computer. The general approaches taken by simulation 
languages designed for modelling natural systems on the 
computer are indicated, and their common characteristics 
described to obtain a basic set of criteria on which to base 
the design of the proposed system. A look at the 
implications of introducing the interactive aspect of 
computing into simulation is undertaken, and the needs of 
the run time simulation environment examined. A set of basic 
requirements is developed to provide guidelines as to what 
facilities and capabilities should be provided by an 
interactive simulation package. Finally, the implications of 
the specialized nature of the particular application are 
discussed. 

The third chapter is concerned with providing a 
general background and understanding of the types of models 


for which the proposed system is designed. It is devoted to 
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describing the nature and characteristics of cellular spaces 
and to reviewing some of the major achievements and 
highlights in cellular automata research. The most 
Significant and active areas of cellular space modelling are 
indicated, and impcrtant work done in each is described to 
indicate the types of application for which the system is 
intended. Brender‘s cellular space simulation system is then 
described. Evidence cf the utility of such a system is given 
by describing two important biological models which were 
implemented and investigated on Brender's system. 

The fourth chapter presents the detailed description 
of the proposed interactive cellular space simulation 
package. First, s description of the hardware and software 
environment for which the package was designed is given; 
second, the cellular space simulation language is described; 
and third, the run time system with which the user controls 
the execution of his simulation is outlined. 

The final chapter is primarily concerned with a 
general evaluation of the proposed simulation package. A 
discussion of how well the proposed system meets the basic 
criteria of an interactive simulation system suggested in 
Chapter if § is presented. The major differences and 
improvements of the proposed system in relation to Brender's 
system are then explained. The chapter concludes with a 
brief discussion cf the current status of the implementation 


and suggestions for the direction of future efforts. 
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CHAPTER IL 


INTERACTIVE SIMULATION 


2-1 Simulation 


The essence of simulation is imitation or role- 
playing. Cne entity, the device perforfing the simulation, 
is made to assume the nature of ancther entity, the 
phenomenon being Simulated. Simulation differs from 
experimentation in that the phenomenon under study is 
usually not physically involved in the simulation. 

Fundamental to the meaning of simulation is the 
concept of a model. A model is the representation of some 
phenomenon; the representation depends upon the desired aims 
of the model builder and the characteristics of the 
phenomenon being modelled. Simulation is the dynamic 
representation of a model, the process of producing 
sequential states in a model. 

Simulation is an experimental technique involving 
the iterative development of models and theories. One of the 


major objectives of simulation in research is to develop a 
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model which will imitate known characteristics of some 
particular natural system. The researcher determines the 
adequacy of a simulation model by comparing information 
regarding the behavior derived from state histories produced 
by the model with corresponding information regarding the 
behavior derived from state histories produced by the real 
system. Once the desired behavior is obtained, the model is 
closely examined in order to determine any unspecified 
Characteristics and relate these back to the real system in 
hopes of discovering some new information about the 
phenomenon. The task of developing the model, generating the 
required behavior, and discovering new aspects of model 


behavior can be made much easier with the aid of a computer. 


221.1 Simulation Programming Techniques 


The use of computers for simulation in scientific 
research and development has been extensive. Programming 
languages of all types have been used in simulation 
routines. In the early sixties, general purpose, 
mathematically oriented compiler languages such as_ the 
various versions of FORTRAN and ALGOL were probably the most 
common ones used for writing simulation programs. Early 
simulation programs were hand-tailored to the problems they 
were designed to sclve and were often quite complex. Such 


specialized programs were usually not flexible enough to 
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aliow model modification with reasonable ease and speed. The 
growing importance of computer simulation and the inadequacy 
of the simulation programming technigues of the time caused 
Simulation programmers to abstract the essential features of 
their prcegrams and create programming aids specifically 
oriented to simulation. 

Several different methods were attempted to provide 
the extra facilities needed for simulation programming. 
Information on the different types of approaches to 
Simulatien programming and their relative advantages and 
disadvantages can be found in [31], [42], and [43]. The two 
Main approaches were: 

1) The subroutine package approach ~- This approach 
consists of providing helpful routines, written in a 
particular general purpose programming language, which can 
be called on by a simulation program written in that 
language. 

2) The simulation package approach - This approach 
consists of creating a special simulation programming 
language and a run time simulation systen. 

The simulation package approach was the more popular 
method, largely because it allowed a much fuller range of 
facilities. A simulation language is usually a better 
vehicle for expressing the concepts arising in a simulation 
study than is a general purpose ilanguage. It contains 


special purpose statements for describing the entities, 


ebiverd) 47 Earqetin Give ole gepsecs () Levies 
oe eee etg, Mol toboesu (ge? betaen, gestelicss? enrE 
oo ~“@edesotgie— mm. aeqrt ‘% ake 2 «a0 oo oto 
fos Qapereavic vetaeies Mast Fae Waprearqary wae 
aev. O47 193 Ome Qi +25 WL! } ca sane’ ~i sen ag us 
. farce oa tosougge a 
40034 3° #2A7 * owen 1 vcard sajtieadis aff fF ~ “. | 
oe Oe weretse <yaattvers ‘vteaesa 5a Trem M 
mum dokds .(CRevetes C18) ree] sane In2ieen-8 
eld af weesdae ¢stigoa solteiudYe. ¢ ve no bedi 
coenggs ail? = Hosopis sywineg wnireiels set (48 
Ce ee Lo. a ee Ss wis urs 2@. 
i : ématzy © mulselenen. ard? ase 2) One 
meiebuiel was. ats doee-vaqe Secrets casei aap 


relationships and procedures of simulation in a brief and 
direct way. Such facilities can save considerable 
programming time and greatly simplify the task of making 
program changes. The function of a simulation language is to 
provide an orderly approach to describing the model and to 
provide facilities for performing simulation experiments 
with the model. 

The simulation system, often considered as just part 
of the simulation language, is the run time environment of 
the simulation. It consists of an executive program which 
controls the routines present during the running of the 
Simulation. The advantage of a simulation system is that it 
provides facilities required as part of any simulation but 
not available with general high level programming languages. 
For instance, the system can handle storage aliocation, data 
generation, output of results, task scheduling, and 
simulated time, the time reference used within a simulation. 

The underlying reason that simulation packages are 
possible is that simulation models and routines arising in 
one simulation study bear certain similarities to the models 
and routines arising in another simulation study. A 
Simulation package exploits these similarities in order to 
furnish a user with facilities that can be applied in most 


Simulation studies. 
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2.1.2 General Modelling Considerations 


There are many difficult problems inherent in 
representing a real world phenomenon on a computer. This 
section examines scome aspects of creating computer based 
system descriptions [3]. 

A description of a natural system may be viewed as 
consisting of twe prime components: (1) a static 
description, and (2) a dynamic description. The first is 
concerned with existence, what components and structural 
relationships between components will be recognized in the 
model? The second is concerned with change, what state 
changes are possible in the model and what are the sequence 
relationships between these changes? 

Computer models of the real world must contend with 
some non-trivial problems regarding the program description 
of natural system characteristics. The major constraints 
imposed on the static description involve the fixed storage 
capacity and organization of the computer. Some troublesome 
aspects of the static model description include: 

1) Transient nature of components - Entities and 
processes are dynamic in the sense that they not only 
change, but aiso enter or leave the system. Systems 
must allocate storage during simulation to deal with 


transient elements within the model. 
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2) Uniqueness of components - All natural things 
are different to some degree; however, any attempt to 
deal systematically with reality requires that 
uniqueness be replaced by classification to avoid an 
unmanageable proiiferation of structures. Data must be 
defined to reflect the significant characteristics of 
each component. 

3) Complexity of organization - Most systems of 
interest are very compiex and are characterized by 
highly organized structures of components. Complex 
relationships must be grouped Or combined into 
hierarchies of organization, and handled by an 
equivalent complexity in data structure. The data must 
be defined to reflect the relationships between 
compcnents. 

4) Numerous components - The great number of 
entities in systems of interest demands the use of 
scaling and aggregation aS a means of representing 


these within acceptable limits. 


The constraints imposed by the computer on the 
dynamic description involve the serial nature of digital 
computation and the use of stored programs. Some troublesome 
aspects of the dynamic model description include: 

1) Uniqueness of history - History is a stream of 


unigue occurrances with no two identical events. 
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Changes must te classified to permit a given program to 
execute the same type of change many times; patterns of 
behavior substitute for the unigueness of the real 
world. 

2) Continuous change - Activity in a system 
involves continuous change, often taking place over 
extended periods of time. In the model such changes are 
handled by a series of discrete changes Which occur 
instantaneously with respect to simulated time. 

3) Interdependence of components - Reai world 
Oca ioe are interdependent in that each event is 
influenced by its predecessors and influences its 
successors. Extensive interdependence in the real world 
implies interaction between program segments in the 
computer model. Either the interaction is continuous 
and a change in any state variable can affect any other 
thereafter, or else the points of interaction must _ be 
specified. In the latter case, computation may proceed 
up to an interaction point, whereupon control may pass 
to another activity. 

4) Concurrent events - Real world activities occur 
in parallel under independent motive forces. Such 
parallel activities give rise to concurrent events 
which must ke correctly sequenced and performed in 


series while simulated time is held stationary. 
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2.2 Simulation Languages 


Simulation languages are generaily categorized as 
either continuous change or discrete change, depending on 
their applicability to the characteristics of the phenomenon 
to be modelled [33]. Continuous change languages are 
appropriate when the subject consists of a continuous flow 
of information or material counted in aggregate rather than 
by individual items. Discrete change languages are used ued 
individual items are studied in relationship to other items 


and to their environment. 
2.2.1 Continuous Simulation Languages 


The continuous change model is generally represented 
mathematically by differential or difference equations that 
describe the rates of changes of the variables over time. 
The model is formulated with certain constraints and 
equilibrium relationships. Problems of this type were first 
solved on analog computers; however, because of the 
limitations of analog computers as to problem size, 
accuracy, and availability, some analog users have turned to 
digital systems. On the digital system, the continuous 
Change model is approximated by choosing some delta step 
size, solving all the equations, stepping the time interval 


by delta, and solving the equations again. This is continued 
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for some definite period of time or until the solution is 
obtained. 

Up until 1967 there were over 23 separate digital 
computer languages developed to perform simulation of 
continuous system dynamics. Some of the more important ones 
include CCBLOC, DAS, DES-1, DSL/90, DYNASAR, DYSAR, MADDLOC, 
MUDAS 0 CHINIC, FeLACTOLUS,. “PARTNER,® and SCADS. A general 
feeling grew that some sort of control and direction was 
needed in the area of continuous simulation languages in 
order to ensure orderly development of the field. The result 
was the specification for CSSL (Continuous System Simulation 
Language), a generalized continuous simulation language 
{14}. The CSSL specification has become the standard of 
evaluation for continuous Simulation languages. 

Continuous change sytems, however, lie outside the 
immediate context of the present study. A survey and summary 
of continuous Simulation languages can be found in [5], [7], 


and [42]. 


2.2.2 Discrete Simulation Languages 


There are several ways in which to approach discrete 
system programming. A few languages, notably IBM's GPSS, 
have used a flowchart approach where flowchart symbols are 
the basic model descriptors. The vast majority of simulation 


languages, however, are statement languages where the models 
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are formulated aS sequences of program sentences or 
statements. 

Statement languages use three principal philosophies 
Or practical approaches to the conceptualization of 
Simulation models and the implementation of these models on 
a computer. These basic system dynamics orientations are 
known as the event, the activity, and the process approaches 
[e372 30,.5234% 

The event languages allow the dynamic creation of 
entities in a system and study the events that these 
entities engage in. Event-centered simulation systems 
segment simulation models into submodels called events. 
Events are programs which specify changes in the system 
Status when these events occur. For example, in a typical 
waiting-line model, there are arrival events, end-of-service 
events, and unserviced-customer-departure events. The 
programmer is responsible for scheduling future events. 
Thus, in a waiting-line model, the presence of a waiting 
customer together with a server-becomes-available event, 
calls for the status of the server to be set to engaged, the 
renoval of a customer from the wait queue, and the 
scheduling of a server-becomes-available event at some 
future time. Future times are usually obtained by random 
sampling from statistical distributions. Whenever an event 
is scheduled, it is put on a calender of future events that 


is maintained by the simulation system executive program. AS 
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event programs are completed this executive removes events 
from the calendar according to their scheduled times and 
executes the events by calling on their respective programs. 

Activity languages consider fixed amounts of 
equipment as basic elements of a system and study the cycle 
of activities in which these equipments engage. The 
activity-centered Simulation systems segment simulation 
models by creating submodels called activities that are 
programs which describe changes of state. They are different 
from events in that they need not represent instantaneous 
state changes such aS an arrival or a departure, but can 
represent time consuming activities like servicing a 
customer. Within activity programs, clocks are set that 
indicate the times at which the various elements involved in 
the activity will change state. Generally at the head of 
each activity program is a set of test conditions that must 
be true for the activity to take place. There is no explicit 
scheduling done in an activity-centered system; a simulation 
executive program scans all activities each time one is 
completed to determine what to do next. AS activities are 
processed, element clocks are set and simulated time 
progresses. 

In some models the timing and scheduling of 
activites is implicit. For instance, in a fixed time step 
environment all activities take place instantaneously and 


concurrently at each time step. The clock is reduced to a 


vad _ / — nin te 
je wont ‘wate ma illness 


" waneapory ret 
‘40 4«6danGaA eels tebinec> salpnapané “pane 


; nin 
ofoes 88) YHGE Daly Harare BX aroyyrs> prada Joon gies 


eét Spat ite aitedaavise  eeure igiAe af Bart ivrsas- 


;Snse »e on evry sivneni 


 eossoioues nhae fice ‘* 
RA 
w2m, SREP WOETEVESS tom 2fehowiity » eEtest> pt eleme 


; 4 : yor 
somiesteh ope yest “a Fe? hs ee | $ ze } pe 7 he donde DaAKIDO “a 


Tr te |e ee , aa it A suova  aots 
an | _ 
amd $. : aie or ae $* Z ipa wePO ea> O78 is 


: Se 
P ua oiv 3 om - P Ju j WYSU0 “y wat seats = a 


sete $2 ‘3 eyoei : ieos ¥3iNL°0% « Latte 3900889, 


o 


mi bevio PRAe ad WARod ; i? @2pfd@ 26 &e4S9 O42 osnott Be 


$e) base od? te gli ez" veanee aueed: dite: vateisos sa 
a 
ee ee es ; 04 @h1208% ‘rttvidsos: doae 
Siobigns On @€S &4 Onis ege7 eG ivites ofF 103 ovata 
god -idele ooinedaye buysrsap-ytivheon we ot ood pad int) 
. et $a% agit 69a ante re if g0eca BA TT | vere 


obhtiwnaos aa oa7*p ue : (one: oainw?sd of teTtet . 


) eke —— yA. 2s ose 3910 digvants 
a ; _ a : snotty : 2 
naw )s 
: - peiwhs 000 es at 
: fi ‘* P71 one eeat soe | os oy 
oe : 7 


7 


> 


16 


simple counter and scheduling is basically absent (although 
some program manipulation may be required for state changes 
from interdependent but concurrent activities). 

The most recent approach to simulation modelling 
involves the process concept, a synthesis of the ideas 
invoived in previous approaches. Basic to this approach are 
the foliowing concepts: 

Process description = A process description 
consists of the distinguishing structural characteristics 
(described by data declarations) and the sequence of state 
changes (described by program statements) which outline the 
behavior of some component of the systen. 

Frocess - A process is a single instance of a 
process description. 

Interaction point - An interaction point is a point 
within a process description at which control may transfer 
to another process. 

A process is a dynamic entity that has a specified 
behavior pattern over time. Thus a process is a series of 
activities joined together by either time-dependent or 
status-dependent conditions. An event is a point in time at 
which processes interact; hence it is a condition point 
separating activities within a process. 

There are any number of processes in the _ systen 
simultaneously. Through scheduling statements within the 


process descriptions, the run time package provides the 
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control necessary to simulate all the processes in the 
system as if all were executing simultaneously. 

Most discrete simulation languages are based either 
directly or on variations of these methods of systems 
modelling. There are simulation languages that use pure 
event systems, such as SIMSCRIPT; there are simulation 
languages that use pure activity systems, such as CSL; and 
there are simulation languages, such as SOL, SIMULA, MMS, 
and MSS that use the process approach. Some other weli known 
discrete simulaticn languages include CLP, ESP, GASP, OPS, 
SIMON, SIMPAC, and SIMTRAN. Furthur information regarding 
these and other discrete simulation languages can be found 


in [17] and [42]. 


2.2.3 Principal Features of Discrete Simulation Languages 


In general, Simulation languages possess_ some 
properties of general purpose languages but trade speed and 
generality for descriptive power and dynamic control. They 
differ mostiy from the more conventional program languages 
of FORTRAN or ALGOL in that they have mechanisms for 
parallel computation, extended data structuring 
capabilities, and convenient notation for random elements 
within expressions. 

Simulation packages attempt to be as general as 


possible. One of their primary goals is to provide a 
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universal simulation system on which nearly any simulation 
problem can be modelled. Simulation languages vary mainly in 
their appearance to the user, the specific programming 
features they offer, and their basic approach to modelling. 
Below is a list of the important features of general purpose 
discrete Simulation programming languages [3, 30]. All 
languages provide some capability in each area, although the 
degree of capability provided may vary considerably. 

1) Status descriptors - Status descriptors provide 
the user with the means for specifying entities of the 
system and their attributes. They ccnstitute the terms 
in which the user identifies the composition of his 
model. The simulation language automatically maps them 
into a set of data structures which constitute the 
static representation of the model. 

2) Descriptive units - Descriptive units are the 
modules from which a description of the behavior of a 
system is constructed, and are always associated with a 
program of some type. Each unit represents a behavior 
pattern - the essence of some type of activity. it is 
employed repeatedly to simulate individual instances of 
the pattern, that is, the processes or dynamic entities 
of the model. Since most situations of interest are too 
complex to allow a single statement of an overall 
pattern, it is often necessary to isolate subordinate 


patterns which are more accessible to understanding. 
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The overail pattern can then be constructed from the 
joint action of these subordinate patterns. This 
modularity is crucial to the ability to construct 
large, complex models. 

5) >) finerrcontroiye—-8" According’ to)» the particular 
systems modeiling approach taken, the system provides 
an executive cr control program which maintains 
Simulated time, manages scheduling, and coordinates 
paraliel computations. 

4) Data structure = Simulation applications 
require complex data structuring capabilities. 
Facilities allowing flexible structuring of the data 
kase, efficient storage of information, and symbolic 
retrieval of highly structured information are 
necessary to represent complex systems. 

5) Data generation - Facilities are provided to 
generate random numbers from theoretical statistical 
distributions, sometimes with a capability for 
generating correiated data streams. 

6) Observationai facilities - These facilities 
provide the means for data output and monitoring of the 
model. This permits the user to extract data at desired 
peints during the simulation for evaluation purposes. 
Convenient descriptive statistical and graphic 
reporting mechanisms and the ability to use secondary 


storage media are valuable. Facilities capable of 


1, vrs 
ere ee ee 
asbivosy werrys 34 <Aveed bosqtage Qed iedap saetaye, 
witetates dete faawery (o3thes. ao ewljwsexe = Ge ~ 
eesedeeeeos tan «pati«ty . ceqaleg~ gente devalenee 
or . .2qedtbao {es rabiaveg 
Saptseciiqgs aoltriyiit - o ab ota a 
mobs bikscay se fib toys atse sagen ae . erupes 
208d O89 So HO Atr Ho Oe tien arivefic. ewisiliost. 
otiutaye His yeOBsapoude: 0 aa we Meelatite voeed 
yie gosene Te ta! GGaVv fea, = sedvbs ed Lev ehptos 
Merete alien Ueteggeit. bro pneeeseg) | 
oF beptro.w, 920 Beart ienet ~ ombipsipneg eres @ 
{arto attate Year tT aeR Gh TF etext. BIA wobukes ssn teaeg 
set “Qéitidee ee 2 uriy «6 kd, ionontotdzieth a 
ra : ore te node beanEeeI04 ea trias 
. aanssbened eet + pet toLina! Tyantsogseeso (a. 
: eta iaas bes tuqeus ciel meee sax oe 


20 


monitoring the dynamic performance of the model allow 
the user to track down logical errors in its 
description and study its operation. 

7) Specification of initial conditions - These 
facilities provide the means for setting the values of 
all system elements to specify the initial state of a 
model, and for setting the model in motion. 

8) File maintanence - File capabilities allow 
storage of data sets according to complex priority 
rules Or logical conditions; efficient search 
mechanisms are used for locating and removing filed 
information. Such facilities are needed to manipulate 


the large amount of data in complex systems. 


2.3 Interactive Simulation 


In its brief history, the digital computer has been 
used primarily as a calculator with the benefit of a memory 
capability. It was recognized that there were very 
attractive possibilities in combining the computer's 
powerful processing and data handling capabilities with the 
reasoning abilities and insight of the human. Man-machine 
interaction of this sort, however, was rather limited, the 
chief reasons being technological and economic in nature. 
Computers were large and expensive, hardware inflexible, and 


machines not very compatible with man in the area of 


aT. 2 &, 7 
ee 
roa 


as seed. -_ 


4 wha , ~~ a 7 - 
. ; 
Bees _ ‘fost 
ioudereuo sti give fas ooltgi2sett 


3 ana : 5 - 
geeel - enetaigsno Ioisiet te jebwoltinegs 4V 


fa, eeolav ea? padtses aes capes cdf Sbbentq. €329412382- 


» ta o7ete bessiae age Vin <4 02 etasaets ecsege ile 


eer z © wai ‘ 2490 ed . i) age | “Ss 30) be > oe Lebeau: 


seotie. omtaiti@eqao 46:5 - » Adoration Ort UI LL” 


qeiaei3- «£afqao™u *&@ @pitvewis stem”-eeas fo a L0cs 


AIZRow 4iet oe)? an in a0 foot wok, i sol Wa 


hoist an. eas i afteool. “ot Beey sae eve .e0iod. j 


~# faucnens == 


" #Fedugén et is GRADES Ss pi*tLigesi- Bene 2 


7 gcicsebyal® av igpe1880: 


— ; 


wed dad vemeqens dein cts. pwaaipag teitd ard ab 
¥ 4 | : 


(20see 6 20 aideaed o4> it > aiendes es ae ci ttanbeq ite 


yay | oe aqae@d rams Sapocus Bhs ra «eed hae 


: rae)» p | j292. ony Ou baiiin 62 : ong zetia $26205 avis 
vit sean wabdbsscann gertiend ebab bas paLemenoay Lut sow 
= 
= 


7 mussel i nad ott In Piven od: co avitifias | git ORS 
~ ss oT: ts ts pata 
= 


a 
wre es ir rerie 
aes ae eae a oa 
hora el on = 


7 oe nite 
7 : rs 7 


Le 


= aa =a 


21 


efficient communication. The result was that computer 
applications largely ignored the interactive aspect of 
computing for many years. 

But rapid and spectacular advances in both hardware 
and software technology were made in the mid-sixties. Much 
smaller and less expensive electronic components were 
developed. As a result, computers were small and inexpensive 
enough so that a user could afford to have his own digital 
computer. A Lets range of modular computers and special 
purpose feripheral equipment became available so that 
hardware systems could be designed around specific 
applications. The development of visual aids such as 
displays reduced the friction at the man-machine interface 
by providing a common level of communication that was 
acceptabie and natural for the human. 

Of far greater significance, however, was the 
development of time-sharing technigues. Time-shared 
operating systems allow many users to sit at various remote 
terminals connected to the same large scale computer, and 
use it simultaneously. The user has apparent access to all 
the computer's power and facilities and sc can feel as if he 
is at the console of his own moderately sized computer. It 
is now possible to take advantage of the computer's 
interactive capabilities. 

One particularly useful application is in the area 


of simulation. Interactive simulation involves allowing the 
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researcher to play an active role in the simulation process 
as it is cccurring. The value of such interaction is that it 
permits the user to cooperate with the program and share the 
task cf performing the simulation. 

Simulation packages to date have rarely allowed the 
researcher to simulate his models interactively. Using time- 
Sharing for interactive simulation is very new; only 
recently have simulation designers begun to take advantage 
of its power. Only a few interactive simulation systems are 
in the working stages. The most nectable are MIT's general 
purpose OFS+4 package [21] and Michigan's special purpose 
Cellular Space Simulation System [4, 16]. One of the major 
goals of this research is to design an interactive 
Simulation package which will combine the advantages of man- 
machine interacticn and simulation languages into a powerful 
simulation tool. 

An interactive simulation package consists of two 
components: 

1) A simulation language with which the user may 
conveniently specify his model. 

2) A simulation system that allows interaction 
with the model. 

In addition to providing an executive routine which 
manages the execution of the simulation, the simulation 
system provides the facilities for the user to interact with 


the simulation program. It not only allows the user to act 
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as a monitor and observe, verify, and record data, but also 
provides facilities enabling him to modify characteristics 
of the model and redirect the simulation when it strays from 
the desired path. The user is given the ability to access 
and modify various elements of the simulation directly from 
the console. At any point during the simulation he can stop 
to examine or change oie tise data controlling model 
behavior, then continue the simulation from the point of 


interruption. Such capabilities as these make an interactive 


Simulation package extremely attractive and valuable. 


2.4 Requirements of an Interactive Simulation Package 


iy 


The following two sections attempt to provide some 


guidelines as to what basic facilities and capabilities on 


should be provided by the simulation language and simulation 


system of an interactive simulation package. 


2.4.1 Requirements of an Interactive Simulation Language 


The first concern of an interactive simulation 
package is an adequate simulation language. The language 
must provide the general programming capability of any 
programming language so the user can properly describe the 
behavior reguired in his model. It must contain the full 


range of data primitives, a data structuring capability, 
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basic logical and arithmetic capability, conditional and 
unconditional branch facilities, an iterative capability, 
and procedural abilities. 

In addition, the language Must provide the basic 
requirements or features of any simulation language as 
outlined in section 2.2.3. Thus, an interactive simulation 
language must have at least the power of any simulation 
language. But the interactive aspect of the simulation 
package has furthur implications for the simulation 
language. For example, a language designed for an on-line 
environment should include Sacilities for allowing the user 
to communicate directly with a model as it is running. The 
type of debugging facilities provided in an on-line 
Simulation system may take, advantage of this communication. 
Therefore, in the interests of a more efficient system, it 
is desireable to design the system and language together 
rather than independently. 

The iateractive aspect of simulation gives rise to 
the following two requirements for the language: 

1) Appropriate input/output facilities are needed 
for transmitting information about the progress and 
‘status of the model during execution. This is essential 
to allow the user to monitor the system properly and 
determine quickly what actions need to be taken. These 
facilities must have access to _ both system and 


simulation data and be powerful enough to be effective 
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for runtime debugging. 

z), A facility to allow the program to be suspended 
in order to interact with the simulation system is 
valuable. Such an interrupt capability would permit 
execution of the program to be automatically suspended 
at a given point in the program logic, and control 
transferred to the user. A message indicating the 
reascn for the suspension could also be passed. Various 
actions may then be deferred until execution time, 
enabling greater user control. For example, this 
facility would allow decisions as to what statistics 
are to be kept, or which of several options to take at 
a crucial moment, to be made while the user has all 
available information at. his disposal. In addition, 
this facility can be used when unexpected conditions 
are detected, to aid user directed run time error 
correction. _ 

The effectiveness with which the language is 
implemented is of prime importance to the user. Such factors 
as speed of compilation, speed of execution, and 
availability and accessability of storage are important. 
Simulation applications tend to be very expensive, so _ that 
every means should be used to minimize the costs involved. A 
compiler implementation is therefore desirable in order to 
generate object modules. In addition, although assembler- 


level pregramming is not always attractive in terms of 
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programming effort, it is the most desirable level of 
implementation due to its versatility and speed. The 
compiler for the language should incorporate comprehensive 
error checking and provide detailed diagnostic explanations 
for errors. Cross-referencing tables and other compiler aids 


should be available to assist in debugging. 


2.4.2 Requirements of an Interactive Simulation Systen 


The interactive aspect of Simulation means 
additional responsibilities for the simulation system. Not 
only does the simulation system assist in the management and 
control of running the Simulation, it acts as a 
communications device through which the user has access to 
his model. A set cof commands, called the command language, 
allows the user to indicate to the system what action it 
should take. 

It is the duty of the simulation system to manage 
the provision and execution of elements which are common to 
every class of simulation model, items such as time control 
and advance, data generation and manipulation, transmission 
of model status, and related mechanics of model transition, 
The simulation system takes the model components as defined 
through the simulation language and constructs the computer 
representation of the model. It then supervises and controls 


the running of the simulation as directed by the user. 
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The most fundamental requirement of a simulation 
system is a concise, easy-to-use command language which will 
allow efficient access to all system facilities. It must be 
both natural and flexible so that the user can easily and 
quickly specify the stream of commands required to smoothly 
run the simulation. It is important that the user does not 
become bogged down in details of communication due to an 
unwieldy and complex command language. He should be able to 
devote full concentration to monitoring the status of the 
Simulation. The command language must not be a factor in the 
time required to carry out a simulation. 

Another important aspect of an interactive 
Simulation system iS a graphics capability. A graphic 
display is desirable for portraying model behavior because a 
visual representation seems to be useful in understanding 
complex systems. Interrelationships between model entities, 
for example, are much easier to grasp if presented in visual 
form. Equally valuable is the speed and readability of a 
display terminal. Consider first a terminal without a 
graphics capability. To determine model behavior during 
execution, the user usually has to rely only on ae few 
characteristic statistics and indicators. In order to 
determine exactly what is happening in his model, he is 
forced to carry out the simulation in slow motion due to 
long waits for the required data tc be printed on _ the 


terminal. Analysis of such output is not very easy because 
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it is not in the mest natural or readable form. A graphics 
display, however, allows fast, selective scanning of large 
amounts of data. The user is able to carry out a simulation 
much faster and still refer to a far greater amount of data 
- The analysis of displayed output is easier due to the 
superior readability and formatting capabilities of a 
display. 

A graphics terminal is also very versatile. The 
display's ability to quickly switch between various blocks 
of information is very valuable in allowing the user to 
quickly refer to diverse aspects of the simulation. The 
light pen capability present on many consoles can be 
constructively used in many interactive procedures. For 
example, it provides a much more convenient and direct way 
of specifying one or more items than listing them on the 


terminal. 


A list of specific guidelines regarding the most 
important facilities for aiding the user in running his 
Simulation follows. 

1) Of utmost importance is the provision of good 
monitoring facilities to determine the status of the 
Simulation. Comprehensive facilities are needed _ to 
quickly and accurately determine if the Simulation is 
proceeding as expected or if a deviant factor has 


influenced its progress and must’ be pinpointed and 
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corrected (or at least watched more closely). As 
indicated, a graphics capability is useful in this 
regard. 

2) The user should have control over the execution 
of the model in terms of simulated time. He should be 
able to specify the length of simulated time for which 
the simulation should proceed, or allow the simulation 
to run continuously. An interrupt facility should be 
available to suspend the simulation at any point in 
Simulated time. 

5) Sethe system should include the ability to 
restart the simulation by reinitializing the model and 
resetting the system to the original configuration. 

4) Facilities are needed to enable the user to 
save the state of the model on secondary storage at any 
point and time by executing a Single command. 
Conversely, a command should be available to restore 
the simulation to a previously saved state from which 
point execution can be resumed. This means a Simulation 
does not have to be restarted each terminal session but 
can be carried out intermittently. In this way, 
interesting behavior can be examined in detail off-line 
to determine causes, possible effects, and an idea of 
what simulation entities and variables to watch before 
resuming the simulation. Such a facility provides the 


ability to save intermediate states of the simulation 
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run for the purposes of rolling back to previous points 
in simulated time. This is valuable in cases where a 
deviant facter enters the simulation, or important 
behavior was bypassed before its presence was 
discovered. 

5) Facilities should be included to allow 
respecification of system characteristics and model 
parameters. Execution time replacement of routines 
controlling key aspects of model behavior should be 
supported. In this way, the model can be changed in 
mid-stream or restarted with a different configuration, 
enabling easy execution of a series of simulation runs. 

€) Comprehensive facilities for collecting 
Statistical measures of a model's performance must be 
present. Easy collection, display, and analysis of data 
produced during runs of a Simulation program are 
essential for proper evaluation of model behavior. 

7) Provision for hard copy output of statistics 
and data produced during the simulation is necessary 
for detailed off-line examination and a permanent 
record of the simulation. 

8) Facilities should be available to automatically 
generate data that may be required during many 
simulation runs. This removes another programming 
burden from the user. 


9) An effective set of debug facilities including 
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the ability to examine and modify crucial data from the 
console is required to allow error detection and some 
form of error correction during execution. The error 
recovery procedures should maintain the integrity of 
the simulation. Immediate on-line diagnostic 
explanations when an error is detected are desirable. 
10) Facilities are required to make the command 
language flexible and versatile. These include the 
ability to stack commands before execution, cancel a 
list of such ccmmands, and define macro commands built 
up of other commands. This last capability permits easy 
specification and execution of commonly used sequences 
of commands. A facility is also needed to allow 
extensions to the command language in order to provide 
operations dependent on the particular model being 


Ssinulated. 


2.5 Focus of the Present Study 


This chapter has provided a brief background in 
simulation and some useful terms of reference, indicated the 
nature and value of interactive simulation, and pointed out 
some kasic requirements for interactive simulation packages. 
The material presented has dealt with simulation in a fairly 


general frame of reference. The purpose of this section is 
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to narrow the frame of reference and focus on the context of 
the present study. 

Discrete simulation packages attempt to be as 
general as possible in order to allow a maximum user market. 
The object of this research, however, is not to develop a 
general purpose simulation package, but rather a general- 
model simulation package, a general system designed around a 
particular type of simulation problem. The intent is to 
design a system around one general model that can be adapted 
to a iarge number of interesting situations in the specific 
problem area. 

The problem area of concern is that of cellular 
automata and the general model in guestion is the cellular 
space. Cellular sepaces are simple discrete change systems 
and will be described in detail in the next chapter. 
Cellular spaces can be considered as_ activity-centered 
systems in a fixed time step environment. The model consists 
of a specified number of identical entities with known 
relationships and a single activity in which all the 
entities unconditionally and simultaneously participate. Due 
to the simple characteristics of cellular models, much of 
the complexity of a general purpose simulation package is 
unnecessary. For example, the fixed time step approach to 
simulated time removes the need for elaborate scheduling 
specifications. Much of the burden of programming cellular 


models can be removed by allowing the simulation executive 
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routine to handle automatically a greater part of the 
specification and operation of the simulation. A powerful 
special purpose simulation package would be very convenient 
and easy to use, a much more valuable simulation tool for 
cellular research than general purpose simulation systems. 

The facilities required in the simulation language 
are directly dependent on the type of model characteristics 
to be specified. The characteristics of the models in a 
general-model system are much more rigidly defined than for 
general purpose systems. Hence many features incorporated in 
general purpose ianguages are not needed, including some 
listed in section 2.2.3. A general-model package can take on 
a larger part of the workload because the system has a much 
better idea of what the model looks like and how it 
operates; less information is needed to construct the model 
for the user, so much of the model specification reduces’ to 
a setting of model parameters. 

The following items indicate some requirements that 
the simulation language for a general-model simulation 
package should meet due to its specialized nature. 

1) The simulation language should be tailored to 
fit the class of models being simulated so the terms, 
concepts, and constructions are familiar to the user, 
making it a much more hatural and comfortable tool in 
his work. The language should contain special purpose 


statements which allow direct and simple specification 
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of all basic model characteristics. As much of the 
definition of the modei as possible should be done for 
the user so he aeed only give the particulars and can 
concentrate on the unique parts of his system. 

Z) The language should also provide programmed 
routines for standard methods of defining common model 
characteristics. One or more routines should be 
available to handle each aspect of model behavior. 
During model definition, these routines could be chosen 


as options to specify the desired characteristics. 


The facilities provided by an interactive simulation 
system are somewhat independent of the type of models being 
Simulated. The major effect of the special purpose nature of 
a general<~model system involves the amount of power the 
executive program has, and hence the degree to which the 
system can remove the programming burden from the user. The 
more general purpose the simulation package, the fewer the 
common elements which are present, and therefore the more 
limited the power, scope, and utility of the simulation 
system. Conversely, the more special purpose the package, 
the more elements which the system can help manage and the 


more valuable the system becomes to the user. 
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CHAPTER III 


CELLULAR SPACES AND MODELS 


3.1 Cellular Spaces 


A cellular space is an idealized model of a finite 
State machine made up of identical components or cells. 
These functionally simiiar cells are connected to each other 
in a regular manner so that the spacial arrangement of cells 
looks the same no matter which particular cell is 
considered. The cells are traditionally square and arranged 
in matrix form, but may also form triangular, hexagonal, or 
other patterned tesselations. Although triangular or 
hexagonal tesselations look different, they are essentially 
the same and can be translated into an equivalent 
arrangement on a square tesselation. 

Cellular spaces are almost always two-dimensional 
due to the problems of complexity introduced with larger 
dimensions. Stanislaw M. Ulam worked with three-dimensional 
cubical tesselations [44, 38], but very little attention has 


been focussed on non-planar cellular spaces. Due to the 
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Similarity of cells and parallel nature of cellular spaces, 
a plane is usually representative of the space as a whole; 
thus, two dimensions are adequate for most applications. The 
shapes of the plane are also usually very simple. 
Practically all cellular models have been rectangular in 
Shape, with the two notable exceptions [13, 20] both having 
boundaries described by simple convex polygons. 

Each cell has associated with it a set of cells 
called its neighborhood. The behavior of a cell can be 
influenced by the cells in its neighborhood. Typically, the 
set of neighbors of a cell is determined in the same manner 
for every cell, and this determination is often closely 
related to the regular manner of interconnection or topology 
of the space. However, cell-specific and dynamically 
changing neighborhcods are characteristic of some models 
fi2g¢ 

Each cell is characterized by one of a known number 
of different states. The determination of each cell's state 
is specified as a function of its own state and the states 
of the cells in its neighborhood at a particular time. The 
value of this function, called the cell state transition 
function, is the new state of that cell. By calculating the 
function for every cell of the space, a new state assignment 
for the entire space is determined. The sequence of cell 
states of a given cell is called the behavior of the cell, 


and the sequence of state assignments for the entire space 
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is called the behavior of the cell space. 

Cell spaces are discrete in time as well as_ space. 
Their highly parallel formulation is inherently closer to 
the fixed time step method of simulaticn than to the next 
event method. As a result, cell spaces are considered to be 
synchronous - a new state is calculated for every cell at 
the same instant based on the current state of the space, 
and the resulting configuration is the state of the space at 
the next time step. Thus, a pattern of states alters in 
discrete time steps according to a set of transition rules 
(defined in the transition function), that apply 
Simultaneously to every cell. The type of tesselation, the 
allowable set of states, the kind of neighborhood, and the 
transition rules together specify a cellular space. 

Some cellular models [12, 13, 29] require some form 
of random or controlled non-deterministic input to the 
Space. This input is known as external input and can take 
part in determining the next state of any given cell. 
External input allows the model to react to input from the 
environment, generalizing the model from a closed system. 

Complex behavior can be modelled in a Simple 
cellular space. An example will help to demonstrate this 
fact and illustrate the nature of cellular spaces. An 
interesting cellular model is John Conway's simulation game 
called "life" [18, 19]. It is so termed because of analogies 


to the rise, fall, and alterations of a society of living 
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organisms. 

The model consists of a two-dimensional array of 
square cells in a checkerboard configuration. The 
neighborhood constitutes the four orthogonally adjacent and 
four diagonally adjacent cells. Each cell must be in one of 
two possible states - empty or occupied. There are three 
behavioral processes incorporated in the model, each of 
which operates over the span of one time step (termed a 
generation). 

1) Survival - An occupied cell remains in its 
occupied state for the next generation. 

2) Death - An occupied cell reverts to the empty 
State for the next generation. 

3) Birth - An empty cell changes to the occupied 
state for the next generation. 

The transition rules for the space are extremely 
Simple: 

1) Each occupied cell with two cr three occupied 
neighbors survives for the next generation. 

Z) Each occupied cell with four or more neighbors 
dies (from overpopulation). 

3) Each occupied cell with fewer than two 
neighbors dies (from isolation). 

4) Each. empty cell with exactly three occupied 
neighbors is a birth cell, otherwise it remains unoccupied. 


From various starting configuations, the cell space 
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undergoes some unusual, sometimes beautiful and unexpected 
changes. in a few cases the society eventually dies out (all 
cells become empty), although this may not happen until 
after a great many generations. Most starting patterns 
either reach stable figures that cannot change, or patterns 
that oscillate forever. Patterns with no initial symmetry 
tend to become symmetrical; once this happens, the symmetry 
cannot be lost, although it may increase in richness. 

The behavior of each initial configuration is 
determined solely from the transition rules. Consider some 
Simple patterns. A single cell or pair of cells vanishes on 
the first generation. The fate is the same for all but five 
configurations of three cells as illustrated in Figure 3.1. 
Note that their orientation is irrelevant. Figure 3.2 
illustrates the life histories of the surviving four cell 
configurations. 

More complex figures produce many fascinating and 
complex patterns. For example, one five cell configuration 
has been tracked for 460 generations without becoming stable 


or oscillating - its fate is still unknown. 
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The fate of the five 3-cell configurations in "life" 


Examples from the "Game of Life" 


Figure ae 
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The life histories of the five 4-cell configurations in 
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More Examples from the "Game of Life" 


Figure 3.2 
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3.2 Cellular Automata 


Cellular automata theory deals with regular networks 
of similar interacting components, an important class of 
models in many disciplines. The theory promises valuable 
applications in such areas as biological modelling, parallel 
computation, and pattern recognition. 

Cellular automata theory began about 1950 when John 
von Neumann undertook to prove the possibility of self- 
reproducing automata [47]. Given proper instructions, such a 
machine would build an exact duplicate of itself. Von 
Neumann first attained his goal by describing a kinematic 
model of a machine that could roam through a warehouse of 
parts, select needed components and put together a copy of 
itself. 

Stanislaw M. Ulam first thought cf using cellular 
spaces for abstract models of automata. He suggested 
cellular spaces to von Neumann as a way to demonstrate the 
possibility of self-reproducing machines. Von Neumann acted 
on his suggestion and defined the general notion of a 
cellular automaton. He defined a particular two-dimensional 
autcmaton with which he proved the existence of a 
configuration of about 200,000 cells that would be self- 
reproducing. 

Von Neumann showed that his space was both 


computation- and construction-universal by embedding a 
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universal computer-constructor in it. A ccmputer-constructor 
is an automaton capable of simulating a Turing machine, and 
therefore able to perform any computation. From the 
descripticn of any automaton, the computer-constructor can 
construct that autcmaton in a designated empty region of the 
cellular space. Self-reproduction follows as a special case 
of universai construction. 

Von Neumann's cellular space consisted of an 
infinite two-dimensional array of identical 29-state finite 
automata. The state of each depended on its own state and 
the states of its four orthogonally adjacent neighbors. A 
particular state was designated as the quiescent state; only 
a finite number of cells in the infinite array could be non- 
quiescent at any bred The states of the model can be 
descriked as follous: 

Cuiescent state - The quiescent state served as the 
unexcited state; it represented unused cells which, 
upon receiving an input, changed state rather than 
becoming excited. A celi remained quiescent if all its 
neighbor cells were also quiescent. 

Crdinary transmission states - There were eight 
states designated as transmission states which served 
as conduction cells for transmitting an ordinary signal 
through the space. The ordinary transmission cells were 
either excited (which resulted in an output at the next 


time-step) or unexcited (no output); there was a_ state 
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for each of the four neighbors (or directions) to which 
the output could be directed. 

Confluent states - There were eight confluent 
states which allowed logical operations on signals, 
time delays, and split signal paths (transmission 
lines). Confluent cells were excited only if every 
immediately neighboring ordinary transmission cell 
whose output was directed at this cell was excited; 
otherwise they were unexcited. After a delay of two 
time steps output was given in every direction except 
those from which input was received. In order to allow 
this extra time-delay, confluent cells passed through 
two states before outputting, output at next time step, 
and output at present time step. 

Special transmission states - There were eight 
special transmission states which were used to transfer 
a cell from any state to the unexcited quiescent state, 
that is, to "kill" the cell. 

Sensitized states - There were eight states called 
sensitized states which were unexcitable in the sense 
that they did not produce an output but rather changed 
state upon receiving an input. These states represented 
the various transformation stages a cell could be in 
when changing from the quiescent to one of the 
confluent, ordinary transmission, Or special 


transmission states. 
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After setting up the above states and the transition 
rules between these states, von Neumann developed design 
algorithms for necessary parts of his self-reproducing 
automaton. He showed how to form fundamental functional 
units by the appropriate assignment of states to a group of 
cells. Some organs he designed were pulse encoders and 
decoders, recognizers, discriminators, a triple return 
counter, and a coded channel for crossing two signal lines. 
Increasingly higher level units were synthesized from such 
basic units. 

Von Neumann planned his universal computer- 
constructor in two parts: (1) a tape unit consisting of an 
indefinitely expandable tape and its control, and (2) a 
constructing unit with a constructing arm and its control. 
To construct an automaton in an empty region of space, the 
universal computer-constructor builds a path or constructing 
arm to the region, sends signals down the arm to effect the 
construction, and then withdraws the arm. The description of 
the autcmaton is represented as a cell configuration encoded 
on a linear array of cells interpreted as a tape. Von 
Neumann showed that if a description of the constructor 
itself was the input, by having the original tape copied 
onto the newly constructed device and activating it, self- 
reproduction was achieved. 

Von Neumann's celebrated existence proof was the 


first major achievement in cellular automata. At his death 
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in 1957, von Neumann had not completed the entire design of 
his universal constructor. Arthur W. Burks has described the 
essentials of the construction of von Neumann's self- 
reproducing automaton, and shown how the original design 
could ke completed [6]. 

Von Neumann's proof has received much attention, and 
considerable effort has been expended in simplifying his 
system. Burks showed that the design of von Neumann's 
universal computer-constructor could be improved by using a 
double-path construction procedure rather than the single 
path design which von Neumann had originally used. James 
Thatcher (40, 41] worked out a complete and improved design 
of the universal ccmputer-constructor based on a double-path 
construction procedure. 

Much greater simplifications of von Neumann's proof 
were also considered. Edgar F. Codd [8] investigated whether 
a cellular space simpler than that used by von Neumann couid 
be found in which a self-reproducing machine could be 
defined. Through heuristic rather than analytical methods, 
he successfully developed an eight state space with the same 
five cell neighborhood as von Neumann's, that was capable of 
self-reproduction. 

The states were divided into two groups. States 0 to 
3 were inactive or.structure states used as components of 
signal paths. State 0 represented the quiescent state 


equivalent to von Neumann's; state 1 represented a bi- 
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directional signal path; state 2 was used asa sheathing or 
insulating state which surrounded a signal path; state 3 was 
used to form uni-directional paths. States 4 to 7 were 
active or signal states used to build and sheath paths, as 
well as ccnstruct and operate various functional units. 
Sequences of these signals were used to define basic 
operations like extending or retracting paths. Such 
operations were concerned with the construction, 
destruction, and reading of arbitrary (0,1) configurations. 
They were used to build sheathed paths to construction sites 
where the appropriate logic circuits or components were 
built by signals from the paths. The transition function 
allowed for conversion of one Signal state into another so 
appropriate signals could he generated when required. 

Codd adopted von Neumann's approach to demonstrating 
that his cellular system had the desired properties. de 
described the basic building components of his system and 
then used these components to construct the various 
conceptual elements (control, memory, executive) of his 
universal computer-constructor. In addition, Codd 
demonstrated that any eight state/five neighbor cell space 
could ke simulated by a two state/eighty-five neighbor space 
so that his results generalized to such a two state cell 
space. 

Other work has been done in attempting to use a 


simplified space to demonstrate a self-replicating 
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automaton. The latest and best reduction is a four state 
cell space developed by Edwin Rk. Banks at MIT (mentioned in 
[t97}) -. 

The difficulties involved in actually providing a 
universal constructor were examined by Michael A. Arbib [1]. 
Von Neumann had discussed his universal constructor at 
length but had not accomplished defining it. Arbib showed 
that a universal constructor could be described with great 
Simplicity by using much more complicated components than 
von Neumann's 29-state automata. Arbib used a tesselation 
model in which his cells corresponded te modules with four 
orthogonally adjacent neighbors. Each module had input and 
output channels for communicating with its neighbors, and 
weld positions by which it could be joined to them. Welds 
allowed formation of rigid sets of cells such that the 
contents of all welded cells would move ‘from module to 
adjacent module as a single unit. Each module had a bit 
register, 22 instruction registers, and four one-bit weld 
registers corresponding to welds to each of its neighbors. 

Arbib described a special instruction code of eight 
‘basic types of instructions for the internal program of the 
modules. ‘hese instructions were stored in the module's 
instruction registers. Variations of the eight types of 
instructions totalled about 20,000 (<215) which meant that 
each instruction register was 15 bits and the state of the 


module was specified in 335 bits. 
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Arbib showed that any Turing machine could be 
Simulated by a tesselation structure of such modules. He 
used a basic set of six programmable operations, each of 
which he defined by a program of module instructions. The 
longest operation needed 22 instructions and hence each 
module required that number of instruction registers to 
perform any operation. Arbib termed a machine programmed 
with the set of operations a CT machine (Construction 
Turing), and described how a CT machine could be embedded in 
the tesselation. He then showed there was a universal CT 
machine and thus demonstrated the desired universal 


construction machine. 


Another kranch of research in which cellular 
automata is playing an active role is the study of parallel 
systems. In particular, cellular research may be valuable in 
the design of parallel computers capable of simultaneously 
executing many different parts of the same task. 

John H. Holland has made significant contributions 
to this area in his study of parallel systems and 
adaptation. He has developed the concept of iterative 
circuit computers (ICC's), a cellular system composed of a 
large number of identical modules operating in parallel. 
Holland's ICC's are somewhat different from the usual 
definition of a cellular space. Instead of permitting direct 


connections between any two modules, Holland's formalization 
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allows direct connections only between immediately adjacent 
cells and provides for the construction of interconnections 
or paths between any two cells in the space. The 
construction process has constraints which prevent the 
existence of cyclic paths. His motive was to provide a 
cellular space in which adaptive processes could be 
conveniently represented by means of programs. Arbitrarily 
long programs couid be stored in a contiguous set of cells 
and moved to another set of celis in one time step. The 
programs can construct, control, and manipulate other 
programs, and duplicate themselves. 

Holland described a particular type of ICC in which 
each moduie contained a relatively powerful finite automaton 
containing a storage register governing itS own operations 
and its interconnections with its neighbors, as well as 
auxiliary registers storing the information to be processed 
(22]. During a computation, programmed paths could be set up 
so that an active cell could take an operand from another 
cell at the end of the path, and carry out the calculation 
in an accumulatcr of still another cell. Control was 
transferred between cells so a program could be distributed 
around the cellular space, allowing many programs to operate 
in parallel. 

Holland gave a mathematical characterization of a 
broad class of such ICC's [23]. The class of ICC's is the 


set of all automata specified by admissable substitution 
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instances of the quintuple (A, A°, X, f£, P) such that each 
substitution instance designates a distinct network 
organization. The components of the quintuple determine the 
following features of organization: 
- The set A determines the gecmetry of the array 
representing the ICC. For an N-dimensional array, A 
must be a finitely generated abelian group having N 
generators, plus a time step generator on which there 
are no restrictions. The subgroup A' generated by the 
set of N generators indexes the location of a given 
module in the “ees 
> The set A° determines the standard neighborhood 
of the modules or connection scheme of the array; it 
indicates the number and arrangement of modules 
directly connected to a given module. A° must be a 
finite set of elements belonging to A‘; A° specifies 
the arrangement of modules directly connected to a 
module at an arbitrary location. 
- The finite set X determines the storage register 
Capacity of a module. The set of internal states of the 
module is the set S = X x Y where 0 Setermines the 
storage register states of a module and Y consists of 
the possible gate configurations for the paths between 
a module and its neighbors. The state of a module 
depends on its storage register states and its gate 


configuration. 
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- The finite function f determines the permissable 
Operations and commands of the module. It specifies the 
instruction set and thus the possible storage register 
transitions. 

- The finite function P determines the path- 
building or addressing capabilities of the modules. It 
takes the current state and yields a new gate 


configuration. 


Holland has shown that his ICC's contain a subset of 
composition-universal compositions - that is, there exist 
Icc*s in which any machine for computing a computable 
function may be embedded, anywhere in the ICC. Holland has 
used his ICC as a kasis for a theory of adaptive systems 
[24]. He has applied genetic theory to populations of 
programs interacting in an ICC, creating a kind of logical 
autcmaton system operating in accordance with the basic 
principles governing an evolutionary genetic system. He also 
applied these ideas to computers organized in hierarchical 
fashion [25]. He described an adaptive system employing 
hierarchical descriptions in which the ferformance of an 


automaton in a given type of environment could be evaluated. 


Other significant contributions have been made to 
the field of cellular automata. Edward Moore has 


investigated the limits on constructibility in cellular 
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automata {34]. He demonstrated that in a broad class of 
cellular automata, there are finite patterns which he termed 
Garden-of-Eden configurations, that can exist only at time 
zero. He also estabiished a sufficient condition for the 
existence of Garden-of-Eden configurations. John Myhill 
reformulated Moore's condition and showed that it was both 
necessary and sufficient for the existence of Garden-of-Eden 
configurations {36]. Myhill also described a hybrid model of 
a self-reproducing machine combining the features of von 
Neumann's kinematic and cellular models [37]. He treated his 
machines in terms of recursive function theory, and defined 
their properties by axioms concerning their computing powers 
instead of primitive elements or states. Myhill showed that 
in this hybrid cellular system, there exists a sequence of 
machines, each of which constructs its successor and each of 
which outputs more arithmetical truths than its predecessor, 
and no falsehoods. Stanisiaw Ulam and his associates have 
investigated growth patterns of simple cellular automata in 
both two and three dimensions [44, 38]. They considered 
successive generations of simple cellular configurations by 
means of simple recursive relations, and carried out 
mathematical analyses of the resulting highly complex 
figures. These explorations led to generalizations and 
conjectures, some of which have been proved as theorems and 


some of which are still open problems. 
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3.3 Biological Models 


A particularly important and active area of cellular 
research is the modelling of biological systems. The concept 
of arrays of interacting cells forming the structural basis 
of many biological processes has received considerable 
attention in natural systems research. Some of the most 
important cellular models of mammalian physiological 


processes are briefly described in the following sections. 


3.3.1 Flanigan's Model of the Atrioventricular ode. oF the 
Kuman Heart 

In a normally functioning heart, the atrial chambers 
contract forcing blood into the ventricles, before the 
ventricles begin to contract at all. The mechanism whereby 
the ventricles Wait until the atria have finished 
contracting before their own contraction cycle starts is 
largely accounted for by the electrical conduction process 
in the atrioventricular node (A~V node). Larry kK. Flanigan 
{13] developed a cellular model of the A-V node of mammalian 
hearts; he investigated the propagation of excitation from 
the atrial to ventrical edge of the node, and the importance 
of this propagation in coordinated contractions of the 
heart. 

Flanigan's model consisted of a two-dimensional, 


six-neighbor, hexagonal geometry with either a rectangular 
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or funnel shape. The number of cells in the space varied 
from approximately 40 to 300 cells but the results did not 
noticeably vary with size. The top of the network was 
considered to be the atrial end and the bottom the 
ventricular end of the network. The cells on these edges 
were either designated as input cells which received 
external input explicitly controlled to provide desired 
effects of external activity, or output cells whose firing 
patterns were studied. In line with electrophysiological 
studies, the lateral sides of the network were assumed to be 
electrically insulated from other excitable tissues and thus 
unable to receive input from external sources. Such cells 
had connections only with internal cells and so had fewer 
neighbors, ranging frem two to five. The cell properties 
included the following: a variable threshold which governed 
the firing of a cell in that the input stimulation or 
excitation had to ke greater than the threshold in order for 
the cell to fire; a local response which was the result of a 
stimulaticn greater than cell threshold and which served as 
a stimulation to neighboring cells; a firing delay which was 
a dynamic time interval between the stimulation of a cell 
and the firing time of that cell; an absolute refractory 
period (AEP) which was a dynamic time interval following a 
cell firing during which the cell was unexcitable; a 
relative refractory period (RRP) which was a fixed length 


time interval following the ARP during which the cell 
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threshold was above, but returning to, its normal resting 
level. All parameters were constant across cells except for 
the absolute refractory curve which governed the length of 
the AERP and varied between cells. 

A cell fired if the local response it received was 
greater than its threshold. Firing occurred after a time 
delay representing the period during which the cell action 
potential was reaching full amplitude. The length of the 
delay was dependent upon the cell state and the strength of 
Stimuius. Cell firing corresponded to the cell giving a 
ijocal response to its neighbors for a short length of time 
(its duration and maximum value decreasing with increasing 
time delay to make stimulus from non-rested cells less 
effective). 

The network was scanned at specific times and cells 
were sequentially updated to a particular time step. Network 
scans were not made at sequential time steps, but were made 
only at those time steps at which some cell action could 
“sei. Each scan determined the next time step at which a 
scan of the network should occur. The cell data structure 
contained information about the cell's particular absolute 
refractory curve, the particular type of cell (input, 
output, Or neither), and the pertinent times which 
determined when cell updates were required. 

Flanigan's programs were written in the MAD language 


(Michigan Algorithm Decoder), and run on an IBM 7090 
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computer. Three supervisor programs were written to 
investigate particular aspects of model behavior: (1) a 
steady state supervisor which investigated the relationship 
between stimulation cycle length and conduction in the cell 
network; (2) a premature beat supervisor which was designed 
to investigate the effect of premature beats upon the 
network; and (3) a flutter-fibrillation supervisor which was 
designed to investigate the effects of flutter and 
fibrillation type inputs to the model network. The 
Simulation was performed in two stages. First, several 
Simulations on small networks were made to obtain a 
satisfactory set of parameter values (although most 
parameters were preset to enforce ratios known or suspected 
to exist in the A-V node). The second stage used the three 
supervisor programs to investigate the behavior of the model 
with the set of standard parameter values obtained in the 
first stage. 

An important aspect of Flanigan's work was that his 
model supported the contention that cardiac behavior is 
explainable on the basis of a cellular system. Flanigan's 
model closely paralleled a normal heart over the frequency 
of heart beats. With the funnel space, he was also able to 
demonstrate directionality of conduction and test the 
hypothesis that the A-V node's geometry was responsible for 
the higher conduction rate in the ventricular direction than 


in the opposite direction. The validity of his model was 
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most strongly implied by the degree to which it mimicked 
abnormal conditions in the heart, such as premature beats 


and atrial fibrillation and flutter. 


3.3.2 Finley's Neural Network Model 


Marion S. Finley {12] is the most recent of a line 
of investigators who have worked with increasingly 
sophisticated models of neural systems in an attempt to 
provide evidence to support hypotheses regarding the 
existence and development of cell assemblies. A_ cell 
assembly is a collection of neurons that comes to operate as 
a functional unit as a result of its stimulus history. 
Finley carried out a series of simulation experiments with 
varying models, hcming in on a model displaying the desired 
characteristics, properties, and abilities by changing 
system parameters and introducing qualifying properties as 
required. 

Finley's basic model consisted of a two-dimensional 
rectangular array of about 400 (but up to 900) neurons, each 
with from 10 to 60 directed inputs or interconnections. Time 
was quantized and proceeded in discrete time-steps. The 
pattern cf interconnections were initially uniformly and 
randomly distributed, but later experiments used a distance- 
biased interconnection scheme where the probability of a 


neuron having a neighbor a given distance from the cell 
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decreased with increasing distance. Fach interconnection was 
Characterized by a synaptic level which determined the 
excitation level of that input. 

Finley's model was governed by a threshold 
calculation which determined whether or not a given neuron 
would fire (output) or not. A threshold function controlled 
the rate of firing of the neurons of the system by referring 
to several state variables of the cell. The threshold 
function was characterized by a fixed length ARP during 
which the threshold was infinite, and a fixed length RRP 
during which a larger than normal but gradually decreasing 
input stimulus could cause it to fire. 

The cell was characterized by a fatigue value such 
that if a neuron fired at a higher than average rate, 
gradually over a large number of consecutive time steps, the 
fatigue value increased to force the firing rate of the 
neuron down. If the neuron then fired below the average 
rate, tke fatigue value decreased slowly back to the resting 
level. 

A given cell fired if the sum of the threshold value 
and fatigue value was less than the sum of the synapse 
values of connections from cells which had just fired, plus 
an external input value. The external input was used as a 
means of introducing environmental effects into the system 
and the means by which training stimuli were applied. fhe 


synapse values of the interconnections were controlled 
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according to the firing relationships of the neurons. If a 
particular neuron was very instrumental in causing another 
tom fire: (that: is, if it fired, the one receiving its output 
also fired at the next time step), then the synapse value 
between the two was increased. Likewise, if the neuron was 
not effective in causing the other to fire, the synapse 
value between them was decreased. Synapse values could be 
either positive and therefore excitatory, or negative and 
inhibitory. 

Finley used an IBM 7090 computer on which to run his 
Simulation program (written in assembly language). Finley 
investigated the concept of recruitment or correlation of 
neurons; he postulated and verified via simulation that an 
isolated neuron receiving inputs from two sets of neurons, 
one providing a continuous source of background noise, and 
the other acting as a source of information (e.g. periodic 
inputs), would correlate with the information bearing source 
and not with the noisy one. 

He also investigated and simulated stable, steady- 
state behavior in his cellular network. Permanent learning 
of the model resided in the synapse levels so the primary 
condition of steady-state behavior was that it not perturb 
existing synapse levels, possibly disrupting an existing 
cell assembly, nor that it give rise to any new cell 
assemblies. Steady-state acted to preserve the structured 


status gue of the network. 
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Finley claimed that he was very close to 
demonstrating the formation of cell assemblies when his work 


was curtailed. 


3.3.3 Kabrisky's Model of Information Processing in the 
Visual Cortex 

Matthew Kabrisky [29] has investigated the human 
visual system and postulated a model for visual information 
processing in the human brain. His model is basically a 
pattern recognition model which deals with the manner in 
which the brain recognizes an image. Kabrisky conjectured 
that the brain is a two-dimensional pattern manipulator - 
that is, that the computation the brain performs involves 
the modification and transmission of two-dimensional 
patterns in accordance with data previously received and 
stored, and that brain variables are not numbers but the 
models of woridly things which are two-dimensional patterns 
in the cortex. Kabrisky designed his model to demonstrate 
that the krain may be structured to perform two-dimensional 
cross-correlation computations and, hence, that this 
powerful technique might be the kasis of some of the 
computations performed in the brain. The model was based on 
the assumptions that localization of activity occurs in very 
small areas of the cortex (which he termed basic computation 
elements or BCE's), and that a reasonable interconnection 


scheme between entire cortical regions exists so that 
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computation of the two-dimensional cross-correlation of 
stored data with input data is possible. 

Kabrisky*s model can be pictured as a three-planar 
system with the first two planes joined by a network capable 
of shifting and scaling, and the last two planes connected 
by a nugber of many-to-one channels. The first plane 
corresponds to the primary visual input area of the cortex, 
and the second plane to the secondary area. The network can 
be considered to be a number of channels connecting points 
in the first plane with those in the second plane which 
allow shifting, scaling, and rotation of images represented 
on the first plane. The third plane can be considered to be 
a number of images cr patterns stored in memory. Every point 
in the second plane is connected by channels to every 
pattern in the memory. An image on plane-1 undergoes some 
transformation in the network (including the identity 
transformation), and the resulting image on plane-2 is 
joined to an image in plane-3. If cross-correlation is 
successful, recognition has occurred; if not, another image 
in memory is tried or the image undergoes a different 
transformation in the network. 

Kabrisky modelled the secondary cortical plane as a 
10 by 10 array of BCE's. BCE's were the basic components of 
Kabrisky's simulation model; they were not individual 
neurons or brain cells, but cylinders of from approximately 


190 to 350 interacting cells. A collection of BCE's arranged 
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perpendicularly to the cortex formed the computational spots 
in the cortex. The BCE's were connected to their immediately 
adjacent orthogonal neighbors. Time was granularized into 
time steps correspcnding to the largest time interval 
between any two attempts at recognition. Connected arrays of 
BCE'S were capable of carrying out a generalized cross~ 
correlaticn calculation. Kabrisky defined a _ system of 
coupled BCE's which describe the behavior of a generalized 
planar pattern manipulator with a memory, and used this for 
his medel of the cortex. 

Kabrisky's program was written in FORTRAN II and the 
Simulation was performed on an IBM 1620 computer with a 
20,000 character memory. Thus, only very simple patterns 
were used to demonstrate the principles of the computation 
Kabrisky postulated for the brain. The basic patterns used 
in the testing of the model consisted of a single bar or 
square of variable size, position, and rotation. The sheet 
of simulated cortex was able to perform two-dimensional 
cross-ccorrelation and recognize stored patterns with 
variations in location and rotation. In addition, printed 


English letters were used, but proved to be toc complicated 


for the sgall array. 
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3.4 A Simulation System for Cellular Spaces 


It has been indicated that cellular automaton 
systems are not only of interest in their own right, but 
provide a framework for investigating many natural systems. 
Ronald Brender [4] undertook to develop an interactive 
graphics system for simulating cellular models. He designed 
the framework for such a system, tased on the hardware 
available to him at the University of Michigan Logic of 
Computers Group, and successfully implemented a working 


versicn of the system. 


3.4.1 Brender's Cellular Space Simulation System 


The hardware for Brender's system consists of two 
small computers connected by a very fast interface: (1) an 
IBM 1800 with a disk, and (2) a DEC PDP-7 connected to a DEC 
338 CRI display. The 1800 is the logical controller of the 
total system. The data base, user defined routines, and 
related information concerning the cellular automaton model 
under investigation are stored in the 1800. This data 
specifies the cellular geometry, the neighborhood relation, 
the transition function, the initial state, and the history 
of the autcmaton. Information relevant to the simulation is 
displayed on the PDP-7's CRT, including the status of the 


model and parameters of the system, The 1800 communicates 
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with the PDP-7 to generate displays of the cellular space 
state and to accept and act on the user's commands. 

The duties of the PDP-7 consist primarily of three 
parallel tasks: (1) a keyboard command ianguage task, the 
master task which accepts user commands input from the 
keyboard; (2) a display maintenance task, which maintains 
display images and accepts input via the CRT's light pen; 
and (3) a communication task, which handles interaction with 
the 1800. 

In designing his system, Brender's primary goal was 
to maximize the usefulness of the system for supporting 
heuristic and interactive exploration of model behavior. 
Brender's system consisted of a cell space simulation 
language (CESSL) and an interactive run time environment. 
The simulation language is a procedure oriented language 
designed to facilitate the writing and simulation of 
cellular models. It provides a full complement of arithmetic 
and logical operators, a facility for defining structured 
data types and operations on those data types, and various 
data elements and special purpose statements intrinsic to 
the simulation environment in which the translated object 
program is executed. Certain data type names and statements 
are included explicitly for the simulation models. These are 
used in defining the data structures required by the model, 
and specifying entry points into the various user-defined 


routines the system uses to carry out a simulation. 
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The run time environment is the application oriented 
subsystem used for coordinating and directing the running of 
a Simulation. It provides facilities for monitoring the cell 
Space data base, invoking the transitica routine, receiving 
and interpreting user commands, and monitoring the displayed 
description of the cell space state. The specifications 
supplied to the system via the language include the cell 
state data structure, the neighborhood relationship, the 
size and shape of the cell space, and the initial state of 
the cell space. The user-defined programs required by the 
system consist of the peeneaeien function, user-defined 
command, statistics gathering routines, input functions for 
entering data, and the mapping or display functions which 
are used to translate the cell space state into a graphic 
representation for display on the CRT. 

Once a Simulation has been initiated, a number of 
command facilities are available to provide the user with 
contrcel over the course of his simulation. These are 
primarily keyboard commands that are usually executed 
immediately on kbeing read in from the keyboard. The 
capabilities available include the following: compute the 
next time step and display the space under the current Map; 
backup tc the previous time step (only once at a time) and 
display under the current map; compute successive time steps 
and display under the current map until haited; halt the 


Simulation; reset the simulation to its initial conditions; 
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select the indicated map for subsequent displays of the cell 
state; recompute the display for this time step (after 
Changing the map); take a picture of the display; call in 
the system debug facilities; save (or restore) the state of 
the space into (or from) a designated file; call the user- 
defined command; define cell inputs; specify the title to 
use on the cell space display; define a micro-program of 
commands. 

Scme aspects of simulation control are handled via 
the light pen in conjunction with the display system. These 
facilities include age pie and changing parameters 
controlling the form of the displayed image, changing the 
state of a particular cell, and changing the state of a 
group of cells. Lists of commands are displayed and the user 
indicates via the light pen what actions are to be taken. 

A fairly detailed description of Brender's system 
can be found in [4] and [16]. However, since his basic 
system was implemented, work has continued on the simulation 
language, and it has been extended to the point where it now 
stands as an entity independent of the simulation system 


[15]. 
3.4.2 Applications of Brender's Simulation Package 


Brender's system has already proven valuable in the 


investigation of cellular automata systems. Two biological 
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models have been itplemented on his system, and detailed 
investigations of the models carried out with the aid of the 
Capabilities of his package. 

James Mortimer [35] has investigated a class of 
neurophysiological models that have been implemented on 
Brender's system. He developed a cellular model of mammalian 
cerebellar cortex in order to investigate the role of the 
principal neuronal populations in determining the spatio-~ 
temporal response of the cortex to natural and electrical 
inputs. 

In designing his model, Mortimer took into 
consideration detailed physiological data about the 
composition and structure of the mammalian cerebellum, in 
particular, the relative size, distribution, 
interconnection, and excitatory or inhibitory nature of its 
neuronal components. He assumed the accordian-like shape of 
the cerebellum had been unfoided into a plane, and the 
density of Roun OnEe: sueeagnent this cortical sheet was 
uniform. Since the model was designed to investigate 
neuronal activity over significant distances, Mortimer 
modelled a fairly large square section of the cortical 
sheet. The large number of neurons involved meant that 
instead of a cell representing only one neuron, it was 
associated with a small square region of the cortical sheet. 
The model consisted of a 20 by 20 array of such regions. 


Within every region of cortex there are six known functional 
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populations of neuronal elements capable of transmitting 
informaticn. Mortimer considered five of these and 
associated a neuron type with each distinct population. Due 
to these multiple elements within the cell, the components 
of the model were represented by vectors or sets of 
descriptors associated with each neuron type. 

Each neuron type was described by an output state, 
an internal state, an input state, and a neighborhood. The 
output state of a neuron type was the mean firing frequency 
of the population cf neurons represented by that neuron 
type, and the output state of a cell was the set of five 
output states of its neuron types. The output state could be 
regarded as a measure of the excitability of the neuronal 
population. Physiological studies indicate there are six 
different axonai plexuses which supply the synaptic 
terminals that end upon the populations of neurons in the 
cortex. Mortimer represented the input state of a cell by a 
Set) of Ssix5 uinput “states of these plexuses. The internal 
state of a neuron type was a vector of the temporally- 
summated amplitudes of each synaptic input to the population 
of cells represented by that neuron type. The internal state 
of a cell was the set of six internal states of all five 
neuron types. Mortimer defined the neighborhood of a neuron 
type to be the set of cells from which axons could have 
originated which terminate upon that population of neurons 


represented by the neuron type. He described a fixed 
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neighborhcod of cells for each neuron type based upon the 
type and size of distribution of connections. The 
neighborhood specification indicated what neuron types 
within the neighbor cells could send input and the relative 
weights of each input. The neighborhood of a cell consisted 
of the set of neighborhoods of its five neuron types. 

Each cell was associated with a finite automaton 
defined by the input, internal, and output state vectors, 
together with a transition and output function. The 
transition function was entirely specified by the 
interconnection scheme id the time course of the various 
synaptic potentials. It gave the next internal state as a 
set of six functicns of the internal states and input states 
of each neuron type on the previous time step. 
Multiplicative constants on the internal states represented 
the rate of synaptic decay in the respective synaptic 
terminations. 

The output chucraed produced the set of firing 
frequencies of the five neuron types as a weighted sum of 
the components of its internal state (at the same time step) 
minus a constant. The weights represented the synaptic gains 
of the various types of synaptic contacts and the strength 
of each individual contact. The constants played a 
relatively minor role and were similar to thresholds (but 
differed by corresponding to a population of neurons instead 


of a single neuron). The constants were used to establish a 
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baseline of spontaneous activity from which the effects of 
other parameter changes could be assessed. 

The dynamic electrophysiological behavior of the 
cortex was approximated by a discrete sequence of changes in 
the firing or output states of all neuron types in the 
network. Each time step the transition function produced a 
hew state configuration which was mapped by the output 
function into a set of firing frequencies which were assumed 
to persist unchanged for the duration of the time step. The 
delay introduced by the time step was assumed to correspond 
to the sum of the conaletnon delay (travelling time of an 
impulse) and the synaptic delay (time between arrival of an 
impulse at a synapse and the production of the post-synaptic 
potential). Testing showed the simulation output to be 
remarkably insensitive to the duration of the time step. 

The simulation runs were organized into five series. 
The first consisted of a general exploration of the 
parameter space of the model with particular emphasis upon 
the stability of the network and the influence of various 
glotal factors on the behavior of the model. In the second 
set, various anatomical constraints were employed to 
determine a parameter set for the cerebellar cortex of an 
anesthetized cat. The third set of runs consisted of 
experiments on the model of the cat cerebellar cortex 
corresponding to actual neurophysiological experiments in 


order to assess the degree of agreement between the model 
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and known data. The fourth series dealt with the systematic 
variation of model parameters in an attempt to investigate 
the influence of various synaptic connections and other 
Characteristics of the cortical network on the global 
Spatio-temporal behavior of the model. The last series was 
concerned with the response of the cortex to the 
distribution and spacial arrangement of one of the types of 


input. 


Brender's package has also been used in the 
Simulation of a cosselation model for expressing and 
manipulating intercell reiationships. The model was based on 
the model of the metabolism of a single cell, the simple 
unicellular bacterium Escherichia coli, developed by eR. 
Wienberg and M. Burkus [48, 49]. Their object was to 
represent the cell in such a way as to simulate its growth 
in real environments of a living cell, as well as during 
changes from one chemical environment to another. The cell 
was represented in terms of the concentrations of 22 pools 
of chemicals, and equations representing the functional 
relationships among pools. The equations were used to 
predict the change in concentration and total amount of 
chemicals in each pool over time. The pool of proteins was 
divided into different enzyme groups; each enzyme group was 
associated with one group Of chemical reactions responsible 


for converting one pool into another pool in the simulated 
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cell. Repression was employed to control the production of 
its enzymes, and feedback inhibition to control the activity 
of the enzymes already present. 

The simulated environment was a chemically defined 
liquid growth medium at a temperature of 37°C with an 
abundance of oxygen. The changes Simulated in the 
environment were changes in the chemical constituents of the 
media. Three different media were fed to the simulated cell: 
(1) a medium containing glucose, ammonium salt, and 
minerals; (2) a medium containing glucose, minerals, and 
amino acids; (3) a coat containing glucose, minerals, 
amino acids, and nucleosides. With successive additions of 
amino acids, and amino acids and nucleosides, the cell grew 
faster since it had fewer molecules to make on its own. This 
agreed with experimental data from real cells; in the three 
media, the simulated and real cells both took the same times 
to reproduce (50, 28, and 25 minutes respectively). 

The simulation was extended so the cell could also 
adapt genetically. A method of evolving DNA was described 
which included simulation of the evolutionary genetic 
operatcrs of crossover, inversion, and mutation. A _ partial 
Simulation of dominance was also included. 

In order to investigate interactions among 
developing populations of cells, E. D. Goodman, R. Wienberg, 
and R. A. Laing embedded this computer simulation of the 


biochemistry of a single cell (with slight modifications to 


ea. 7. a ~ vs 


ee _ 7 


phohceoe sda lonteo> oF goniidhdnk ; 
vei Qah  2iSeatond. Sem thesaos vas cotolumiz odd 
5 (ete 2°Tf 20 e@Epasteqmes » te @Uiben tyvotp 
) «+ @wFetdate  sepaeis oT .epYR. 2D 
Wi 1S sttons ttetwo>Lapiaqade ott @ aepkal> Srey 
feo, lm slimtaomis of 692 ey~e vbbam gaegw9tth.oord? -stbow 
tle autores ,seemelp. pimisseon oviboa s ‘4 
‘ jaisieein yedoonrt). o¢iedesaem wa Lh ww & {E)- 
isiagia  «augpul iddakatedy aczane a (ct) sebtogo ane 
> “aqolsl bbe etin-oone ecae -oe52g08fSu0 ens .abtos onkas 
wutp thd Ad. wiowe (te GA kee etiee Bam) ,ahbos> embas: 
vag? sa) Sy!) go -i/s of a8) (lade aes Bee os enate aetast 
werd? wt at ioflyo ines wond, oth bkphethaegze dete hasape 
Fh TMGM wile 9009 aed efter Lenaq0ee Apes tonte ant yaid 
| af¥fevi-costdeus aodvata- es be. 338 ¥0¢) anghosgya) 
ihe Oers Thao. oo) of [4tge2 em, 268 apiselumice edt - a 
beviveegh, Gee Ate -petv ove AR yen ne 
Sheep yaewobtdlowe off po gotee hye Tie 
isivvey,, & .«tObTAbY Goo , 072A ceped .sevoqaots & 
| | hiatus bebe ald ‘ven ined FOF 


j 


%, ean stotsouamsat Brags Fan igh, te pera 


75 


allow interactions with other celis) into a _ tesselation 
model programmed on Brender's system [20]. The simulation 
was designed to model the growth of cells on a solid medium 
(e.g. Petri dish), represented by the bottom edge of the 
space. In order to model large colonies of cells using as 
few simulated cells as possible, the simplifying assumption 
of cylindrical symmetry about the centre of a colony was 
used to represent a three-dimensicnal colony ina two- 
dimensional space. Each point in the space represented all 
the cells in the colony at a particular distance from the 
centre and height from tie medium. This resulted in the 
Space being triangular in shape, with a fixed bottom and 
left side (the centre of the colony), and a diagonal edge 
representing the maximum extent of cell growth (upwards and 
to the right) that could be modelled. 

Initial experiments were conducted with a pattern of 
31 cells (which represented several hundred cells in the 
expanded space). The bottom row represented the layer of 
cells resting directly on the growth medium. When the 
Simulation was started, only the lower-left hand point 
represented cells, the other points coming to represent 
cells when adjacent cells divided and daughter cells assumed 
positions at the new points. The interactions among the 
cells in the colony were restricted to the four orthogonally 
adjacent cells which formed the cell neighborhood, and _ the 


other cells with the same radius, since they were all 
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represented by that cell itself. 

Experiments with the simulation mcdel were able to 
demonstrate slow monolayer growth by normal cells, and rapid 
multilayer growth by abnormal cells. Relatively minor 
changes in the metabolic pathways produced markedly 
different growth patterns in the colonies. The biochemical 
equations responsible for such effects were completely 
explicit and could be varied at will, so the biochemical 
basis of global effects in cell growth could be examined in 


detail. 
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CELSIM - A CELLULAR SPACE SIMULATICN PACKAGE 


4.1 The Proposed Simulation Package 


The major goal Of this research is to design an 
interactive simulation package along the lines of Ronald 
Brender's system, for use in cellular research at the 
University of Alberta. The hardware configuration for the 
proposed system consists of two basic components: (1) a 2741 
terminal connected to an IBM 360/67 Digital Computer running 
under the MTS (Michigan Terminal System) time-shared 
operating system; and (2) a CDC Remote Graphics Display 
Terminal Subsystem (hereafter referred to as _ the GRID, 
Graphic Remote Interactive Display) connected to the IBM 
360/67 [9]. 

(1) The IBM 360/67 Under MTS 

MIS is a multi-terminal/batch operating system for 
the IBM 360/67 [45, 46]. An MTS user may operate in 
conversational mode from a typewriter (or CRT) terminal or 


in batch mode (card input). MIS supports a comprehensive and 
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powerful file storage system with public and temporary or 
permanent private user files. The system proposed here makes 
use of koth the batch and conversational modes available, 
and relies heavily on the MTS file storage system. 

Programs and data in MTS are either stored in files 
or else input from hardware devices. To indicate the 
location of a prcegram or data set, the user specifies an 
fdname, the name cf the file in which the program or data is 
stored or the name of the device which is ready to input the 
program or data set. In order to specify the location of 
data at execution time, a Ser enle name called a logical I/0 
unit is used in a program to specify the source of input 
data or the destination of output data. When a program is 
run, it becomes necessary to first specify, for each logical 
I/O unit used by the program, which actual file or device 
should ke used. The MTS commands used by the proposed system 
for compilation of the simulation program and initiation of 
execution of the simulation system make use of logical I/0 
units to provide a high degree of flexibility. 

Many useful facilities are provided for the user as 
an integral part of the MTS system. For example, the file 
editing facilities are valuable for entering and changing 
programs cr data, while a debug package is available for 


program debugging. 
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(2) The GRID 

The graphics terminal consists of a software 
refreshed CRT display and a display controller, interfaced 
to the multiplexcr channel of the IBM 360/67. The working 
area of the CRT screen is regarded as a grid of 1024 by 1024 
addressable positions. The controller is a small general 
purpose computer, modified to include commands to display 
points, symbols, or vectors on the CRI. A sequence of 
display commands is used to define a "display file" which is 
repetitively executed to maintain a picture on the screen. 
The display hardware includes a light pen and an 
alphanumeric and function keyboard. 

The display terminal is coupled to the IBM 360/67 
through a 40.8 kilcbits per second communication line. A 
package of subroutines called GRIDSUB is available to 
assemble in the 360 a display file representing the picture 
the programmer wants to display on the screen [28]. GRIDSUB 
was designed for ering the display file by means of a 
FORTRAN program; the subroutines are written in 360 
assembler language but called from the FORTRAN program. 
Another subroutine allows transmission of the display file 
to the terminal to be displayed. During display, the display 
file is held in the core store of the display controller. In 
addition to the display file, the display controller 
contains a supervisor. The supervisor determines the 


behavior of the terminal in all respects other than _ the 
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content of the screen picture. Initially, the supervisor is 
bootstrapped into GRID core and awaits a message from _ the 
360 to assemble a message for it. Once control has been 
passed to the GRID from the 360, the supervisor assembles 
the sequence of actions made by the display operator into a 
message for transmission to the 360 for processing by the 
user's program (the applications program). The meaning of 
the actions performed by the display operator is determined 
by the nature of the applications program supporting the 
display. The program effectively defines a command language 
for use by the display Eee tare 

In the proposed simulation system, the GRID is used 
primarily as an output device on which to display a visual 
representation of the state of the cellular space, and as an 
input device on which to accept command input from the user. 
These commands are relayed to the 360/67 for execution. All 
command handling is done in the 360/67 where the system 


routines reside. 


The software for the package consists of two major 
components: (1) the compiler for SLANG, the Simulation 
LANGuage (to be identified only with the language defined in 
this paper), and (2) the Monitor (run time control program). 

(1) The Compiler for SLANG 
The user's program is written in the simulation 


language SLANG. The fregram serves two fun3tions: (1) vias 
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defines the required model parameters specifying the size, 
Shape, neighborhood specification, cell data structure, 
method of handling external cells, and other characteristics 
of the model; and (2) it defines the set of routines the 
Monitor requires tc perform the simulation. These routines 
are needed by the Monitor tc perform such basic operations 
as (1) cell space initialization, since the initial values 
of the data structures must be specified, (2) cell space 
transition, since the method of time stepping the cell space 
must be defined, and (3) cell space display, since the type 
of visual Pepressaeeeion of the state of the space to be 
shown on the CRT must be indicated. Other facilities are 
available to the user to allow additional control 
capabilities (e.g. user defined commands and special purpose 
routines). 

The major portion of the total set of routines 
required for performing a Simulation are already provided 
and available to the Monitor. The .user need program only 
those routines defining the behavioral characteristics 
unique to his model. 

The compiler for SLANG is a one pass 360 assembler 
language program [{ 26, 27] which is stored as a set of object 
modules in a file called SLANG. In addition, the compiler 
makes use of a file called ERRFILE which contains the error 
messages used by the compiler for reporting problems to the 


user. 
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When processing the source program, the compiler 
uses the parameter specifications to define the required 
data structures and pointers for use by the Monitor, and 
allocates storage for them. The information supplied about 
the system and data structures (pointers to the symbol 
table, data descriptors, data values, systems flags, 
Simulation routines, etc.) is placed in a block of data 
called KCMMON. The user-defined routines are compiled and 
the resultant object modules placed in a designated file 
along with the data structures and KOMMON. Compilation is 
done separately and reduces a file containing the object 
modules the Monitcr requires. 

(2) The Monitor 

The Monitor is the applications program resident in 
the 360/67 during the simulation. It consists of a set of 
object mcdules residing in a file called CELSIM. These 
modules consist of the code for processing user commands 
from GRID. fhe Monitor controls the execution of the 
Simulation as directed by the user through the command 
language, and prcvides the framework in which a simulation 
model is embedded. The user defined routines are employed in 
running the simulation. 

The major simulation calculation is the procedure 
used for time stepping the cell space. This is an iterative 
process which steps through‘the array of cells maintained in 


the space. The Monitor uses the user defined transition 
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function to change the state of a single cell. The states of 
the neighbors of the cell which is being time stepped are 
made available to the transition function. The cell whose 
State is being updated one time step is called the current 
or central cell. After the central cell has been time 
stepped, its state is transformed into a graphic symbol for 
display on the GRID. This transformation is specified by a 
user defined mapping function. 

These operations are repeated for each cell in the 
space. After ali the cells in the space have been time 
stepped, the graphic evanowe representing the state of each 
cell are displayed in their relative positions on the 
display screen of the GRID. To accomplish this, the Monitor 
uses a subroutine which employs the GRIDSUB routines to 


build the display file and transmit it to the GRID. 


The procedure for initiating compilation and 
execution of the simulation program is straightforward. In 
MTS, execution of a program is initiated by a $RUN command. 
This command requests loading and execution of the program, 
and provides MTS with the name of the file or device where 
the program can be found and the fdnames to be used for the 
references to the various logical I/0 units used by the 
program. The $RUN command used to initiate execution of the 
compiler has the form: 


$RUN <objfdname> <iospecs> PAR=<parameter list> 
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The <objfdname> specifies the object program which 
is loaded. The <iospecs> contain the specifications for the 
assignment of logical I/O units used by the program. The 
logical IyvO units used by the proposed system and which can 
be specified by the user are: (1) SCARDS, which specifies 
the input source file or device (i.e. the user's SLANG 
program); the default values are the input card deck (card 
reader) for batch runs, and the typewriter terminal at which 
the user is signed on for conversational runs. (2), «SPRINT, 
which specifies the output file (i.e. the program listing); 
the default values are cheap eentee for batch runs, and the 
typewriter terminal for conversational runs. (3) SPUNCH, 
which specifies the output punch file (i.e. the object 
modules produced by the compiler); for batch runs, this unit 
defaults to the card punch; it must be specified for 
conversational runs as no default value is assumed. (4) 
SERCOM, which specifies the error message file to which 
error reports are directed: the default values are the 
printer for batch runs, and the typewriter terminal for 
conversational runs. 

The PAR parameter is optional and is used to pass 
data to the pregram on initiation of execution. The 
<parameter list> is terminated by a blank. Options on _ the 
$RUN command can be specified to generate a load map and 
cross reference table if desired. The comment: 


EXECUTION BEGINS 
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is printed if there are no fatal loading errors, and control 
is passed to the entry point of the progran. 

Compilation can be done in either batch OL 
conversational modes; the actual simulation, however, should 
in general be done in conversational mode from a typewriter 
terminal located near the GRID. 

Before initiating execution of the simulation, the 
user must prepare the GRID to communicate with the 360. This 
is accomplished ty entering a bootstrap from a teletype 
attached to the GRID; execution of the bootstrap causes GRID 
to await communication fou the 360 (which will read the 
supervisor into the GRID). The user may then begin execution 
of his simulation. The $RUN statement for loading the 
Monitor with the user program (in object module form), and 
initiating execution of the simulation is very simple. This 
is due to the fact that there is no need to specify any 
logical Iy0O units (since SCARDS, SPUNCH, and SERCOM are not 
used, while SPRINT defaults to the terminal as desired), and 
there is nc parameter specification accepted by the run time 
system. The following example illustrates typical program 
compilation and simulation initiation commands: 

$RUN SLANG SCARDS=INFILE SPUNCH=OBJFILE SERCCM=ERFILE 

This command causes loading and execution of the 
compiler, which exists as a set of object modules in the 
file called SLANG. The input source program is found in the 


file called INFILE, while the output object modules are _ to 
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be placed in the file OBJFILE and the error messages in the 
file ERFILE. The default specification for SPRINT is used 
for the program listing. 

If the compilation was successful, the simulation 
could be started (from the typewriter) by the command: 

SRUN CELSIM+OBJFILE 

Input and output file specifications are not 
required since such files are specified directly by the 
SLANG i170 statements. An error file specification is not 
needed because the Monitor assumes all error reports are to 
be directed to the GRID. The statement loads CELSIM, the 
file containing the object modules for the simulation system 
(the Monitor, Simulation routines, and GRIDSUB routines) 
along with OBJFILE, the file containing the compiled user 
defined routines, KOMMON, and the data structures. If the 
load is successful, control is transferred to the Monitor 
and initialization procedures are started. These consist of 
housekeeping operations Within the Monitor. On completing 
initialization of the system, the Monitor sends the 
supervisor (contained in a file called SUPFILE) to the GRID 
which then awaits input from the user. The message: 

INITIALIZATION COMPLETED ~ CCNTROL RESIDES AT GRID 

is printed on the typewriter terminal by the Monitor which 
then awaits communication from the GRID. 


A block diagram depicting the total system is shown 


in Figure 4,1. 
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Single line connections are unidirectional (downward) 
Double line ccnnections are bidirectional 
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Figure 4.1 CELSIM System Configuration 
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4.2 SLANG : The Simulation Language 


The source statements of the language are read in 
from a file in free format. The basic unit recognized by the 
compiler is the construction, a string of characters 
followed by a semicolon. Program statements consist of one 
or more ccnstructicns. The following table specifies the 
Characters and character groups recognized by the compiler. 
Alphaketics weceeeceeee ABCDEFGHIJKLMNOPQRSTUVWXYZ$! 
NUMELTICS: eccecvecscceece 0123456789. 

QUOTE: wiatic cece ceeccves 

SPACCs sw etestatel es ess eeiee  ' (ise. blank) 

Special Characters .eow +-*/, 33 () =<> 
Alphanumerics ...eee.e-e- alphabetics and numerics 


Non-alphanumerics ..... quote, space, and special characters 


Other graphics are recognized only in character 
string constants fone 1¢{E7#%?2?0 ). The basic lexical unit 
is the atom. An atom is one of the following: (1) any 
special character, (2) any string of alphanumerics delimited 
by non-alphanumerics, or (3) any string of characters 
enclosed in quotes (quotes are represented in a string by 
quote pairs). Space is not an atom but is used as a 
delimiter between atoms. A construction is a sequence of 
atoms ended by a semicolon. A statement consists of a given 


number cf constructions concatenated together and satisfying 
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certain ccnstraints. 

There are several types of atoms. A numeric atom is 
(1) an atcm consisting only of numeric characters and no 
period, treated as a decimal integer constant, (2) an atom 
consisting of numeric characters including a single period 
(decimal point), treated as a real constant, or (3) a 
sequence of characters in FORTRAN-like E-notation. An 
alphabetic atom is an atom of only alphabetic characters; a 
keyword is a predefined alphabetic atom. An anatom is an 
alphanumeric atom beginning with an alphabetic which is not 
a numeric atom nor a keyword. Anatoms are used as data type 
Names, variable names, labels, entry point names, function 
Names, etc. Numeric atoms may contain a maximum of 15 
Characters, alphabetic atoms and anatoms a maximum of 8 
characters, and character strings a maximum of 256 
characters. 

Each atom has one or more attributes associated with 
at which enable the compiler to interpret statements 
properly. Anatoms may have at most one of the attributes 
data type name, entry point name, function name, label name, 
variable name, format name, or equate name. Variable atoms 
(variables) may have other attributes including data type, 
formal parameter, etc. The language allows only one use for 
an atom, which applies at every occurrence of the atom. An 
exception is the atom "-" which can be a monadic or dyadic 


operator, as inferred from the context. 
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Every variable or constant must have as an attribute 
some data type. There are five predefined or primitive data 
types for which the operators of the language are defined. 
(1) A BOCL is a logical variable which occupies one byte and 
may take on the value "F" (for false, represented internally 
by X'0C', where X means hexadecimal) or "T" (for true, 
represented internally by X01"). (2) A CHAR is a one byte 
text value containing one EBCDIC character. (3) An INT is a 
whole-valued number in the range (-127, +127) which occupies 
one byte and serves as a short integer. (Its main use is for 
specifying relative coordinate values of neighborhood cells 
Within the cell data structure; see section 4.2.2.1, the 
NBRS array.) (4) An INTEGER is a whole-valued number in the 
approximate range (-23!, +231) which occupies one word (four 
bytes). (5) A REAL is a fractional-vaiued number in the 
approximate range (1/1079, 1075) and with a precision of 
about 7 decimal places, which occupies one word. 

Constants BERS above types are recognized by their 
lexical properties. A BOOL constant is the atom "T" or “F", 
An INT constant is a numeric atom without a period occuring 
in a context requiring an INT (e.g. in an assignment or 
comparison operation). All integer constants not required to 
be of type INT are translated as INTEGER constants. An 
INTEGER constant is a numeric atom without a period. A REAL 
constant is a numeric atom with exactly one period, or an 


alphanumeric atom in E-notation. A CHAR constant is a string 
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of characters enclosed in quotes. A CHAR constant may 
contain two or more characters although a CHAR variable has 
only one. As many bytes are used for the constant as are 
needed. 

The compiler does not directly provide full scale 
string manipulaticn facilites. However, the system provides 
a special data type calied CHARARAY which describes arrays 
of up to 256 characters. A variable of this type has special 
properties in that in many instances it acts as if it were a 
primitive data type so that the entire string can be 
Manipulated as a unit instead of as an array of one 
Character elements. This capability, along with the ability 
to refer to individual items in the array, makes it possible 
to perform many useful string manipulations. 

There are several types of anatoms with specialized 
properties. (1) LAEEL names take the value of an address in 
the program. A LABEL constant is any anatom that is the 
first atom of a statement followed by a colon (and which has 
not been declared to be a FUNCTION or ENTRY POINT name), or 
any atom following an unconditional branch statement (see 
section 4.2.1.2, the GOTO statement). (2) ENTRY POINT names 
take the value of an address in the program and are 
associated with entry point numbers used by the Monitor as a 
means of subroutine identification. ENTRY POINTS are 
identical to labels except that they must have been 


previously declared as entry point names (see section 
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4.2.1.1, the ENTRY statement). (3) FUNCTION names take the 
value of an address in the program and are associated with a 
parameter list and result description. FUNCTION names are 
identical to labels except that they must have been 
previously declared as entry points to functions. (4) EXTFUN 
QMames are identifiers used for associating an external 
program with the name used as the external reference to 
identify it. These names are used to call an external 
program from a SLANG program and thus, allow access to 
programs written in other languages. EXTFUN names are 
Similar to FUNCTION names except that they are not defined 
Within the SLANG scurce program, but are defined when the 
external programs are made available at loading time. (5) 
FORMAT names are identifiers used for associating a format 
list with an arkitrary output statement. (6) EQUATE names 
are identifiers used to associate a primitive constant with 
a compile time mnemonic name. Direct substitution of the 
value for the aS ae hame is made at compile time. 

Each of the above classes of anatoms requires a 
distinct type of symbol table entry. The total set of anatom 
types requiring entries in the symbol table is: 

Tata Type Names 
Variable Names 
LABEL Names 

ENTRY POINT Names 
FUNCTION Names 
EXTFUN Names 


FORMAT Names 
EQUATE Names 
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In the statement descriptions which follow, two 
conventions are used: 

(1) Curly bracket pairs ( {3} ) are used to 
indicate a repetitive item in the statement. A strict 
maximum number of repetitions is indicated as a superscript 
on the closing bracket. If no superscript is present, there 
is no strict maximum (and as there is no maximum statement 
size, the implementation does not impose an arbitrary 
restriction). The minimum number of repetitions is 1 unless 
the superscript specifies a maximum of 1, whereupon the 
minimum is 0. 

(2) Multiline sguare brackets are used to enclose 
several items to indicate that one (and only one) of the 


items must be used in the statement. 
4.2.1 General Programming Statements 


The sparc wonte of the language are of three types: 
declaration statements, executable statements, and comments. 
The declaration statements assign attributes to atoms or 
give instructions to the compiler or Monitor. The executable 
statements specify operations to be performed. Comments are 
statements which are ignored by the compiler except to be 
printed out in the program listing. Their sole purpose is to 
provide documentation. A comment is a statement which begins 


with an asterisk and ends with a semicolon as follows: 
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* This is a comment statement ; 

Due to the compiler being one pass, a somewhat rigid 
program format is required. The program is divided into two 
sections; (1) the declaration section, and (2) the 
executable section. Only declaration statements and ccmments 
May appear in the declaration section, while only executable 
statements and ccmments may appear in the executable 
section. These two sections are logically divided by the 
BEGIN statement (a declaration statement) : 

BEGIN ; 

This statement allows the compiler to do some 
housekeeping chores (e.g. set up symbol table entries based 
on the cell data structure). 

The first source statement read in is assumed to be 
a declaration statement. The executable section (and the 
source program) is terminated by the ENDPROG statement: 

ENDPROG ; 


An illustration of the structure of a SLANG source 


program follows. 
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SS ee 

{ {— 

| { { 

{ { | 

{ Declaration | { Declaration 
| Statements { | Section 
i { { 

{ | { 

| | | 

BEGIN ; | A 

|- —--------——- | 

{ Co Sata 

\ | { 

| | | 

{ Executable | { Executable 
| Statements | | Section 
| | | 

| { | 

4 | | 

{ ENDPROG ; { --A 
| 


4.2.1.1 Declaration Statements 


The DECLARE Statement 

The DECLARE statement assigns the attribute of a 

particular data type to each of a series of anatoms: 
DECLARE <data type> : <anatom list> ; 

The <data type> must be an anatom which is either a 
primitive data type, a system defined data type (e.g. 
CHARARAY) or a complex data type which has been previously 
defined by the user (see DEFINE statement immediately 
following). The <anatom list> is a sequence of variable 
hames separated by commas to which the data type attribute 


<data type> is assigned. Some examples of DECLARE statements 


are given in Figure 4.2. 
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The DEFINE Statement 

This statement is used to define new data types with 
the option of assigning the newly defined type to variables. 
Two complex data structuring cperations are available. The 
first is the fixed length array: 

DEFINE <anatom> ARRAY <data type> SIZE <integer constant> ; 

The miteroreracian is that <anatom> is defined as a 
data type name which identifies a data structure consisting 
of a fixed number of elements as specified by the <integer 
constant>, all of which are of type given by <data type>. 
The <data type> must be either a primitive data type, a 
system defined data type, or a complex data type previously 
defined ty the user. The second operation provides for the 
definition of a block of contiguous data whose elements may 
be of diverse types. The form of the statement is: 

DEFINE <anatom> BLOCK <data type list> ; 

The interpretation is that <anatom> is defined to 
consist of a Soererirs of contiguous data elements in the 
order they are given in the <data type list>. The data type 
hames specified in the <data type list> are separated by 
commas and each must be either a primitive data type, a 
system defined data type, or a complex data type previously 
defined by the user. 

Arbitrarily compiex data structures can be 
constructed in hierarchical fashion. However, a given data 


structure cannot be given more than one name. AS a 
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convenience, the type defined in the DEFINE statement may be 
assigned to a list of anatoms, as in the DECLARE statement. 
This is accomplished by inserting, before the terminating 
semicolon, a colon and a list of anatoms. Thus the full form 


of the define statement is: 


<array definition> 


DEFINE <anatom> 


<block definition> : <anatom list> ; 
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Some examples of statements defining complex data 
structures can be found in Figure 4.2. 

The language allows for a total of 255 data types 
(including both primitive and user defined data types). Each 
data type name and variable name occupies a four word symbol 
table entry which describes the anatoms' attributes. The 
size of the symbol table governs the number of declarations 
possitle; the programmer can specify the size of the symbol 
table in the parameter field of the $RUN statement which 
initiates execution of the compiler. 

Components of a complex data type may be identified 
by a subscript list (a parenthesized list of expressions 
separated by colons) following the variable name. Subscripts 
are interpreted left to right as identifying successively 
lower data type components in the hierarchical description 
of the data structure. Items in a subscript list may be of 
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truncation) only. All subscript items must be positive and 
the value may not exceed the number of elements in the array 
Or block. Variakle subscripting of arrays is allowed 
(although array bound checking is not done), but 
subscripting of klocks is restricted to constants. This 
requirement removes the need for expensive run time data 
type checking. 

In order to refer conveniently to primitive elements 
of a data structure, fields are used to denote the primitive 
data type substructures of a complex data type. Fields are 
ordered linearly according to the subscripting notation as 
illustrated in Figure 4.2. Each declared variable is 
assigned storage on a fulliword boundary. For the most 
efficient code, whclie word data types (INTEGER, REAL) should 
start on fullword boundaries. The ordering of data types 
within blocks in most cases can and should be specified to 


allow such alignment of fulliword fields. 
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Example Variable Declarations 


DECLARE INTEGER : COUNT, LOOPCNT, TOTAL, PARTSUM ; 
DECLARE BOOL : FLAG1, SWITCH, DECISION ; 

These statements specify each of COUNT, LOOPCNT, 
TOTAL, and PARTSUM to be INTEGER variables, and FLAG1, 
SWITCH, and DECISION to be BOOL varaiables. Storage is 
allocated for each variable. 


Example Data Structure Definitions and Declarations 


(1) Define a data structure consisting of an array 
of 20 real numbers: 


DLEFINE ROW ARRAY REAL SIZE 20 : ANYROW ; 


This statement defines ROW to be a 20 element 
array of REAL numbers, and declares ANYROW to be a 
variable of type ROW. 


(2) Define a 10 by 20 matrix of real numbers: 
DEFINE MATRIX ARRAY ROW REAL SIZE 10 : DATARAY ; 


This statement defines MATRIX to be a 10 by 20 
Matrix of real numbers, and declares DATARAY to bea 
variable of type MATRIX. Note that the array is stored 
in row major crder rather than column major order. 


(3) Define a 10 entry table of 3 columns where 
column 1 contains a real data value, column 2 contains 
3 small integer data values, and column 3 contains a 
flag: 


DEFINE INT3 ARRAY INT SIZE 3 ; 
DEFINE TABENTRY BLOCK REAL, INT3, BOOL ; 
DEFINE DATATAB ARRAY TABENTRY SIZE 10 : OUTCOME ; 


These statements (1) define INT3 to specify an 
array of 3 INTs, (2) define TABENTRY to specify a block 
of data containing a REAL followed by an INT3, followed 
by a BOOL, (3) specifies DATATAB to be a 10 element 
array of TABENTRYs, and (4) declares OUTCOME to be a 
variable of type DATATAB. 


(continued on next page) 


Figure 4.2 Data Structuring Facilities 
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A Detailed Example 


Consider a table of row entries of the type 
defined in the last example as follows (only 2 entries for 
simplicity): 

DEFINE TABLE ARRAY TABENTRY SIZE 2 : RESULTS ; 


Subscripting: 


RESULTS (1) and RESULTS (2) are undefined 
RESULTS (1) and RESULTS (2) are of type TABENTRY 
RESULTS (1:1) and RESULTS (2:1) are of type REAL 
RESULTS (1:2) and RESULTS (2:2) are of type INT3 
RESULTS (1:3) and RESULTS (2: 3) are of type BOOL 


RESULTS (122: 1), RESULIS(1s 222) ;, 
RESULTS (1: 2:3) , RESULTS (29221) , 


RESULTS (2:2:2) , RESULTS (23233) are of type INT 
Structural Diagram: 
TABLE 
i 
{ 
¢ iS as a a eee | 
I i 
| | 
TABENTRY TABENTRY 
{ [ 
{ | 
( SSeS eee 
{ { ! { { i 
j { 1 1 I l 
REAL INT3 BOOL REAL INT3 BOOL 
: atl : : | x 
. { ° ° | ° 
° Ct ° . nas al ° 
° | { i ° ° { { | . 
s { | | . | i { 
e INT INT INT e INT INT IND ° 


e 
e 

e ° e e e e e e ° e 
e ® e @ e e 
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NUMBER 

FIELD 0 & 5 6 qi 8 iz 13 14 15 
DISPLACEMENT 


(IN BYTES) 


Figure 4.2 Concluded 
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The ENTRY Statement 
The ENTRY statement assigns the attribute ENTRY 
POINT to each of a series of anatoms. The form of the 
statement is: 
ENTRY { ( <anatom> <entry point number> ) } ; 
This defines each <anatom> to specify an entry point 
to one cf the user defined routines used by the simulation 
system. The associated <entry point number>'s are assigned 
to each entry pcint name. These numbers are the means by 
which the user will identify his routines at run time. They 
will be used when he is specifying which routines are to be 
used by the simulation system to carry out various fportions 
of the state transition (see sections 4.2.2.2 and 4.3.3.3, 
the SET statement and Set command). There are a total of 50 
such entry points available to the user (numbered 1 through 
50). Each <anatom> defined aS an entry point will appear as 
a label on the first statement of a routine in the 


executable section of the progran. 


The FUNCTION Statement 

The FUNCTION statement assigns the attribute 
FUNCTION name to an anatom and describes the parameter list 
and result of the function associated with that anatom. The 


form of the statement is: 
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RESULT 


FUNCTION <anatom> { ( <parm> <data type> ) } ; 


NORESULT 


ee =e 
[See eas 6) 


The <anatcm> specifies the function name and is the 
label which appears on the body of the function in the 
executable section of the program. There are two classes of 
functions, those which return a value and are used in 
expressions, and those which do not and are executed by 
means cf the CALL statement (see section 4.2.1.2). The 
keyword RESULT/NORESULT determines the type of the function. 
The list of ( <parm> <data type> ) pairs specifies the data 
types of the function variables. If the function returns a 
result, the first element pair of the list specifies its 
hame and type. The remaining values specify the parameter 
names and their types. 

If the function has no result, it acts like a 
separate program; if the function has a result, it acts like 
an operator. The calling sequences for the functions are the 
same: 

<anatom> { < <expression list> > }! 
where <anatom> is the function name and <expression list> is 
a sequence of expressions of the correct data types 
separated by commas (and enclosed in angle brackets), which 
will be used as the initial values of the function 
parameters. Functions may have no parameters and so the 


parameter list may be null. An example follows: 
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Declaration: 


FUNCTION DOTPROD RESULT (RESNAME RESTYPE) 
(PARM1 TYPARM1) (PARM2 TYPARM2) 


; 
Definition: 
CCTPROD ;: 
Statements defining DCTPROD using parameter 
identifiers PARM1 and PARM2, and specifying a 
value for RESNAME 
RETURN ; 
Call Form: 
CCTPROD<PARAMH1,PARAM2> 
where PARAM1 is of type TYPARM1 and PARAM2 is of 
type TYPARM2 (Note that they could be expressions 
rather than simple variable names) 
A function must be defined before it can be referred 
to in a statement. It is recommended that all function 


definitions be grouped together at the beginning of the 


executable section. 


The FCRMAT Statement 

This statement is used to specify a format list for 

the WRITE statement much like in FORTRAN. It has the form: 
FORMAT <anatom> ( <format list> ) ; 

The <anatom> is the format identifier which a WRITE 
statement may include to specify a particular format for the 
output record. The <format list> is a sequence of format 
items or codes separated by commas. The format items are 


analagous to FORTRAN format codes and have the same form 
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{10]. A iimited subset of FORTRAN format codes is supported. 
These consist of the control character code for Spacing 
control, the Hollerith code (character strings enclosed in 
quotes only) for placing alphabetic infcrmation into the 
output record, the xX code for the output of blanks, the A 
code for the output of CHAR and CHARARAY values, the I code 
for the output of INT or INTEGER values, the F and E codes 
for the output of REAL values, and the L code for the output 
of BOCL values. Each of the L, I, F, and E codes can be 
preceded by a field count c (e.g. cIw) which specifies c 
consecutive fields of the type indicated. If the field count 
is omitted, a default of one field is used. 

The processing of the FORMAT statement sets up a 
FORTRAN type format list. The IvO statements are processed 
by FORTRAN 1/0 routines (with some modifications), and thus 
receive the same type of format list. There are no implied 
loops allowed, and only primitive data elements can be 
referred to in the data list. The WRITE statement is used to 
specify one output record and thus, the format list should 
not specify an output string longer than the length that can 
be handled by the output device (e.g. a line file, 255 
characters; a typewriter terminal, 133 characters, including 
carriage control). 

The processing of a WRITE statement with a format 
specification uses the format codes in sequence from left to 


right, matching codes with variables to be printed. If a 
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format code requiring an entry in the print list is 
encountered after all items in the frint list have been 
processed, the action is terminated. This allows Hollerith 
or X codes to be placed in the output record after the last 
item of the print list has been processed, 

If no. format list name is specified in the WRITE 
statement, a default format is used where INTEGERS and INTs 
are output in I6 format, REALS are output in E16.7 format, 
BOOLS are gutput in L6 format, and CHARS and CHARARAYS are 
output as they appear, with no delimiting blanks. 

There is no formatted input statement in the 


language. 


The EXTFUN Statement 

The purpose of this statement is to allow the SLANG 
program to call programs written in other languages. It is 
used to set up access to the external procedures the user 
may wish to execute. The form of the statement is: 

EXTFUN { ( <anatom> <result type> ) } ; 

Each <anatom> in the ( <anatom> <result type> ) pair 
is specified to be an EXTFUN name used within the SLANG 
program to call the appropriate external program. The 
associated <result type> is the data type of the result 
returned by the external program. On the initial call to an 


external program, storage for the address of the external 


symbol is reserved in order to effect branches to the 
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indicated program. The calling procedure format is exactly 
the same as for internal user defined function calls; 
external programs and internal functions are treated in the 
same manner. There is no parameter nhame and type 
specification because there is no parameter type checking on 
calls to external frograms. Standard 360 calling conventions 
are used {27]; a list of the addresses of the parameters is 
constructed and made available to the external program. The 
address of the result is passed back on return from the 
external program. It is entirely up to the programmer to 
ensure that the parameters passed to the external program 
are as expected, and likewise that the result is of the 
correct type. All external programs used must accompany the 
SLANG program in object module form when starting the 
Simulation. Referring to the example initiation procedure 
given in section 4.1, the object modules of the external 
programs could be placed in a file (called PROGFILE in this 
example) and the simulation started by the MTS command: 
$RUN CELSIM+OBJFILE+t PROGFILE 

This facility can be used for many purposes and adds 
greatly to the flexibility and generality of the system. For 
example, it can be used to obtain access to statistical 
routines for generating random numbers or special 


distributions of numbers. 
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The DATA Statement 

The DATA statement allows variables to be preset to 
specific values prior to execution of the program. Variables 
which do not appear in DATA statements have an indeterminate 
initial value. The form is: 

DATA. { ( <anatom> <data. constant> ) } ; 

Each ( <anatom> <data constant> ) pair assigns the 
value BE <tata constant> to the variable <anatom>. The data 
value must be of the same type as the variable name with 
which it is paired. For primitive data type variables, <data 
constant> is just a constant, but for structured variables, 
it is given by a list of constants separated by commas which 
specifies successive fields of the variable. Text constants 
(CHARARAY) are exceptions in that a character sequence can 
be given which the compiler will assign to successive 
locations. 

As an example, for NUMBER of type INTEGER, PI of 
type REAL, MESSACE of type CHARARAY, and ROW an array of 6 
INTEGERS, the form of the DATA statement could be: 

DATA (NUMBER 5) (PI 3.14159) (MESSAGE "SYNTAX ERROR") 
(ROW 10,10,10,20,20,20) ; 

The same constant may be assigned to successive 
fields by the use of a multiplicative factor. A 
specification of the form N*C where N is an integer and C a 
constant, is equivalent to NC's separated by commas. For 
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DATA (ROW 3*10,3*20) ; 

A DATA statement which does not provide enough data 
items to preset all fields of a variable will result in the 
setting of all unspecified fields to machine zero (i.e INT 
and INTEGER zero, BOOL false, REAL true zero and CHAR 
undefined), including the entire data structure. For 
example, 


DATA (ROW) ; 


sets all the entries of ROW to zero. 


The EQUATE Statement 

The EQUATE declaration directs the compiler to 
substitute any occurrence of one atom specifying a primitive 
data type by a primitive constant. Its primary use is to 
allow mneumonic names to be used in place of integer 
constants (especially when used as subscripts). The form of 
the statement is: 

EQUATE { ( <anatom> <data constant> ) } ; 

For each ( <anatom> (data constant> ) pair, the 

compiler substitutes at compile time the value of the <data 


constant> for any occurance of the associated FQUATE name. 


A list of the declaration statements of SLANG 5 


given in Figure 4.3. 
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BEGIN ; 


DATA { ( <anatom> <data constant> ) } ; 


DECLARE <data type> ;: <anatom list> ; 


<array definition> : 
DEFINE <anatom> 


<block definition> : <anatom list> ; 


fr nn a | 
| ee 


ENTRY { ( <anatom> <entry point number> ) } ; 


9 


EQUATE { ( <anatom> <data constant> ) } ; 


EXTFUN { ( <anatom> <result type> ) } ; 


FORMAT <anatom> ( <format list> ) ; 


RESULT 
FUNCTICN <anatom> { ( <parm> <data type> ) } ; 


NORESULT 


P= os oe | 


bee cme oe ee od 


SET <keyword> <keyword specification> ; 


The various forms of the SET statement are shown 
in Figure 4.9. 


Figure 4.3 Declaration Statements of SLANG 


we a = 
7 


- 7 : 
ms Coaeds ger soe ase? a 
si s <Peld gotnne> % & Le ee abe y THAT 
a A a 
‘ 7 oT ‘ a 
\ © Pay cay? ti at soS aetia| 
f tt ‘SOF Bite 
; : Cierl wodena>- 2 | 7 6 aaah AOL 4 
‘ J i er 7 
; |}. ( <ueuusa toLoq, Yraite® <aozeus> )- j <— 
- i “wor —eh> <eedsas> ) } Sh 
: | ( <> (yo ot iqees> <iepane? L ] denies 
; sik qne208> b: <eodeay: pando . 
4 
| waa” - 
its Pp 1S vy foah> nea? ) J ~} a3 i : ae 
} 7.02 a az : 
‘ z A _ ; 


an weg 


- 
ae 
si 


cone wea ES iH Sates shaowye a 


s 
* 


WOls Sk Toate SB aie se aio’ axes OW ade 

ae of 1 us 
- 
7 


7 _ i 
- &. p 7 
i > i 
= - - = ~ ‘. a 
7 * 


351 sisné de Eee €.8 expr . 


4.2.1.2 Executable Statements 


The most important executable statements of the 
janguage are the assignment statement, the conditional 
branch, the conditional statement, the iteration statement, 


and the input/output statements. 


The Assignment Statement 

The basic statement of the language is the 
assignment statement. The most concise way to present its 
acceptable form is via the producticns of a grammar. This 
description is found in Figure 4.4, The following 
observations are made about the syntax. 

1) The nonterminal symbols <variable anatom> and 
<fname anatom> designate the class of variable names and 
function names respectively. 

2) The terminal symbols M$ and $D$ designate the 
set of monadic operators and the set of dyadic operators 
respectively. System defined functions are treated as 
monadic or dyadic operators. The operators available to the 
user are listed in Figures 4.5 and 4.6. 

3) The operator precedence is resolved by using a 
right to left precedence rule whereby the statement is 
parsed and executed from the right. Parentheses are used to 
force calculation of ex pressions forming the left hand 


argument of a dyadic operator. 
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4) Embedded assignment statements of the form 

Roa AS eee = (Ce (DS 2) ).) 

are allowed. The embedded assignment statement must be 
enclosed in parentheses. 

Operators have well defined meanings only when they 
are applied to operands which are of the proper data types. 
Figures 4.5 and 4.6 contain the mode combinations for which 
the operators are defined, and the mode of the result of 
applying the operators. The language maintains sone 


automatic primitive data type conversions as outlined below. 


The follcwing conversions are done automatically by 
the compiler. Note that the internal specification of the 
BOOL constants F andT are equivalent to the INT values 0 
and 1 respectively. 


BOOL -——> INT -——> INTEGER 


The following conversions are done automatically 
only where indicated by the system description. However, 
they are available to the user through the system operators 
FixXS and FLOAT$. Note that REALS are truncated when 
converted to INTEGERS. 


REAL ——> INTEGER 
INTEGER ——> REAL 


CHAR variakles are not convertable. 
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<assignment 
statement> ——-? <left descriptor> = <expression> ; 


<left descriptor> -——> <variable anatom> 


-—-—> <variable anatom> ( <subscript list> ) 


<subscript list> -—> <subscript list> : <expression> 


-——> <expression> 


<expression> ——> <expression> $D$ <expression> 
-—~+> MS <expression> 


—> <descriptor> 


<descriptor> -——> <left descriptor> 
-——> ( <expression> ) 
——> <function call> 


——> ( <assignment statement> ) 


<function call> ——> <fname anatom> 


——> <fname anatom> < <expression list> > 


<expression list> —> <expression list> , <expression> 


——> <expression> 


Figure 4.4 Assignment Statement Syntax 
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Legend: 


OPERATOR 
NAME 


NOTS$ 


BITNOTE 


EXP$ 
ALOG$ 


SQRT$S 
FLOATS 


FIX$ 


CcCOSs$ 
SINS 
ATANS 


TANHS 


{| TYPE {| TYPE 
+--+ 

| Ss) | S 
{ a { a 
| R. { R 
| { 

| 5 | S 
i i { AB 
| R | R 
] ! 

{ B { B 
| | 

| S | S 
| I { rs 
{ { 

| R { R 
| { 

{ R | R 
| | 

{ S i S 
| xt i L 
{ R \ R 
| | 

{ af | R 
| | 

{ | 

| R | 3 
{ { 

| { 

{ R | R 
| | 

| R { R 
| | 

| R { R 
{ | 

| R | R 
| | 


B = BCOL 
Ss = INT 
R = REAL 


OPERAND { RESULT 


"ol 


INTEGER 
The Operand 


FUNCTION OF OPERATOR 


a ee eee 


Figure 4.5 Monadic 


Unary minus (negative of X) 


‘Absolute value of X 


S 


Logical complement of X 


Bitwise complement of X 


e to the power X 
Natural log of X 


Square root of X 


INTEGER to REAL conversion 
of X 


REAL to INTEGER conversion 
raw age 4 


Cosine of X (in radians) 
Sine of X (in radians) 
Arctangent of X (in radians) 
Hyperbolic tangent of X 


(in radians) 
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Legend; B = BOOL CA = CHARARAY 
Ss. = INT CC = Character Constant 
I = INTEGER X = Left hand operand 
R = REAL Y = Right hand operand 
C = CHAR A = An arbitrary data type 
OPERATOR | OPERAND {| RESULT | FUNCTICN OF OPERATOR 
{ ceo! { { 
tis | Sais ] S fee se Plus ny X minus Y 
7 { iia e { zu {| XxX times Y X divided by Y 
SPE { R R | R | X to the power Y 
{ ge | R I 
| Rit i R J 
| i | 
= { A A { A { xX is set equal to Y 
{ Reo L: { R | (a strict data movement 
| ck | I ( for all cases except 
{ Loo es) { Cc | REAL = INTEGER and 
| Sin¢ | S | INTEGER = REAL where 
| CA CGY »| CA | conversion occurs first) 
| \ | 
SMODS { ar pts) | 3 | X modulc Y 
| ae ee | I | 
{ i i 
$EQS FNEF | A A | B { is ,% = Y¥ ? iseie a= 11h? 
S$GTE FGCEF | i f-R | B | Sake 04s, ES ko es, 
$LT$ F$LES | Ra | B ES ge SCL? eo ee ie 
{ Come | { B { Comparison operations 
{ TiC i B | 
{ CALCG { B | 
i EGICAy =| B i 
{ | i 
SANDS | B B | B { Logical operations: 
S$OR$S { { | X "and" Y xX “ort Y 
$XOKF { | luetatexclusive or" Y 
| | { 
$BITANDS | SuGs | S {| The bitwise logical 
S$BITORS | Gas i C { operations: 
S$BITXCHRE ] SHC { S X "and" Y I "ort- x 
| cf | cr ij xX “exclusive or" Y 
| | 
$LS$ Ss s | S { Logical left shift X by Y 
$RSS | Gc =s { e { ‘Logical right shift xX hy 7% 
{ Tyee i I i 


Figure 4.6 Dyadic System Operators 
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The assignment operator causes a movement of data if 
the modes of both operands are the same. However, 
INTEGER/REAL or REAL/INTEGER operand combinations are 
acceptable; truncation or conversion of the right hand 
operand is undertaken in such cases. For CHARARAY/character 
constant combinations, the character constant need not be 
the same length as the CHARARAY variable; truncation of the 
constant cccurs if it is too long, and left justification in 
the variable field if too short. In addition, a CHAR 
variable will accept a right hand operand of type INT. 

The relational operators $EQ$, SNES, $GT$, SCES$, 
$LT$, and $LE$ are defined between operands of non-primitive 
data types and produce a BOOL result. These are performed 
one byte at a time until the result is established. 


Character compariscns are done by a multiple byte compare. 


The Unconditional Branch 
The form of the unconditional branch is 
GOTO <label name> ; 
where <label name> is a LABEL name. Control is passed to the 


statement with the designated label. 
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The Conditional Statement 

The conditional statements are of the following 
forms: 

IF <logical expression> ; 
(Statement Body) 

{ ORIF <logical expression> ; 
(Statement Body) } 

ci Spee. 
(Statement Body) 

ENDIF ; 

The <logical expression> of the IF and ORIF 
statements must yield a BOOL value. The IF statement begins 
the conditional and can occur only once in ae given 
conditional. The ORIF statement is optional but can be used 
Many times in a given conditional. The ELSE statement is 
Optional and can only occur once at the end of the 
conditional. These statements denote mutually exclusive 
sequences of Berceisyts of which the first true condition 
will enable its corresponding body to be executed. If no 
previous IF or ORIF statement yields a true condition, the 
ELSE body (if any) will be executed. The ENDIF signifies the 
scope of the conditional statement. As a convenience for 
identifying ENDIFs, all characters between the keyword ENDIF 


and the terminating semicolon will be ignored. 
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The LCOF Statement 
The LOOP statement is of the forn: 


LOOP <anatom>=<expression> : <increment 
expression> :; <logical expression> ; 


(Statement Body) 
ENDLOOP 3; 

The left hand side of the assignment statement must 
be an anatom of type INT, INTEGER, or REAL, and specifies 
the controlled variable for the loop. The assignment is 
performed and then the <logical expression> is evaluated. If 
true, the iteration is terminated by a transfer to the first 
statement after the ENDLOOP statement; ctherwise, the loop 
body is executed. It is possible that the loop body may not 
be executed at all. The ENDLOOP statement returns control to 
the loop header statement where the controlled variable is 
incremented by the result of the <increment expression>, and 
the termination test is performed again. The LOOP statement 
given above is exactiy equivalent to: 

<anatom> = <expression> ; 
GOTO LABEL2 ; 
LABEL1 : <anatom> = <anatom> + <increment 
expression> ; 
IABEL2 : IF <logical expression> ; 
GOTO LABEL3 ; ENDIF ; 
(loop statement body) 


GOTO LABEL1 ; 
ILABEL3 : CONTINUE ; 


Note that this equivalence implies that the loop 


control variable and increment may be changed in the middle 
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of the loop, that the loop control variable retains its 
value when the loop is satisfied, and that control may be 
transferred into the body of a loop. As a convenience for 
identifying ENDLOOPs, all characters between the keyword 
ENDLOOP and the terminating semicolon will be ignored. 

For example, the loop header: 

LOCP COUNT=X : COUNT-X/Y : TESTVAL $LT$ 100 ; 
indicates that CCUNT is the loop control variable with an 
initial value as specified by the variable X. COUNT is 
incremented upon the completion of each loop by the value of 
COUNT-X/Y (where Y is another variable). The loop will be 
terminated when the variable TESTVAL is less than _ 100 
(presumably the loop statement body will operate on TESTVAL 


according to the value of COUNT). 


The InputyCutput Statements 

Input or output is performed on a particular file or 
device specified een the IyvO statement and consists of a 
stream cf EBCDIC characters directed to or from the 
specified file or device. General I/0 to or from the GRID is 
not supported by these statements due to interfacing and 
display problems. However, character string output to GRID 
is possible by means of the MESSAGE$ function (see section 


Hee we aS) 


= 


byowga® ant or acleasvint Lis pl 1 

shoveuns @@ 1e)4 | Galcthid prtsdnteamse ode saa bai a 
jadbeed Qunt aft catia mt 

= OOF DPE GATTute 3% i-eenod a eeuMID FIO 

sc dots Ridarsey Townes geoi =o eh BeOS, ands ~ woz 

§ WOT Ae Shitgtiny cfs vi etthowgq= ans autor aes 

to wpiev sA9 Fi-qool fige. > sud ead yee we nog 


nay 


m Lkin So0t oeT-.felidsicav * pone Bf ¥ evedv) FEXi-Ta 


- 
‘ 


= 
—— 


~< in 


_ 
oT 


_ a _ 


_ 


= 
> 
oo 
a 


aot sane «eet et Thies sitn hee? eddy aoew - 


’ 


Pe | 
— 


1h0gewr id efyrbao Liliw yho. 70 9uetareS qoot on? os 
(14002 So Godse say od Da 


at Die 
2 


aa 8 


| . _ aIReND RATE auqennsems 
wy erit poitasy & ho. ipa o1rseg, ab daiitao 20 toqna ” 
& toretaianoo bus Snouernse OCF aa bok tionge: ee 


aig woot 6° of ts"xateh azote = wD “7 
boegu ‘i 


- 


_ 7 a 
a 


Sc Giegd ode nom? a0 joe (UNE sii a 
hae oboe Sa Hine as vu an : _ a « . eae 4, 7 
assy oF GUNG Rte DMR RS auton ei gill 

onissen 662) ecisnsinu 2S 


The READ Statement 
The form of the read statement is: 

REAL <idname> <anatom list> ; 
where <fdname> denotes a file or device (an fdname or 
logical IyC unit) specifying the source of the input record, 
and <anatom list> is a list of variables and/or subscripted 
variables separated by commas which specify primitive data 
types (or CHARARAY strings). The elements in the input 
record are recognized as data constants as defined in 
section 4.2. The input routines expect the primitive data 


constants in the order indicated by the variables in the 


<anatom list> and delimited by blanks. 


The WRITE Statement 
The WRITE statement has the form: 
WRITE ( <fdname> { <format name> }! ) 
<expression list> ; 

where <fdname> eee the file or device specifying the 
destination of the output record, and <format name> is a 
FORMAT name and specifies the format list which defines the 
format of the output record. The format specification is 
optional and default format specifications are used if it is 
omitted (see section 4.2.1.1, the FORMAT statement). fhe 
<expression list> is a list of expressions separated by 
commas. Each expression in the list must denote a primitive 


data type or a charater string. The WRITE statement will use 
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Slightly modified FORTRAN output routines. The list of 
output variables must be able to be accomodated within one 


output record. 


The Function/Simulation Routine Control Statements 

The executable section is logically a sequence of 
user defined routines ended by an ENDPROG statement. The 
last statement of the program must be an ENDPROG statement 
to indicate when compilation terminates. Control must not 
flow into the ENDPROG statement at execution time. The user 
defined routines consist of simulation routines and function 
definitions. Each simulation routine has an entry point 
label on its first statement and a RETURN statement in its 
body to return control to the calling program. Each function 
definition has a function entry point label on its first 
statement, contains a RETURN statement in its body for 
returning control to the calling program, and is terminated 
by an FEND statement. Simulation routines are recognized by 
a label (previously declared to be an ENTRY POINT name) on 
the first statement of the routine. Functions are recognized 
by a label (previously declared to be a FUNCTION name) on 
the first statement of the routine. Both types of routines 
can contain any legal executable statement (except for the 
ENDPROG statement). The last statement of a simulation 
routine or the statement preceding the FEND statement of a 


function definition, must be either a RETURN or a GOTO 
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statement which fasses control back into the body of the 
routine or function. Execution cannot flcw past the last 
Statement of the simulation routine nor into the FEND 
statement of a function. 
The form of the RETURN and FEND statements are: 
RETURN ; 
FEND ; 
Due to the amount of code involved in passing 
control ketween programs according to adepted conventions, 
only one RETURN statement per routine or function is 


reccommended. 


The CCNTINUE Statement 

This is essentially a do-nothing statement which 
serves as a convenient place to introduce a label. It has 
the form: 


CONTINUE ; 


The CALL Statement 
The CALL statement has the form: 
CALL <anatom> < <expression list> > ; 
It is used to initiate execution of a user defined function 
which does not return a result (see section 4.2.1.1, the 


FUNCTION declaration). 
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The BASEORG Statement 

This statement is used to respecify a base register 
origin. Only two base registers are assigned for each user 
defined routine sc that a maximum of 8K cf code per routine 
is allowed. Compilation is terminated if the code for a 
routine exceeds the space controlled by the base registers. 
This statement has been included to allow long routines to 
be written by having the user indicate when he has reached a 
logical point in his routine where previous code is no 
longer needed (i.e. will not be branched back into). At this 
point the base registers can be reassigned and the user 
again has an 8K region of core with which to work. The 
BASEORG statement has the form: 


BASEORG; 


A list of the executable statements of SLANG is 


given in Figure 4.7. 
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<assignment statement> (see Figure 4.4) 
BASEORG ; 
CALL <anatom> < <expression list> > ; 
CONTINUE ; 
ENDPROG ; 
FEND; 
GOTO <label name> ; 
IF <logical expression> ; 
ORIF <logical expression> ; 
ELSE ; 
ENDIF ¢ 


LOOP <anatom>=<expression> : <increment expression> 
<logical expression> ; 


ENDLCCE ; 
READ <fdname> <anatom list> ; 
RETURN ; 
WRITE ( <fdname> { <format name> }! ) 


<expression 1list> ; 


Figure 4.7 Executable Statements of SLANG 


 €HolePn yD) vt 


ee ek ls Bh>) 
d coune 
ain: 


TAL - 
2 ws? a se a4 


124 


4.2.2 Simulation Oriented Aspects 


Certain data type names, variables, and declaration 
Statements are available in SLANG explicitly for defining 
the Simulation model. The SET statement is used for 
specifying the simulation oriented aspects of the program. 
The SET statement has many different forms, but all fit the 
general form: 

SET <keyword> <keyword specification> ; 

The <keyword> is a system keyword which indicates 
the particular type of information supplied in the <keyword 
specification>. The following sections will indicate the 
exact form of the SET statement for each type of 
specification. The forms of the SET statement are listed 
near the end of this section in Figure 4.9. 

The declarations concerning the nature of the 
Simulation and the manner of its execution are divided into 
two parts: (1) the definition of the physical 
characteristics of the cell space, and (2) the declaration 
of a series of entry points specifying the computational 


procedure for the simulation. 


4.2.2.1 Physical Description of the Space 


The definition of the size and shape of the cell 


space to be simulated is given by specifying the vertices of 
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a polygon in two space according to the form: 

SET SIZE { ( <coordinate> ) }15 ; 
where <coordinate> is a pair of integers. At least three 
coordinates are required to define the polygon and at most 
15 are allowed. Coordinates of the vertices must be given in 
clockwise order around the figure starting at any point. The 
acceptable shapes of the space consist of all polygons 
descrited by such a sequence of coordinates whereby the 
sequence of xX values is monotonically increasing in the 
range from the minimum X value to the maximum X value, and 
is monotonically decreasing in the range from the maximum X 
value to the minimum X value, and whereby an identical 
condition holds for the sequence of Y values. This group of 
figures includes the set of all convex polygons. For 
example, 

SETSSIZE (10720) (100-20) (— 107-20), -(—10 20); 

specifies a rectangular space of 20 by 40 cells. If no SIZE 
declaration is given, the default is a 20 by 20 square array 


of cells with the lower left hand coordinate at (1 1). 


The CELL Data Structure 

The anatom CELL is reserved as a data type naming 
the data structure for the cells of a simulation. The form 
of the SET statement used to define the exact structure of 


the cell type being simulated is exactly analagous to the 


form of the DEFINE statement: 
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<array definition> 
SET CELL 
<block definition> 


Pe — | 
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If the cell is to be a single primitive data type, 

the user can specify this in the following way: 

SET CELL SAMAS <primitive data type> ; 

This avoids the awkwardness of having to define CELL 

as a one element array, forcing every reference to a cell to 

be subscripted. The user must define the cell data 


structure; there is no default specification. 


The Cell Neighborhood 
The declaration: 
SET NBRHOOD { ( <relative coordinate> ) } ; 
is used tc define the general neighborhood relation common 
to all cells in the space. Each <relative coordinate> is a 
pair of INTs specifying the coordinates (relative to the 
central cell) of a cell which is in the general 
hneighborhcod. For example, the standard five cell 
neighborhood used ky von Neumann and Codd could be defined 
as 
SET NBRHOOD (0 0) (1 0) (0 ~1) (-1 0) (0 1) ; 
If no general neighborhood specification is given, 
the default is a null general neighborhood. 
The neighktors of a cell are referred to as elements 


of the system defined array variable NBRS. The neighbor 
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cells correspond with subscripts of NBRS (e.g.- NBRS(1), 
NBRS(2), «-. ) in the order given in the neighborhood 
definition. Furthur subscripting may be used to refer to 
substructures of the neighborhood cell data states. The size 
of the array NBRS is determined from the general 
neighborhood declaration and from the number of local 
neighkors. 

A tecalenetgihor is a neighbor of a cell which is 
specified by the cell itself. Local neighbors allow for the 
specification of a different neighborhood for each cell in 
the space. Local neighbors are specified in the CELL data 
structure. The appearance of a predefined data type called 
LNEIGH in any substructure in the definition of the data 
type CELL is interpreted as specifying a iocal neighbor of 
the cell. LNEIGH describes a two component array of type INT 
which are interpreted as relative coordinates in the 
coordinate space of the cell space. A variable of type 
LNEIGH specifies a direction vector (with base in the 
current cell), pointing to another cell in the space (the 
local neighbor). These relative coordinates are in the range 
(-127, +127) and thus determine the maximum size of the 
Space, a 128 by 128 array. Since the local neighbor 
specification is used at run time and is part of the cell 
state specification, the transition function can dynamically 
Change a cell's neighborhood. 


The local neighbors specified in the central cell by 
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components of type LNEIGH are automatically accessed and 
made available as neighbors through subscripts of the 
variable NBRS. The local neighbors fcllow the general 
neighbors and correspond with successive subscripts of NBRS 
in the crder in which they occur in the specification of 
type CELL. The NBFS array must always be subscripted due to 
the manner of implementation (instead of pointing to an 
array of cell values, NBRS is a set of pointers to the cell 


values). 


External Input 

External input to a cell at run time is available by 
means of a variable called INPUT. This variable may be 
accessed to obtain input from a particular input stream. An 
input stream 1S a Sequence of integer values defined by the 
user at run time. If the predefined data type EXTINPUT 
(equivalent to an INT) is included in the definition of data 
type CELL, the INPUT variable is assigned the next value of 
the particular input stream specified by the INT value in 
the EXTINPUT field of the current cell. There may be only 
one EXTINPUT field specified in the CELL data structure. 
There are five input streams available to the user (see 


section 4.3.3.3, the EXT command). 


The Cell Space Data Structures 


From the specification of the CELL data structure 
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and the size definition, the cell space data structures are 
created. The system maintains two separate cell space state 
storage areas since all the cells must retain their current 
state vaiues until the next state has been calculated for 
ali the cells. Thus, the transition calculation proceeds by 
accessing neighbcrs and the central cell from one area, 
calculating the new value for the central cell, and storing 
the new value in the other area. These state storage areas 
are called $OLD and $NEW and can be considered to be linear 
arrays of cells beginning with the leftmost cell in the 
bottom row of the cell space, and then proceeding to the 
right and upwards. The coordinates of the cell have no 
bearing on the ordering. 

At run time, the linear equivalent of a cell whose 
coordinates are known may be calculated by means of the 
system function COORDS with a parameter of type COORD 
containing the xX and Y coordinate values. COORD isa 
predefined data type denoting a two element array of 
INTEGERS. The function returns the linear equivalent of the 
cell with those coordinates. For example: 

CELFLD3 = $OLD(COORDS<INDEX>:3) ; 
would store in CELFLD3 the value of the third field of the 
cell in $0LD given by the coordinates specified by the 
values of INDEX (a two element integer array). 

$CLD and $NEW always address different storage 


areas. $NEW always addresses the storage area currently 
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being calculated (during a transition calculation), or the 
area which was last calculated (at all other times). $OLD 
refers to the other area containing the state calculated one 
time step previous to the contents of $NEW. SOLD has no 
meaning, however, immediately after initialization, after a 
Backup command (see section 4.3.3.1), or after a Restore 
command (see section 4.3.3.5), because at these times it 
does not contain the state of the space for the time step 
previous to the one contained in $NEW. To allow the 
programmer to determine when one of these conditions holds, 
a system variable called BACKSW of type BOOL is available 


which is "I" when $0LD has meaning and "F" otherwise. 


The NEWSTATE Variabie 

NEWSTATE is a system defined variable which is 
associated with the cell space data structures. It is pre- 
declared to be of type CELL and takes various meanings at 
different stages of a simulation as follows: 

1) At cell initialization time, it is to be set by 
the user's cell initialization routine to the value desired 
for that cell. 

2) At external cell access time, it contains the 
value of the central cell. 

3) At transition function time, it contains the 
current value of the central cell, and is to be changed by 


the transition function to be the new value of the central 
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cell. 

4) At map time it contains the value of the cell 
for which the map function is to provide a graphic display 
code. 

Thus, NEWSTATE represents an arbitrary cell on which 
the user's routines operate as the simulation routines 


iterate through the cell space time stepping each cell. 


The User Command Variables 
By means of the user command, the experimenter can 
invoke one of the routines which he has defined to carry out 
his personal commands. The user can supply a text string as 
input to these routines. The input string has the following 
general form: 
"<text string>:<data constant list>" 
The <text string> section of the character string is 
a sequence of any legal characters except ":", and can be 
used by the routine as data for a command interpreter in the 
routine. This text string also allows input of CHAR and BOOL 
constants. The <data constant list> section is used solely 
for specifying INT, INTEGER, and REAL data constants. These 
constants are separated by commas in the list and have the 
following form: 
S=<INT constant> ........ for INTs 
I=<INTEGER constant> .... for INTEGERS 


R=<REFAL constant> eceeseee FOL REALS 
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The data string is collected and made available to 
the user routine via the four system defined arrays: 

USERARAY ~ This is a Character array of up to 100 
Characters which will contain the EBCDIC codes for each 
Character in <text string>. 

USERREAL - This is an array of up to 10 REALS which 
will contain the REAL values specified in the <data constant 
list> in the order of their appearance in the string. 

USERINT ~- This is an array of up to 10 INTs which 
will contain the INT values specified in the <data constant 
list> in the order of their appearance in the string. 

USERINTR - This is an array of 14 INTEGERS which 
will contain the INTEGER values specified in the <data 
constant list> in the order of their appearance in the 
string. USERINTR(11) to USERINTR(14) wili contain the length 
of the <text string> section, the number of REALS placed in 
USERREAL, the number of INTs placed in USERINT, and the 
number of INTEGERS placed in USERINTR, respectively. 

This type of input allows a reasonably complex 
command interpreter to be written. These commands might do 
such things as set parameters, or output data to a file. To 
allow for command checking and reporting of ill-formed user 
commands, a user command routine must set the system defined 
BOOL variable USEROKAY either to "TI" to indicate the command 
was successful, or "“F" to indicate it was rejected. A 


rejected command will be reported to the experimenter via a 
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message on the GRID. 


The State Display Variables 

There are two variables concerned with the graphical 
description of the state of the space on the GRID's CRT, 
GRAPHIC and DISPLAY. GRAPHIC is a variable of type INT which 
represents a graphic code corresponding to one of the 
characters which can be displayed on GRID. Figure 4.8 (page 
143) lists the symbols available and their corresponding 
codes. GRAPHIC is set by the mapping routine to the code for 
the symbcoi which is to represent the current cell. DISPLAY 
is an array of INTs, one for each cell in the space, which 
is used by the system display routines when building the 
display file. The value of GRAPHIC set up by the user's 
Mapping routine is entered into the position of DISPLAY 
corresponding to the cell just time stepped. The position of 
DISPLAY corresponding to a given cell can be determined by 


means of the COORD$ function. 


Miscellaneous System Variables 

A system variable called LOCATE is defined to be of 
type COORD. LOCATE takes the coordinates (LOCATE(1) is the 
X-coordinate, and LOCATE(2) is the Y-coordinate) of the cell 
which is specified by NEWSTATE at all times except during 
external cell access time when it contains the coordinates 


of the cell external to the space. 
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NTRANS, NUMCELLS, and NUMNEIGH are all variables of 
type INTEGER. NTRANS is the time step indicator specifying 
the current time step just performed. NUMCELLS is the total 
humber of cells within the boundaries of the space. NUMNEIGH 
is the total number of neighbors in the cell neighborhood 


(both general and iocal neighbors). 
4.2.2.2 Entry Point and Variable Name Specification 


The computational procedure for effecting the 
Simulation of a cell space is controlled by the Monitor. The 
Monitor treats portions of the SLANG program aS subroutines 
controlling one of the options available within the scope of 
the simulation. These options include such things as the 
transition functions, control of external cell accesses, 
initialization of the space, and so on. 

Control cf these options is specified by (1) 
declaring certain labels in the program (by means of their 
corresponding entry point numbers) as entry points to 
simulation calculaticn subroutines to be called when 
required by the Mcnitor according to the function they 
perform, and (2) specifying certain variables as sources of 
special data to be accessed when required. This is 
accomplished by means of the SET statement. The general form 
of the SEI statement for entry point specification is as 


follows: 
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SET ENTRY { ( <keyword> <entry point number> ) }8 ; 

Each ( <keyword> <entry point number> ) pair 
specifies the <entry point number> to correspond to an entry 
point to a particular routine whose function is indicated by 
the <keyword>. There are a total of 8 entry points which the 
user Can specify. All entry point numkers appearing in the 
list must appear in an ENTRY statement to define the 
associated entry point names (see section 4.2.1.1). A 
Similar command is available at run time by which the user 
can Change these specifications (see section 4.3.3.3, the 
Set command). 

Tke general form of the SET statement for variable 
Mame specification is as follows: 

SET VALUE { ( <keyword> <anatom> ) }@ ; 

Each ( <keyword> <anatom> ) pair denotes the 
<anatom> to be a variable name specifying a data value to be 
used under special conditions as indicated by the <keyword>. 
There are only 2 variables which the user can specify. Each 
variable name appearing in the list must be previously 
declared to be a variable of type CELL. 

In the descriptions of the entry points which 
follow, the meaning of the pertinent global variables at the 
time of the calls to these routines is given, along with the 
expected action of the called subroutine. These routines 
have no explicit parameters but use the global variables 


defined by the system, as described in the previous section. 
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Global variables not mentioned in a particular entry point 
description in general have no meaning at that time and 
should not be used. 

Simulation entry points should not be treated as 
labels in the sense that control should not flow into then 
and they should not be addressed by GOTO statements. Control 
is returned from these routines via the RETURN statement 
which must be specified. No return value is required by the 
routines; rather, they must set up the global variables 


indicated. 


Global Initialization 

Upon first entry to the Monitor and after every 
Restart command (see section 4.3.3.5), the Monitor will call 
a global initialization routine (if one is specified). fhe 
purpose of this routine is to allow the user to set up any 
internal variables he wishes, and to do whatever else makes 
sense (such as set up the initial or external cell values) 
before the cell initialization procedure begins. The global 
initialization routine entry point is specified by: 


SET ENTRY ( GLBINIT <entry point number> ) ; 


Initial Cell Specification 
After a global initialization, every cell in the 
Space is set to an initial value. The user may specify a 


variable which will be used to set the initial value of 
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every cell by the following statement: 
SET VALUE ( INITIAL <anatom> ) ; 

The <anatcm> is to be a variable of type CELL. It 
may appear in a DATA statement or be set up by the GLBINIT 
entry. If a more complex setting for the initial value of 
the space is desired (e.g. a different initial value for 
each cell), the user may specify an entry point to a routine 
analagous to a transition function which will be called once 
for each ceil in the space (in sequence from left to fright 
and bottom to top). This is specified by the statement: 

SET ENTRY ( INITRANS <entry point number> ) ; 

On entry to the cell initialization routine 
specified by the <entry point number>, LOCATE will contain 
the coordinates of the cell which is to receive an initial 
value. The initial value calculated by the user program is 
to be stored in the variable NEWSTATE. 

If neither of the above options is specified, the 


default initial value for the space is machine zeros. 


External Cell Specification 

In accessing neighbors of a cell, it is possible 
that a neighbor lies outside the boundary of the polygon 
defining the space. The user is able tc specify the value 
given to cells which lie beyond the physical limits of the 
space. If a constant value is desired for all neighbors 


outside the space, the following statement may be used: 
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SET VALUE ( EXTERNAL <anatom> ) <s 

The <anatom> is to be a variable of type CELL, and 
every time an external cell value is needed, the value of 
that variable wili be used. If a more complex method for 
specifying the value of the external cell is desired, the 
following option is available: 

SEIT ENTRY ( EDGENT <entry point number> ) ; 

in this case, the <entry point number> given must 
specify the entry point to the routine defined to perform 
the desired action. At entry, NTRANS will have the time step 
number, LOCATE will contain the coordinates of the cell 
which is external to the space, and NEWSTATE will contain 
the contents of the central cell. 

There are two possible actions that the external 
cell procedure may take. The first action is to change 
LOCATE to point to a cell which is within the space, thereby 
declaring the neighbcr to be that cell. Otherwise, if on 
return from the edge routine, LOCATE is still not within the 
Space, the neighbor is taken from the EXTERNAL cell 
variable. The external cell routine may directly specify the 
external cell value by storing it in the named external cell 
variable and by not changing LOCATE. 

There is another option available. The specification 

SET ENTRY ( EDGENT <keyword> ) ; 
specifies the external cell routine to be one of three 


system defined subroutines indicated by the <keyword>. The 
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<keyword> may take one of the anatoms XWRAP, YWRAP, or 
TORUS. These routines specify the following mappings of the 
external coordinates into the coordinate range of the space: 

XWRAP - This causes the right edge (increasing X 
direction) to be wrapped around and made adjacent to the 
left edge (decreasing xX direction). Due to the possible 
irregular shape of the space, each row is wrapped 
independently of the others. The xXWRAP option will take 
place only if the Y coordinate is within bounds, otherwise 
LOCATE will not be changed and the external cell value will 
be used. 

YWRAP - The intent for YWRAP is similar to XWRAP 
but due to the data structure used to define the space 
internally, the transformation is not well defined unless 
the shape of the space is vertical or concave to the right 
and left, and the top and bottom edges of the space have the 
same length and X-coordinate range and are parallei to the Y 
axis. Rectangles, the commonest form of space, are well 
defined. If the X coordinate is out of bounds, the change in 
LOCATE is not made and the external cell value is used. 
Spaces for which YWRAP is undefined might achieve the 
desired wrapping by performing a 90° rotation of the space 
and using XWRAP. 

TORUS - The action of TORUS is essentially that of 
performing both wraps at the same time (first YWRAP and then 


XWRAP). ‘he restriction on the shape of the space is the 
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same as for YWRAP. 

A cell external to the space is therefore mapped by 
these routines into the space according to its distance 
beyond the space boundary. The space is wrapped so that the 
cell is mapped relative to the opposite side of the space 
(with respect to a given row or column). If the space is 
rectangular, these options will give the usual cylindrical 
or toroidal geometries in a Cartesian space. If none of the 
above external cell specifications are made, the external 


cell takes a default value of machine zeros, 


User Lefined Commands 

The user way wish to define personal commands with 
which to carry out his particular simulation. By the 
statement: 

SEI ENTRY ( USERCOM <entry point number> ) ; 
the user specifies the <entry point number> to correspond to 
an entry point intc a routine which performs the action he 
wishes. This entry point is entered at run time as a result 
of entering the User command at the GRID (see section 
4.3.3.4). The User command will usually be called between 
transition steps when $OLD and $NEW refer to the previous 
and current time steps respectively. NTRANS will contain the 
number of the transition step last calculated. For a 
description of the type of input to the routine, see section 


4.2.2.1, the User Command Variables. 
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Roane on Functions 

The specification of the transition routine entry 
point is made by: 

SET ENTRY ( TRANS <entry point number> ) ; 

The routine specified by the <entry point number> 
will be called once for each cell in the space. The order of 
calculation is by rows, left to right and bottom to top. At 
the time of the call, NTRANS contains the time step number 
for the transition being performed, LOCATE contains the 
coordinates of the cell to be calculated, the NBRS array 
contains the neighbors cf the central cell, $0LD contains 
the previous time step*s state, $NEW is not useable, and 
NEWSTATE has already been set to the old state of the 
current cell. This means a default state value of the same 
state as the previous time step is given to NEWSTATE. The 
task of the transition routine is to change those parts of 
NEWSTATE which are different for the next time step. Since 
NBRS is an array of pointers only, the values of the 
neighbor cells can be accessed but not changed. 

There are two other routines which may be associated 
with the transition function. Before calculation of a 
transition step, the Monitor can be directed to call a pre- 
transition routine; similarly, a post-transition routine can 
be called once the transition has been completed. These 
routines enable parameters to be set or statistics to be 


collected independently of the transition calculations. The 
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declarations are: 


SET ENTRY ( PRETRANS <entry point number> ) 
( POSTRANS <entry point number> ) ; 


where the <entry point number>s identify the entry points to 
these routines. Upon call to the transition routine, $0OLD 
and $NEW correspond to the previous and about-to-be- 
calculated time steps respectively; NTRANS has just been 


incremented. 


Mapping Functions 

In order to monitor the behavior of the space, an 
image of some aspect of the state of the cells in the space 
is displayed on the GRID. The user must supply subroutines 
to map the state of the cell into information which the 
system routines will use in creating such a display. The 
form of the statement to specify one ot these routines is: 

SET ENTRY ( MAP <entry point number> ) ; 

Each time the display is to be updated, the map 
function is called once for each cell in the space. On 
entry, NEWSTATE has the value of the cell for which a 
mapping is needed, NTRANS has the number of the current 
transition, LOCATE has the coordinates of the cell which is 
to be calculated, and $O0LD and $NEW refer to the previous 
and current time steps respectively. The only action 
expected of the map function is to calculate an integer in 


the range 0 to 61, and to store it in the system variable 
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GRAPHIC. This value will be used tec specify the graphic 
image corresponding to the cell state. The display images 
corresponding to the particular character codes can be found 
in Figure 4.8 below. 

Typically cne map will be wanted for each field of a 
cell's state, but any others may be used as desired (e.g. 
LOL grouping cells according to sets of state 


characteristics). 


Symbol Code Symbol Code Symbol Code 
0 0 K 21 ( 44 
1 1 L 22 ) 42 
Z 2 M 23 < 43 
3 3 N 24 > 44 
4 4 O Zo = 45 
5 2 Pp 26 # 46 
6 6 Q Zul, < 47 
7 7 R 28 = 48 
8 8 Ss 29 A 49 
S 9 rT 30 Vv 50 

blank 10 U 31 + 51 
A 11 V 32 = 52 
B 12 W 33 * 53 
C +3 X 34 ye 54 
D 14 Y 35 ~ 55 
E 15 Z 36 * 56 
F 16 . a > a 
G 1% P 38 ¥ 58 
H 18 : 32 _ a9 
if: 19 : 40 % 60 
“i 20 $ 61 


Figure 4.8 GRID Symbols and Corresponding Codes 
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<array definition> 
SET CELL 


-_——— + 


<block definition> 


SET CELL SAMAS <primitive data type> ; 
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SET ENTRY { ( <keyword> <entry point number> ) }8 ; 
The acceptable forms of the ( <keyword> <entry 
point number> ) pairs are: 


EDGENT <entry point number> ) 
EDGENT TORUS ) 

EDGENT XWRAP ) 

EDGENT YWRAP ) 

GLBINIT <entry point number> ) 
INITRANS <entry point number> ) 
MAP <entry point number> ) 
POSTRANS <entry point number> ) 
PRETRANS <entry point number> ) 
TRANS <entry point number> ) 
USERCOM <entry point number> ) 


FE ee ee cee ete ee ee ee ce a 


SET NBRHOOD { ( <relative coordinate> ) } ; 


SET SIZE { ( <coordinate> ) }15 ; 


SET VALUE { ( <keyword> <anatom> ) }@ ; 


The acceptable forms of the ( <keywWord> <anatom> ) 
pairs are: 


( EXTERNAL <anatom> ) 
( INITIAL <anatom> ) 


Figure 4.9 Forms of the SET Statement 
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4,222.3 Special Simulation Functions 


There are two simulation oriented additions to the 
executable repertoire of the language. These are two 
Simulation functions available to the user to be called by 
means of the CALL statement. 

1) The MESSAGES Function - This function is used 
to send a character string message to be displayed on the 
GRID. The MESSAGES function has one operand, a character 
string of up to 86 elements long. The character string will 
be displayed in the message area of the display until either 
the next message is sent or else until the next command is 
entered from the GRID. Control passes to the next statement 
in the program as soon as the message has been sent to the 
GRID. 

Z) The SUSPENDS function ~- This function is used 
to suspend processing and pass control to the GRID. At GRID 
the user has the full resources of the command language at 
his disposal. This is particularly valuable for debugging. 
The SUSPEND$ function has one parameter, a character string 
of up to 86 elements which is displayed in the message area 
of the CKI. The message should indicate to the user that he 
has been passed control and point out the reason for being 
given control. The SUSPENDS function can be used to indicate 
to the user that an unexpected event has happened and allow 


him to possibly correct it or make amendments to neutralize 
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it before the integrity of the simulation is destroyed. Two 
special ccmmands have been included in the command language, 
one to return control to the GRID to allow processing to be 
continued from where it left off, and the other to abort the 
calculation being undertaken. In addition, several 
indicators are availiable to the user at GRID to indicate 
when a suspension is in effect, and, if the suspension 
occurred during a transition operation, which cell in the 
Space was being time stepped (see section 4.3.3.7, the 


Parameter command). 


Scme of the material of the last few sections has 
been condensed in tabular form. Figure 4.10 contains a 
summary of the entry points, the variables available to each 
routine, and the action required by each. A list of all the 
reserved atoms used by SLANG is given in Figure 4.11. A very 
simple SLANG program corresponding to the "Game of Life" (as 


descrited in secticn 3.1) is given in Figure 4.12. 
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System 
Variable 
Name 


BACKSW 
DISPLAY 
GRAPHIC 
INPUT 
LOCATE 
NBRS 
NEWSTATE 
NTRANS 
NUMCELLS 
NUMNEIGH 
USERARAY 


USERINT 


Figure 4.10 


Variable 
Variatle 
purposes 
Variakle 
Variakle 
Variable 
probably 
Variable 
probakiy 
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may be read 


and used, but not changed 
may be read 


and changed for special 


may be read and changed 

must be set by the user 

is available for reading only; 

not useful at this time 

is available for reading and changing; 
not useful at this time 


"—- A blank entry signifies the variable has no 
meaning at the entry point 


Entry Point 
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Variable/Entry Point Relationships 
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Figure 4.10 
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Keywords Monadic Operators Dyadic Operators 
ARRAY GLBINIT = + 
BASEORG GOTO ABS$ = 
BLOCK iy ALOGS x 
CALL INITIAL ATANS y 
CELL INITRANS BITNOT$ = 
CONTINUE LOOP COSs$ SANDS 
DATA MAP EXP$ SBITNOTS 
DECLARE NBRHOOD FIX$ $BITOR$ 
DEFINE ORIF FLOATS SBITXORS 
EDGENT POSTRANS NEGS SEQS 
ELSE PRETRANS NOoTSs $GES 
ENDIF READ SINS $GT$ 
ENDLOOF RETURN SQRT$S SLES 
ENDFERCG SAMAS TANH$ SLS$ 
ENTRY SET S$LT$ 
EQUATE SIZE SMOD$ 
EXTENT TORUS SNES 
EXTERNAL TRANS Logical Constants $OR$S 
EXTFUN USERCOM $PS 
FEND VALUE aE S$RSS$ 
FILE WRITE F SXORS 
FORMAT XWRAP 
FUNCTION YWRAP 
System Entry Points 
Name Action of the Associated Routine 
EDGENT Change LOCATE to be in the space, or set up 
the external cell value 
GLBINIT Up to the user (e.g. data initialization) 
INITRANS Set up the initial value for the cell at 
LOCATE into NEWSTATE 
MAP Set GRAPHIC to a whole number in the range 
(0,61) for use in setting up the state 
display for the cell in NEWSTATE 
POSTRANS Up to the user (e.g. collect statistics after 
a transition) 
PRETRANS Up to the user (e.g. set switches prior to a 
transition) 
TRANS Change the values in NEWSTATE which are to be 
different for the next time step. 
Use the current NEWSTATE and NBRS. 
USERCCM Up to the user; any personal commands 


(continued on next page) 


Figure 4.11 Reserved Atoms 
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System Functions 


Name Parameter Type Result Type 
COORD$ COORD INTEGER 
MESSAGES CHARARAY ~none- 
SUSPENDS CHARARAY -none- 

System Variables 
Name Type Name Type 
BACKSW BOOL NUMNEIGH INTEGER 
DISPLAY #DISP USERARAY #UARAY 
GRAPHIC INT USERINT #UINT 
INPUT INTEGER USERINTR #UINTR 
LCCATE COCRD USEROKAY BOOL 
NBES #NBRS USERREAL #UREAL 
NEWSTATE CELL SNEW #SPACE 
NTRANS INTEGER SOLD #SPACE 
NUMCELLS INTEGER 
System Data Types 
Name Data Structure Description 
BOOL Primitive data type 
CHAR Primitive data type 
CELL User defined 
CHARARAY Array of up to 256 CHARS 
COORD Array of 2 INTEGERS 
EXTINPUT Identical to INT 
INT Primitive data type 
INTEGER Primitive data type 
LNEIGH Array of 2 INTs 
REAL Primitive data type 
#DISP Array of (NUMCELLS) INTs 
#NBRS Array of (NUMNEIGH) CELLS 
(really just pointers to them) 
#SPACE Array of (NUMCELLS) CELLS 
#UARAY Array of up to 100 CHARS 
#UINT ArrayvoOf7 10"aNTs 
#UINTR Array of 14 INTEGERS 
#UREAL Array of 10 REALS 
Figure 4.11 Concluded 
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*The follcwing statement allows the terms “alive" and "deaa" 
to represent the associated cell state values ; 
EQUATE (ALIVE 1) (DEAD 0) ; 

*The following statements declare the program variables ; 
SET CELL SAMAS INT ; 


9 


DECLARE INT: CCUNT, CROWD, I ; 
DECLARE INTEGER: NUMALIVE, LASTALIV, CHANGE ; 
DECLARE CELL: STARTVAL, OUTSIDE ;: 

DATA (LASTALIV 0) (STARTVAL DEAD) (OUTSIDE DEAD) ;: 


*The next statement associates an entry point number with 
each routine ; 

ENTRY (GENESIS 1) (MAP1 10) (MAP2 11) (STATS 20) ; 

*The following statements specify the space as a 25 by 25 
Square array with the appropriate neighborhood ; 

Seresi zee (0) 0) (0.25), (25.25) (25) 0): 
Sua NBRHOOD] (—1 > =1)> (-1°0)) (=-1) 1)" (0° 1)) UF 1) 7 0) 
CISS ty (OR an); 

*The next statement specifies the initial value of all cells 
in the space and any external cell as "dead" ; 
SET VALUE (INITIAL STARTVAL) (EXTERNAL OUTSIDE) ; 

*The next statement specifies which routines are to be used 
during the simulation calculations ; 

SET ENTRY (TRANS 1) (PRETRANS 20) (MAP 10) ; 
BEGIN ; 

*The user defined routines include the transition function 
(GENESIS), two mapping functions (MAP1, MAP2), anda 
statistics routine (STATS) ; 

GENESIS: CROWD=0 ; LOOP COUNT=1 : 1 : COUNT $GT$ NUMNEIG: ; 
CRCWD=CROWD+NERS (COUNT) ; ENDLOOP ; 

IF (CROWD $EQS 3) $ORH (NEWSTATE SEQS ALIVE) SANDS 
CROWD $EQS 2 ; 
NEWSTATE=ALIVE ; ELSE ; NEWSTATE=DEAD ; ENDIF ; RETURN ; 

*The following mapping routines represent living cells by 
"x" and dead cells by "-" or blank, respectively ; 

MAP1: IF NEWSTATE $EQ$ ALIVE ; GRAPHIC=53 ; 

ELSE ; GRAPHIC=32 ; RETURN ; 

MAP2: IF NEWSTATE $EQ$ ALIVE ; GRAPHIC=53 ; 
ELSE 3; GRAPHIC=10 ; RETURN ; 

*The following routine is used to collect statistics 
regarding increases and decreases in cell populations ; 

STATS: NUMALIVE=0 ; LOOP I=1 : 1: I $GT$ NUMCELLS ; 
NUMALIVE=NUMALIVE+$OLD(I) ; ENDLCOF ; 
CHANGE=NUMALIVE-LASTALIV ; 

WRITE (STATFILE) NTRANS, NUMALIVE, CHANGE ; 
LASTALIV=NUMALIVE ; RETURN ; 
ENDEFROG ; 


Figure 4.12 SLANG Program for the "Game of Life" 
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4.3 The Run Time Control Systen 


Once a simulation has been initiated and control 
transferred to GRID, a number of commands (the command 
language) provide the user with control over the course of 
the simulation. The run time control system interfaces with 


the user through the GRID. 
4.3.1 The Facilities at GRID 


The commands for controlling the progress of the 
Simulation are entered via the GRID's alphanumeric keyboard. 
Individual commands are not executed immediately upon being 
entered, ktut are stacked in a command list. The command list 
has a ccunterpart called the message vector which is the 
encoding cf the user's commands used by the supervisor to 
send to the Monitor. This message vector is a sequence of 
entries defining the commands specified by the user's 
actions at the GRID. An end of message indication is given 
to send the assembled message to the 360 for processing by 
the Monitor (i.e. for execution of the encoded commands). 
Normally, if a command requiring transfer of control from 
the 360 to the GRID is executed, the remainder of the 
commands in the message vector are ignored. The next command 
executed by the Monitor is the first command in the message 


vector received from the GRID when control is again 
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transferred to the 360. If the message vector is empty upon 
transmission to the 360, control is simply returned to the 
GRID by the Monitor. 

Commands consist of three character alphabetic 
hames, optionally followed by a sequence of one or more 
parameters, and terminated by striking the HERE key on the 
alphanumeric keybcard. The action of the HERE key is to 
delimit the commands within the message vector; it serves to 
Start each command as a new (and thus separate) character 
string in the message vector. 

Parameters consist of either alphanumeric characters 
or light pen picks. Command keywords and aiphanumeric 
parameters are delimited by blanks. Unreccgnizable commands 
or commands with improper parameters result in error reports 
and the deleticn of the remainder of the message vector, 
with control returning to the GRID. 

The elements forming the message vector represent 
the following user actions: 

(1) Typing strings of alphanumeric characters, 
interpreted as a sequence of commands with accompanying 
parameters where appropriate. Alphanumeric input is 
displayed as a string of characters on the screen in a 
reserved portion of the CRT denoted as the command input 
area (see section 4.3.2 for a description of the display 
screen). This area can seceanoidtd two lines of up to 86 


Characters of alphanumeric input each. If the line length on 
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the screen is exceeded, characters are displayed from the 
left edge on the same line. To avoid cverlapping, the RETURN 
key on the keyboard can he used to shift the display of 
alphanumeric input to the beginning of the second line. A 
cursor (underline character) is displayed on the screen to 
indicate the next available position in the command list. As 
each Character is typed, it is added to the message vector 
as well as being displayed in the command list at the 
position on the screen indicated by the cursor; the cursor 
is moved ahead one character position. A positive indication 
of delimiting a command with the HERE key is given by having 
a semicolon displayed in the next available character 
position in the command list. 

(2) Picking displayed items with the light pen, 
interpreted as specifying parameters of commands. Light pen 
input is generated by picking items displayed on the screen. 
A light pen pick is indicated by a blinking arrow appearing 
above the object chosen. Appropriate information regarding 
the displayed item is placed in the message vector to allow 


identification of the item by the Monitor. 


In additicn to these types of inputs, the user has 
available certain keyboard keys which have some special 
meaning within the supervisor. These are used to supply 
command editing features concerned with providing a flexible 


method of command input. (1) The RESET key causes the entire 
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message vector ccmposed so far to be deleted without 
transmission to the 360. The command list displayed on the 
screen is also removed. (2) The Delete key causes deletion 
of the last component in the message (i.e. a character 
string or light pen pick). It can be used to delete an 
entire alphanumeric command from the message vector and 
delete its counterpart on the screen. (3) The BACKSPACE key 
causes deletion of the last unit of the current component of 
the message vector (i.e an alphanumeric character or light 
pen pick). The deleted character will also be removed from 
the screen. (4) The SEND key causes the message vector 
composed te the user to be transmitted to the 360 for 
execution by the Monitor. When control is returned from the 
360, the command list and message vector are emptied and the 
cursor is positioned at the beginning of the first line of 
the command input area. 

The message vector has two restrictions: (1) tables 
associated with the message can accomodate a maximum of 20 
entries (character strings or light pen picks); and (2) the 
sum of the elements in the message vector cannot exceed 200 
entries (each light pen pick constitues about 6 entries, 
each alphanumeric string constitutes 2+[N/2] entries where N 
is the number of characters). If these limits are exceeded 
the message must tke deleted and a shorter one re-entered. 
These restrictions are somewhat artificial and can be made 


less restrictive; as a rule, however, a limited number of 
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commands will be stacked before being executed. The only 
time when the limits may be approached is when long 
parameter strings are required. 

The Monitor does comprehensive error checking with 
respect to analyzing the command list and carrying out the 
requested actions. If an error condition is detected (e.g. 
an unrecognizable command keyword or parameter, a syntax 
error, cr a request for an illegal ccmmand), the Monitor 
reports the situation to the user at the GRID via an 
appropriate message in the message area of the screen. If 
the situation can ke neutralized by the Monitor without user 
interaction, this is done (e.g. for a syntax error, ignore 
the command altogether, or if the error is discovered during 
execution of a command, abort furthur processing of the 
command avis the status of the system will not be 
compromised). 

For machine level error conditions (e.g. data 
overflow cr underflow, divide exceptions, etc.), the error 
interrupt is trapped, and control is passed back to _ the 
Monitor. An error report indicating the error condition is 
generated and sent to GRID, transferring control to the user 
there. The full resources of the simulation system are 
available to the user and he is free to do whatever he deems 
appropriate. If the error was not crucial to the simulation 
(e.g. the error occured acu execution of a non-essential 


User command routine), the user may continue the simulation 
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as if nothing happened. However, if the seriousness of the 
error was not apparent, or a serious error threatening the 
integrity of the simulation could not be corrected from the 
GRID, the user still fas the ability to react to the 
Situation. For example, he might undertake procedures to 
help pinpoint the cause of the error, or save the status of 
the simulation for continued processing in the future. AS a 
debugging aid, the Dump command is available to cause a dump 


of the simulation program. 
4.3.2 The Display Screen 


The display screen is a 1024 by 1024 raster unit 
matrix of addressable positions. It is divided by the systen 
into three distinct areas separated by horizontal line 
dividers: 

(1) The Command Input Area - This is a narrow 44 
raster unit area along the very top of the screen. It is the 
area where the alphanumeric command list is displayed. In 
addition, an indicator specifying the state of the GRID is 
shown in the upper right hand ccrner to indicate the 
activity being executed within the system. The indicator may 
take the following values: (1) STATE (blank), which 
indicates that control resides at the GRID and the user can 


input a message of commands for the Monitors’ (zy STATE 2607 


which indicates that information is being transmitted to the 
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360 or that the Monitor is carrying out the actions 
requested. In this state, all user actions at the GRID are 
ignored by the supervisor. (3). -STATE RECEIVE, which 
indicates that the GRID is receiving information from the 
360 (e.g. the display file). In this state, all user actions 
at the GRID are ignored by the supervisor. (4) STATE TOO 
LONG, which indicates that the message being composed by the 
user has exceeded the number of allowable components. The 
user must use the RESET key to delete the message and 
compose another shorter message. 

(2) The Message Area - This area occupies a narrow 
10 raster unit portion of the CRT at the very bottom of the 
screen. It is used for the display of diagnostics frem the 
Monitor and messages resulting from the execution of the 
MESSAGE$ or SUSPENDS functions. 

(3) The Cell Space State Display Area - The 
remaining centre portion of the screen (970 raster units) is 
used to display the state of the space. 

The maximum size of space that can be displayed with 
a regular square geometry is a 40 by 40 cell array. A 
windowing facility is available for large cell spaces to 
enable various portions of the entire space to be displayed. 
If the window is smaller than its maximum size, it is 
centered in the state display area. 

The following Perec rier is provided to describe 


the particulars of the current cell space state display: (1) 
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The title (optionai) is displayed above the centre of the 
state display. (2) The time step number is displayed on the 
right side of the cell state display and in about the middle 
of the screen. (3) The coordinates of the lower left hand 
corner of the display window are displayed immediately below 
the time step number (the X-coordinate first and the Y- 
coordinate beneath). Due to the limited area available to 
display these values, a maximum of six characters can be 
displayed; the least significant digits of coordinate values 
containing more digits than can be accomodated will be used 
(including the sign if negative). However, as a rule, the 
coordinate values will be small enough to be accommodated. 
(4) The transition function entry point number is displayed 
immediately below the window coordinates; and (5) The map 
entry point number is dispiayed immediately below the 
transition function entry point number. This information is 
particularly valuable when hard copy output of the contents 
of the screen is desired. 

An example of the screen configuration is given in 


Figure 4.13. 
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Figure 4.13 Display Screen Configuration 
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4.3.3 The Command Language 


The commands with which the user carries out his 
Simulation are described in the following sections. The 
three character alphabetic command keywords used to 
represent the ccmmands in the command list are placed in 
parentheses following the name of the command. The form of 
all commands which require parameters is illustrated. All 
commands whose descriptions do not evedieiviry indicate that 
parameter input is required can be assumed to have none. 
These commands are effected by entering the command keyword 
and striking the HERE key (represented by a semicolon in the 
command illustrations as well as in the command list 


displayed on the screen). 
4.3.3.1 Transition Calculation Commands 


These commands are used when time stepping the cell 


Space, and managing the progress of the simulation. 


The Step Command (STE) 

This command is used to step the cell space a given 
number of time steps with the current transition function, 
display the resulting cell space state under the current 
map, and return ccntrol to GRID. The Step command has one 


parameter which specifies the number of time steps the cell 
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space is to be stepped before the _ state display is 
calculated and displayed. If omitted, the parameter defaults 
to a value of 1. The Step command has the form: 


STE <integer step value> ; 


The Simulate Command (SIM) 

This command is used to place the simulation ina 
special mode (the simulate mode) used to easily compute 
successive time steps with the current transition function, 
and display them under the current map. The Simulate command 
causes the next time step to he calculated, the cell space 
State to be displayed on the screen, and control to be 
returned to GRID. However, on each subsequent transmission 
from the GRID, after processing the commands in the message 
vector, the Monitor automatically time steps the space once 
and displays it under the current map. In this way, the user 
is able to step through his simulation quickly by just 
hitting the SEND key whenever he is finished examining the 
state display. The simulate mode is terminated when the Halt 
command is processed as one of the commands in the message 
vector. The space will not be time stepped for the 


transmission containing the Halt command. 


The Halt Command (HAL) 


This command is used to terminate the simulate mode 


initiated by the Simulate command. The Halt command restores 
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the system to normal operation so that the Monitor no longer 
automatically time steps the space at the conclusion of 


processing each message received from the GRID. 


The Backup Command (BAC) 

The Backup command is used to restore the state of 
the simulation to that of the previous time step. This is 
possible because of the duplicate data structures of the 
cell space state, one containing the current state, and the 
other the previous state. The Backup command may only be 
performed when there is a valid preceding state data 
structure from which to recover the immediately preceding 
time step. Thus, Backup may not be executed from the initial 
state, immediately after a Restore operation (see section 
Ge3~e3~5,)°the oativatien State Control Commands), or after 
another Backup. The value of the BACKSW flag, the indicator 
which specifies whether backup is possible or not, can be 
determined by means of the Parameter command (see section 


4.3.3.7, Miscellaneous Commands). 
4.3.3.2 Display Control Commands 


These commands are used to allow display of the cell 
Space state under different mappings, and to provide a 


windowing capability for large cell spaces. 
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The Draw Ccmmand (DRA) 

This command is used to recompute the display file 
for the current time step (presumably after changing the nap 
entry point number), and display the cell space state on the 


screen. 


The Size Command (SIZ) 

This command is used to reset the dimensions of the 
display window. The initial setting for the display window 
is the maximum unless the cell space dimensions are smaller, 
whereupon the window is made just big enough to display the 
entire space. The window size can be defined to be any size 
up to 40 by 40. The command has two integer parameters - the 
first specifies the vertical extent of the window, the 
second specifies the horizontal extent. The Size command has 
the form: 


SIZ <vertical extent> <horizontal extent> ; 


The Position Command (POS) 

This command is used for repositioning the display 
window cver desired portions of a cell space which is larger 
than the maximum size of the window. This is accomplished by 
specifying the ccordinates of the lower left hand corner of 
the window. The command has two integer parameters - the 
first specifies the gccooraunace value, and the second the 


Y-coordinate value, of the lower left hand corner of the 
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window. The Position command has the forn: 
POS <X-coordinate> <Y-coordinate> : 


4.3.3.3 Entry Point and Data Specification Commands 


These commands are provided to allow dynamic 
specification of data the Monitor will require during the 
Simulation, including indicating which user defined routines 


are to be used in running the simulation. 


The Set Command (SET) 

This command is used to reassign the user defined 
routines that are used by the Monitor when carrying out the 
various simulation calculations. The form of the command is 
analagous to that of the entry point specification version 
of the SET statement of SLANG. Each of the 8 entry points 
can be reset by assigning the entry point number of the 
desired routine to the entry point in guestion. An entry 
point is denoted by the first three characters of its name. 
The Set command has the following form: 

SET { (<entry point keyword> <entry point number>) }® ; 
An <entry point keyword> can be one the following: EDG, GLB, 
ENE) CMAP? “POS; PRE, “TRA, Or USE; an <entry point number> is 
an integer between 0 and 50 (where "0" indicates the entry 


point is unassigned). 
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The Define Command (DEF) 

This command is used to provide the user with an 
easy way Of making execution time changes to the cell space 
state. The Define command takes a sequence of parameter 
sets, each of which allows the specification of a new 
primitive data type value to a particular field of a series 
of cells. A parameter set consists of three components: (1) 
An integer value specifying the field number of the celi 
data structure which is to be changed. (2) The primitive 
data value which is to be entered into the field specified; 
and (3) A sequence of light pen picks of the cells displayed 
on the screen whose state is to be changed. Each light pen 
pick adds the cell picked to the set of cells which are to 
undergo the indicated change. The primitive data values take 
the following forms (for the indicated types): S=<INT 
constant>, I=<INTEGER constant>, B=<BOOL constant>, K=<REAL 
constant>, and C=<CHAR constant>. This encoding is used to 
allow checking of the type of the constant against the data 
type of the field it is to occupy. If types do not match, 
the Define command is aborted (as are the remainder of the 
commands in the message vector) and an error report is 
given. Each parameter set is enclosed in parentheses. The 
form of the Define ccmmand is: 


DEF { (<field number> <data value> { <light sprck>+})i 2}; 
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The External Input Command (EXT) 

This command is used to define an input stream of 
external cell inputs. There are five such input streams 
available. The External Input command allows specification 
of integer input only (but this input can be used by the 
transition function as a data value specifying in some 
manner a value of any data type). The command has the form: 

EXT <input stream number> <data specification> ; 

The <input stream number> is an integer value in the 
range 1 to 5 specifying the particular input stream being 
defined. The <data specification> is a character string 
which defines a sequence of integer values. The syntax of 


the string is defined as follows: 


(fe fs) 
| eee 
<data specification> --—> { <repeat group> } | | 
| 
pe) 

<repeat group> —> (<repeat count> <repeat group>) 


——> { <integer> } 


R denotes that when the input stream is exhausted it 
will restart at the beginning. This is equivalent to an 
additional layer of parentheses with an infinite count. T 
denotes that when the input is finished it will not restart; 
if a non-circular input stream is specified as the source of 
external input and it has been exhausted, the last integer 
value in the stream will be used. A <repeat group> is a 
string of integer eaites optionally surrounded by 


parentheses with a repetition count <repeat count> following 
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the left parentheses. Nesting of <repeat groups> is 
permitted. With this command, complex input streams of 
integers may be defined. For example, the command 

BRT de (2 (492). (3 259) ):R 
defines input stream 1 to consist of a circular set of 
integer inputs as follows: 


2282 2 Sy Se 8e9o BSL9 22 Cee o 9 78-(9 4g. 9 
4.3.3.4 Macro and User Defined Command Facilities 


The purpose of these commands is to provide the 
command language with flexibility and generality by allowing 
the user to define macros of commands, and to execute his 


Own personal commands. 


The Macro Command (MAC) 

This command is used to define the sequence of 
commands which form the macro. A maximum of five macros can 
be defined. The form of the command is: 

MAC <macro number> "<command list>" ; 
The <macro number> is an integer in the range 1 to 5 which 
is used as an identifier for denoting the macro when it is 
executed. The <command list> is a sequence of commands in 
the same form as they occur in the command list for 
immediate execution upcn ye Raedeaion of the message vector 


to the 360. Commands reguiring light pen input for parameter 
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values are not legal in macros, In addition, macro commands 
cannot be included in a macro definition. The <command list> 
is stored in the 360 and used as the source of command input 
to the Monitor when the macro is executed. 

There are two pseudo-commands which can be used 
Within a macro frogram (but not from the GRID). These are 
the "(" and ")" commands. These pseudo-commands must occur 
in balanced pairs and may be nested to a depth of 5. The 
pair of parentheses specify a substring that is to be 
repeated a given number of times as specified by a 
repetition count which immediately follows the "(", The 
remainder of the string constitutes the sequence of repeated 
characters. For example, the macro 

MAGCV™ ITS" (1G5STE 5 3° SET MAP 22 3 RED s SSET CMAP 21)" = 

defines macro number 1 to consist of the command list 
enclosed in quotes. When By ecete Mane set of commands is 
iterated through 10 times. The parenthesized list of 
commands (less the repeat count) causes: (1) the cell space 
to be time stepped 5 times and the cell space state to be 
displayed under the current map; (2) when control returns 
from GRID, the map entry point to be set to 22 and the space 
state recalculated and displayed on the screen; and (3) when 
contrcl sreturns from GRID, the map entry point to be set to 
21. 

If a command reguiring transfer of control to GRID 


is encountered, execution of the macro resumes when control 
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is passed back to the 360 via the SEND key. However, a check 
is made of the message vector to see if it is empty. If not, 
and if the first command in the vector is an AMP command, 
the macro program is terminated, and the remainder of the 


message vector is processed. 


The Execute Macro Program Command (EMP) 

The Execute Macro Program command is used to 
initiate execution of a command macro. It has a single 
parameter, an integer in the range 1 to 5 specifying the 
particular macro to be executed. A macro program terminates 
when (1) the macro is completed, (2) the user issues an AMP 
command, or (3) an error condition arises. The form of the 
command is: 


EMP <macro number> ; 


The Abort Macro Program Command (AMP) 

This command is used to terminate execution of a 
macro program. If control is passed to the GRID as a result 
of normal processing of a command in the macro (e.g. to 
display the current state of the cell space), the user can 
terminate execution of the macro by entering the Abort Macro 
Program command at the beginning of the command list. It 
must be the first command in the message vector in order to 


have any effect, otherwise the Monitor ignores the message 


vector and continues processing the macro. 
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The User Command (USE) 

The User command is used to provide the user with a 
means of executing his personal commands by initiating a 
cali to the current user routine. A character string 
parameter list enclosed in quotes can accompany the call; it 
is decoded and made available to the user routine through 
the variaktlies USERARAY, USERINT, USERINTR, and USERREAL. The 
description of these variables and the form of the parameter 
list is described in section 4.2.2.1, the User Command 
Variables. The form of the command is: 

USE { "<character input string>" }!1 ; 

One of the important aspects of the user command 
facility is that it gives the user a means by which he can 
change the values of system variables, as well as personal 
data, at run time. (In this respect, special consideration 
has been given to easily changing the cell state data 


structures by means of the Define command.) 
4.3.3.5 Simulation State Control Commands 


These commands are concerned with providing the 
ability tc restart a simulation at the beginning or from an 
arbitrary intermediate time step, and to allow a simulation 


to be executed at various sessions at the GRID. 
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The Clear Command (CLE) 

This command is used to reset the simulation to its 
initial conditions in order to allow the simulation to be 
restarted, At the beginning of the simulation, a sequential 
file called RESTART will be created and the crucial 
information about the status of the simulation will be 
placed in this file. This information consists of pointers, 
flags, and system data values located in a contiguous block 
in KOMMON. The data is written as one lcng record into 
RESTART. 

When the Clear command is executed, the contents of 
the RESTART file are read into the appropriate locations in 
KOMMON, thus resetting all pertinent parameters to their 
initial value. The user specified entry points are not 
included in this bliock of information so that they can be 
changed as desired before the simulation is restarted. This 
allows the user to employ different initialization 
procedures and transition functions with which to rerun his 
Simulation. 

Once the system has been reset, control is passed to 
the global initialization routine specified by the GLBINIT 
entry point, thereby restarting the simulation at the 


beginning. 
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The Save Command (SAV) 

This command is used to save the state of the cell 
space and status of the system at an arbitrary time step. 
The entire set. of pointers, flags, and data values in 
KOMMON, along with the cell state data structures are stored 
in a sequential file. If the file already exists, it will be 
emptied f£inst; !ethervise!’ at “will be’ created:  Tf£° “the 
specified file exists but is not sequential (and therefore 
will nct accept the large length records which are to be 
placed in the file), the command will be aborted and _ the 
reason reported to the user. The largest record length 
possible is used for fast and efficient data transfer. The 
purpose of this command is to allow the user to stop his 
Simulation and resume it at another time, and to allow 
recovery of a given cell space state. The command takes one 
parameter, the name of the file in which the cell space 
state and system status is to be stored. The form of the 
command is: 


SAV <filename> ; 


The Restore Command (RES) 

This command is used to restore the state of the 
cell space and status of the system to an arbitrary time 
step from data saved in a file by the Save command. This 
command can be used ge restart a Simulation from an 


arbitrary time step. All the data necessary to completely 
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restore the simulation is preserved by the Save command. If 
the simulation is stopped and restarted at another time by 
use of the RESTORE command, the object modules containing 
the mcdel description and user defined routines should be 
the same as those used for the initial simulation from which 
the stored file was generated. The command takes one 
parameter, the name of the file from which the cell space 
state and system status is to be read. The form of the 
command is: 


RES <filename> ; 
4.3.3.6 SUSPEND$ Function Commands 


The SUSPENDS function of SLANG is provided to allow 
the user to gain control at GRID during execution of a user 
defined routine. This provides the user with a means of 
error recovery and run time decision making. When the 
SUSPEND$ function is executed, the point of execution is 
saved along with the remaining portion of the message vector 
(Or macro) and control is passed to GRID. 

The simulation operations in which user defined 
routines take part include the initialization procedures, 
display calculations, user commands, and transition 
calculations. The most crucial suspensions occur when the 
Space is being time efouved and an unexpected situation has 


been discovered. The principal factor controlling these 
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Situations is the state of the cell being time stepped or 
the state of one of its neighbors. The user has access to 
the cell state data structure ($0LD) via the Define command, 
and thus can recover from errors by changing the state of an 
offending cell without endangering the integrity of the 
remainder of the simulation. 

When contrcl is passed to GRID, the user has the 
full srescurces of the Simulation system at his disposal to 
examine the conditions leading to the situation and to 
handle them in the appropriate manner. The commands for time 
stepping (or backing up) the simulation must not be used 
until the suspension has been cleared. 

The user has the following two commands available to 


clear the suspension and resume the simulation. 


The Return Command (RET) 

This command is used to remove the suspension and 
resume execution where it left off. The remainder of the 
command being executed at suspension is completed, and _ the 
remainder of the ccmmands in the message vector (or macro) 
at suspension are executed. For example, the Return command 
can be used when a suspension is made to make a decision 
(say by setting up some variables by means of a user 
routine); once the variables are set, execution should 
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proceed from where the user routine was interrupted. 
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The Cancel Command (CAN) 

This command is used to remove the suspension, and 
abort execution of the user defined routine, the command in 
progress at suspension, and the remainder of the commands in 
the message vector (or macro) at suspension. The commands 
which follow the Cancel command in the message vector are 
then executed. For example, the Cancel command can be _ used 
to clear a suspension when an undefined transition is 
discovered and the state of the cell at fault is changed (by 
the Define command); the transition should be aborted and 


the time step tried again with the new cell state. 


4,3.3.7 Miscellaneous Commands 


These commands are used for various unrelated 
actions which might be required during the course of a 


Simulation. 


The Parameter Command (PAR) 

This command is used to display the values of 
pertinent system variables which determine the current 
status of the simulation and system. The system variables 
appearing in the list are: 

NTRANS, the current timestep. 
XDIM, the horizontal extent of the space. 


YDIM, the vertical extent of the space. 
XDIS, the horizontal extent of the display window. 
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YDIS, the vertical extent of the display window. 

XLLC, the X-coordinate of the lower left hand corner of 
the window. 

YLLC, the Y-cocrdinate of the lower left hand corner of 
the window. 

NUMCELLS, the number of cells in the space. 

NUMNEIGH, the number of cell neighbors. 

BACKSW, the backup flag indicating if backup is possible. 

EDGENT, the current edge routine entry point number, 

GLBINIT, the current global initialization routine entry 
point number. 

INITRANS, the current initialization routine entry point 
number. 

MAP, the current mapping routine entry point number. 

POSTRANS, the current post-transition entry point number. 

PRETKANS, the current pre-transition entry point number. 

TRANS, the current transition routine entry point number. 

USERCOM, the current user routine entry point number. 

SUSEFEND, an indicator specifying if there is a suspension 
in effect. 

XSUS, the X-coordinate value of the cell being time 
stepped (meaningful only if a suspension is in effect 
which occurred during a time stepping operation). 

YSUS, the Y-cocrdinate value of the cell being time 
stepped (meaningful only if a suspension is in effect 
which occurred during a time stepping operation). 


In addition, the full set of command keywords are 
listed (without explanation) to provide the user with a 
quick reference. All data is displayed in the cell space 
state display eee of the screen. The list of command 
keywords is displayed along the left edge, and the systen 
variakles are displayed in the centre. Figure 4.14 
illustrates the display configuration generated by the 


Parameter command. 
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Command Keywords 


AMF 
BAC 
CAN 
CLE 
DEF 
DRA 
DUM 
EME 
EXI 
EXT 
HAL 
MAC 


eee a a ee ee Se ee 


PAR 
PLO 
POS 
RES 
RET 
SAV 
SET 
Sih 
SIZ 
STE 
TTL 
USE 


Figure 4.14 
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System Variables 
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Parameter Command Screen Configuration 
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The Plot Command (PLO) 

This command is used to obtain hard copy output of 
the picture displayed on the GRID screen. The Plot command 
causes the display file to be output to a file. A program is 
available [28] to transform the file into input compatible 
with a model 770/663 CalComp Plotter which is available at 
the University of Alberta. The command has one parameter 
specifying the size of the display screen output. The form 
of the command is: 

PLO <integer scale factor> ; 
The plot is a square picture which depicts the contents of 
the display file as they appear on the screen (this does not 
include the command list which is generated by the 
supervisor and therefore not contained in the display file). 
The <integer scale factor> is an integer in the range 5 to 
28 which specifies the length of the edge of the square 
plotted in inches. A scale factor of 10 is approximately 
equivalent to the actual screen size, while the maximum size 
that can be accomodated by the plotter is a 28 inch square. 
The Plot command places all display file data into a file it 
creates called PLOTFILE. As each plot is requested, it will 


be added to the file. 
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The Dump Command (DUM) 

This command is used to cause a core dump of the 
user's simulation program to aid in debugging when a machine 
level error is encountered during execution of the 
Simulation or _for whatever other reason such data is 
desired. The Dump command places all output dump information 
into a file it creates called DUMPFILE. As each dump is 


requested, it will be added to the file. 


The Title Ccmmand (TTL) 

This command is used to specify a character string 
Which will appear centered above the display of the cell 
space state. It is used solely for documentation purposes. 
The command has one parameter, a character string of up to 
79 elements enclosed in quotes. It has the form: 


TIL "<character string>" ; 


The Exit Command (FXI) 

This command is used to end the simulation. Upon 
encountering this command, the Monitor terminates execution 
of the simulation and indicates this on the 2741 by printing 
the message: 


EXECUTION TERMINATED 


A list of the commands in the command language is 


given in Figure 4.15. 
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The GRID Commands 


AMP 
BAC 
CAN 
CLE 
DEF 
DRA 
DUM 
EMP 
EXI 
EXT 
HAL 
MAC <macro number> "<command list>" ; 

PAR 3; 

PLO <integer scale factor> ; 

POS <X-coordinate> <Y-coordinate> ; 

RES <filename> ; 

RET 3 

SAV <filename> ; 

SET { (<entry point keyword> <entry point number>) }® ; 
SIM 3 

SIZ <vertical extent> <horizontal extent> ; 

STE <integer step value> ; 

TTL "<character string>" ; 

USE { "<character input string>" }!1 ; 


(<field number> <data value> { <light pick> }) } ; 


macre number> ; 


input stream number> <data specification> ; 


wo Ae A oe we pu we 2 ws oe 


Special Macro Pseudo-Commands 


( and ) 


Figure 4.15 The Command Language 
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CHAPTER V 


EVALUATION AND DISCUSSION 


5.1 Evaluation of the Proposed Package 

The task of this chapter is to attempt to determine 
in some general ways, "how the system measures up". The 
terms of reference set forth earlier give rise to two 
approaches to this_ problem: (1) how well does the systen 
fulfil the basic requirements of an interactive simulation 
system as outlined in Chapter II; and (2) how does the 
system compare to its predecessor, R. F. Brender's cellular 
Simulation system (Section 3.4.1) ? 

The proposed package uses the powerful simulation 
language/simulation system approach to simulation. The 
Simulation language has Une only a basic programming 
Capability but also specialized Capabilities designed 
strictly for interactive cellular space modelling. Included 
in the language is an essential core of programming 
facilities which provide a capability powerful enough to 
build the types of programs for which the system is 


designed. For instance, the language has an especially 
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powerful data structuring facility. No attempt has been made 
to provide an exhaustive and elaborate set of facilities 
with the notion of creating truly general purpose 
programming capabilities. The frills of most general purpose 
languages have been omitted leaving only a basic programming 
capability with which a novice programmer can describe the 
behavior cf his model. 

Of primary importance however, is that the language 
provides special means by which the user may easily and 
efficiently specify his model description. Special 
statements are included to handle the description of each 
physical characteristic of cellular spaces. General methods 
of specification allow the user to define the details of 
these aspects of the model himself or to choose system 
facilities which handle the specification of some common 
cellular model characteristics or problems (e.g. the 
external cell wrapping routines). 

At run time, the simulation supervisor embedded in 
the Meniter provides the execution time management routines 
needed for a simulation. The supervisor maintains stepwise 
execution of the simulation as directed by the user. For 
example, the tasks of cell space data structure creation and 
referencing, timing management, state display management, 
transition management, and exreenal cell handling have all 
been provided for within the simulation system. Thus, while 


only a basic programming capability has been provided for in 


_ bails tt a mf 


onan wis we shisbhc ai peas ee 
solsitivat %6 Bak OUNHoMsls Sus ‘ovbstapedte” 
ssorusg fesarep yams twitesw> %o ' wofton 
yeoQg Tn ‘fsxeah" ror 40 effin? off wmetvilideqs2 po 
tebnw erp bead. = ¥Lbd pwivest Seetéas wand evad re: 
edt saiasb aeewededs yoo, «etven & Wakde dste vibitdeqes: 
| tebon aid to sobvaded 
aneveest eft feds BL .ievevnd wonetroqet yYaair4 to or ¥ 
Soo _¢itaee Yew sees 213 fodde qi easpaoLaioege esbivord 
fs raand shoo? tase vb fsbo 7 oe (iisege  yltoetortis 
dyer 36 Wettalence’ edd elbaad os febelanl ots 2fis¢esssa 
pteetton EPiiewcesonjs selulive “to Divelmeso81445 inoteyag 
t) difteteh G49 «sitet 69 oer ode wollte hobsootiioeqa: 30 
wadeys aeoois ne to tleestd ishoe edt 90 asooqad ‘enoas 7 
Weydoo Meshes fo solrs-iiteeqs 40s ofSapd Motew sekeRitegs | 
4) Geie)> G#eiosq x0 nnd sa Fees peneen ave <ehniteeG 
| .(sealtvog pnégqs1y leo Lansetee, 
CP BONN HEL resco svt ruluede ade gomLe aor ah bie Sal 
‘ dbie200d Saeebereee ani+ wobtucwes 949029076 —— 
oelvysta, athbeetes soviersqes ott, mokghtomte-n ® aT) 
Wet Peer Geyseutt a» Hols eiodde bode <2 prion 7 
“ae nooo ata 903e ath conga flue" eo naest oto Gory 
afesaspeate yakgath eters tesa genem parnet 
ar tbaad 2 p Lawrene chegeaete ae 


184 


the language, specially tailored facilities for cellular 
Space modelling are available in the system as a whole. 

For the types of models for which the package is 
designed, the proposed System provides adeguate and 
generally substantial facilities in each area of capability 
displayed by general simulation languages (as listed in 
Chapter II). The power of the MIS system under which the 
Simulation is run adds greatly to the facilities available 
to the programmer. For example, access has been provided to 
the powerful file handling facilities and the broad range of 
statistical and mathematical packages available (via the I/0 
and external program statements in the language). 

In addition, the system has numerous features which 
take advantage of the interactive simulation environment. 
The simulation language provides many facilities to allow 
the user to access information and data relevant to his 
Simulation. The Iyv0O routines allow data transfer to or from 
any file or device (other than GRID), with special 
facilities available in the simulation system to transmit 
certain types of data (such as character strings or the cell 
space data structures) for display on the GRID. The SUSPEND$ 
facilities provide a run time communications link with the 
user from within his frogran. 

The design of the package has also taken into 
account basic ere te ea bare like efficiency and 


user aids. The entire package is being written in 360 
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assembler language, and the compiler produces its own object 
modules instead cf requiring an assembly step to do this. 
Comprehensive error detection procedures are being built 
into the compiler, with detailed error diagnostics supplied 
to the user. The compiler will provide for special 
alphabetic listing cf the entries in the symbol table. Cross 
referencing tables and maps are provided at loading time. 
All these types of facilities combine to produce an 
efficient, user-oriented systen. 

The simulation system makes use of a very large and 
powerful general purpose computer and a functional graphical 
display computer connected by a very high speed interface. 
The value and versatility of the graphics aspect of the GRID 
have been utilized in the system. The run time facilities 
with which the user carries out his simulation are both 
powerful and easy to use. A concise command language and 
versatile command entry facility are used to direct the 
Simulation. Command editing features, a macro command 
ability, and user defined command facilities have all been 
provided. This last facility in particular is valuable in 
that it introduces a great deal of flexibility into the 
Simulation system. The user can use these personal commands 
for error detection, output or specification of data, and a 
host of other purposes. The capabilities available via the 
command language provide the user with complete control over 


the execution of the simulation and with monitoring 
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facilities to view its progress. The system allows the user 
to restart his Simulation, or save (and thus later restore) 
the status of his simulation at any point in simulated time. 
The system design is functional in that the user defines the 
routines which Carry out the crucial aspects of the 
Simulation (as well as optional routines for special 
purposes), while the simulation executive handles’ the 
Management end of the simulation. These routines can be 
respecified at run time to allow flexibility within the 
system and give the user control at ali levels while the 
Simulation is under way. Special routines can be specified 
to do initialization procedures, statistics gathering, 
inputyoutput, and data analysis. Run time debugging features 
have been incorporated into the system to give the user the 
ability tc detect errors as they occur, save the status of 
the simulation when it is threatened, and the ability to 
help pinpoint and recover from errors. 

The system description should have illustrated the 
utility, versatility, flexibility, and power of the proposed 
package for cellular space modelling. M. R. Lackner has 
written: "In conversational mode, simulation - the deduction 
of a sequence of system state descriptions - must be 
interruptible, restartable, modifiable, and above all 
displayable" [32]. An attempt has been made to embody each 


of these gualities in the system design. 
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Consider now the second yardstick against which the 
proposed system can be measured. Brender's system has proven 
to be very valuable for cellular space modelling; the value 
of the proposed system may be partly judged by how it 
compares with its predecessor. 

Much of the top level design regarding the 
appearance and programming facilities provided within the 
package has been adapted from Brender's system description. 
Although the methods of implementation were developed 
independently, both systems have similar objectives and deal 
with the same types of models. Thus the capabilities and 
general framework are Similac. There are, however, 
Significant differences in the two systems. 

It is apparent that differences in hardware are 
responsible to some degree for differences in both the power 
of the systems and the methods of providing some types of 
facilities. For example, the command language and display 
facilities are approached differentiy. Brender's system 
makes greater use of the display screen and light pen 
facilities. Most commands are entered by means of a light 
pen and executed immediately upon entry (even for entering 
arithmetic data, numeric characters are displayed on the 
screen and light penned in sequence rather than being 
entered from the keyboard). This approach is made possible 
from a response point of view due to the fact that the user 


has dedicated hardware at his disposal. In the time-sharing 
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environment of the proposed system however, the GRID must 
compete with other remote terminal (and batch) users for the 
360's resources. Response time degrades during periods of 
peak usage. This is partially responsible for the approach 
taken in the _ proposed system, which uses a much more 
keyboard oriented command input facility (since each display 
change requires transmission to and from the 360) and stacks 
sequences of commands before transmission for execution. 
Additionally, it would seem that many types of data input 
are as easy and natural from the keyboard as with the light 
pen. 

Many differences between the systems are of 
borderline importance. For example, the proposed system has 
a wider range of data primitives with a greater value range, 
a Simpler operator/function priority system, a more 
efficient compilation procedure, but a slightly more rigid 
program structure. Other differences deserve more attention. 
The interactive communications capabilities provided by the 
SUSPEND$ and MESSAGE$ functions are much more powerful and 
useful than the interactive facilities provided in Brender's 
language. Likewise, the 1/0 facilities and data management 
capabilities of the proposed system are better. The 
debugging facilities are more sophisticated and user- 
oriented. The manner of producing hard copy output of the 
display screen is au eres in)’ that a scaled. plot, is 


available whereas a picture must be taken by camera in 
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Brender's systen. 

A look at the commands in the command language 
reveals that although many are quite similar, those in the 
proposed system are more general and versatile. For example, 
the commands for time stepping the cell space allow 
calculaticn of more than one step between display images if 
desired. The list of system variable values available for 
display on the screen is larger, giving the user immediate 
access to more information. 

Although in many ways the two systems are very much 
alike, the proposed system surpasses Brender's system in two 
important aspects: generality and flexibility. These 
qualities have been explicitly emphasized in designing the 
proposed system. A particularly obvious example is that 
(other than mapping functions) only one of each type of user 
defined routine is allowed in Brender's system; the proposed 
system, however, allows multiple definition of all user 
defined rcutines and execution time respecification of the 
particular ones used in the simulation calculations. A total 
of 50 such user defined routines are allowed (the number 
could easily be expanded). Since there is no other limit on 
the number of each type of routine, the user has 
considerable flexibility in meeting the needs of his 
simulation. Another example of the increased generality is 
that the maximum size of the cell space is considerably 


larger at 128 by 128 cells compared to the 25 by 25 maximun 
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possible on Brender's system. Much larger data bases can be 
accommodated in the proposed system. 

In conclusion, the system proposed in this study 
Should be at least as valuable a research tool in cellular 


Space modelling as Brender's system. 
5.2 Status of the Implementation 


implementation of the proposed system has. been 
started. The top level design has been completed and many of 
the system facilities are in varying stages of development. 
Much high level fiowcharting work regarding the simulation 
calculaticn procedures has been completed. Algorithms have 
also been developed for such routines as the assignment 
statement parser. The framework of the compiler has_ been 
programmed with respect to such procedures as_ source 
statement type analysis, program listing generation, and 
error reporting procedures. The routines handling symbol 
table entry and referencing, the data structure facility, 
and cell space data structure referencing have been largely 
written. 

The program structure of both the compiler and 
Monitor is modular in design so that the implementation of 
the total system can proceed in a stepwise manner. 
Individual facilities aan then be programmed and fitted into 


the overall system structure with a minimum amount of 
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difficulty. 
even simple 
first. The easiest way to 
implemented 


the package. These facilities 


user greater control over 
auxiliary aids. 
The programmer can 


deletions from the simulation 


(1) The formatted 
Since the default formats are 

(2) The DATA 
assignment statement used in 


do an equivalent job. 


A basic core of facilities is required to 


statement 
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allow 


applications and therefore should be implemented 


delimit the core which must be 


first is tc describe the nonessential aspects of 


generally involve giving the 


his simulation and valuable 


manage with the following 


language: 
output facility is not needed 
unsophisticated but workable. 

be and the 


can omitted 


the initialization routines to 


(3) The external function call facility is not 
vital; its omission forces the user to define all the 
programs he requires himself. 

(4) The BASEORG statement is valuable for large 


routines, but can te replaced 


sections which can be 


sequence. This statement, 


(5) The external cell routines 


by modularizing a routine into 


defined as functions and called in 


however, is trivial to implement. 


XWRAP, YWRAP, and 


TORUS could be left for the user to progran. 


(6) The 


monadic operators for bitwise complement, 


exponential, natural log, square root, and the trigonometric 


functions are not essential. 


Likewise, the dyadic 


operators 
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for modulo, bitwise logical operations, and shifts are not 
necessary. 

The user can manage with the following deletions 
from the command language: 

(1) The Simulate command facility could be omitted 
Without preventing control over time stepping the 
Simulation. However, this might make running the simulation 
very tedious. 

(2) The windowing capability could be temporarily 
omitted and thus have the restriction of a maximum size of 
cell space at 40 by 40 celis. 

(3) The Define command could be left out with the 
result that the user could not easiiy change the cell space 
data structures from the GRID, but would have to employ a 
user defined ccmmand to do so. 

(4) The macro command facilities could be omitted 
entirely with very little effect. These were included mainly 
aS a convenience and for generality. 

(5) The Title, Plot, and Dump commands can be left 
out with little effect on the actual simulation, but with a 


loss of valuable user aids. 


Other than these two sets of facilities, the 
simulation package should remain intact to provide enough 
programming and run time monitoring capability to actually 


describe a model and run a Simulation. As a final drastic 
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measure, the user defined command facility could be omitted, 
but this would mean a great loss of power and control for 
the user. 

The implementation should take into account the 
omitted facilities in order to make it an easy task to add 
them later. For example, ali the omitted statements and 
commands should be valid and undergo preliminary processing. 
They should be analyzed to the point where control is passed 
to the routine which is to handle execution of the actions 
required by that statement or command. This routine, 
however, could be temporarily replaced by another which 
would generate a message to indicate the facility was not 
yet available. In this way, when a facility is implemented, 
it can more easily be fitted into the system without the 
need f0r major reprogramming, since a place has been 


reserved for it. 


In summary, this research has examined two areas, 
interactive simulation and cellular space modelling. An 
attempt has been made to create a system specifically 
designed for describing and interactively running cellular 
models. A detailed proposal for a cellular space simulation 
system of this sort has been presented, and an 
implementation started. In addition, a set of guidelines has 
been developed to aes in the design of interactive 


Simulation systems. 
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